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LIFETIMES OF FAST IONS IN THE PLASMA OF DCX 


C. F. Barnett, J. L. Dunlap, R. S. Edwards, G. R. Haste, H. Postma, 
J. A. Ray, R. G. Reinhardt, W. J. Schill, R. M. Warner, and E. R. Wells 
Oak Ridge National Laboratory, * Oak Ridge, Tennessee 
(Received May 4, 1961) 


Previous publications’* have outlined the gen- 
eral construction of the DCX machine and have 
described various measurements of the plasma 
ring formed in DCX by injecting 600-kev diatomic 
hydrogen ions in the median plane and dissociating 
them with a vacuum carbon arc. A plasma ring of 
300-kev protons concentric with the magnetic axis 
can also be formed without an arc, since dissoci- 
ating collisions occur between the diatomic ions 
and residual gas molecules. In this Letter we re- 
port the results of recent lifetime studies for pro- 
tons trapped in DCX via the latter mechanism, gas 
dissociation. These studies show mean containment 
times in the order of seconds at initial proton den- 
sities of 10’ cm™°. 

The present experiments were conducted in the 
DCX 2:1 mirror geometry with an axial midplane 
magnetic field of 10 kilogauss and an H,* beam of 
about 1.4 ma. With the H,* beam entering the 
system, oil diffusion pumps and an evaporated 
titanium film produced a pressure of 5x10™° mm 
Hg within the water-cooled copper liner which 
surrounded the plasma region. The steady-state 
fast proton density was about 10’ cm™~°, as de- 
termined from the expression 


n =In 0 ,L7/V, 


0 


where J is the injected H,* current, m, is the 
neutral gas density, og is the dissociation cross 
section for proton production, L is the path length 
of the H,* for dissociations that result in effective 





trapping (approximately 10 cm), 7 is the measured 
lifetime of the fast protons, and V is the plasma 
volume (about 10*° cm‘). 

A major loss mechanism for the circulating 
protons is electron-capture collisions with the 
background gas. The lifetimes of the protons 
were deduced from measurements of the ener- 
getic neutrals resulting from these collisions. 
This neutral flux is a measure of the charge- 
exchange reaction rate within the plasma volume, 
N.N90-,v, where o,, is the electron capture cross 
section and v is the velocity of the circulating 
protons. The rate at which this neutral flux de- 
cays after the H,* beam is interrupted is then 
the rate at which the fast proton density disap- 
pears, since ”90,,v is nearly constant during 
the decay interval. If charge exchange is the 
only loss process, the proton density and the 
neutral flux should decrease exponentially with 
a time constant, tT =1/("g0¢,v). 

In these experiments a portion of the energetic 
neutrals was passed through a 10~°-inch thick 
nickel foil and then monitored by a Faraday cup. 
This foil thickness is sufficient to insure a sta- 
tistical charge state equilibrium of the ions and 
neutrals escaping the back surface. For energies 
of 300 kev the emerging beam is greater than 
ninety-nine percent protons. Particle beams 
with energies greater than 150 kev are not at- 
tenuated by this thickness. The current from 
the Faraday cup was amplified by a transistor 
circuit and then displayed on an oscilloscope. 
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FIG, 1. Reciprocal lifetime as a function of hydrogen, 
helium, and residual gas pressure. 


Photographs of the oscilloscope trace decay after 
interruption of the H,* beam provided curves of 
the decay of the fast proton density. 

The pressure responses of the decay traces to 
helium and hydrogen were determined in the 
range from base pressure to 1x107° mm Hg. 

The decay curves were exponential through this 
pressure range. Figure 1 shows the reciprocal 
lifetimes observed during the helium and hydro- 
gen runs, along with the reciprocal lifetimes 
expected from the loss of 300-kev protons by 
charge-exchange collisions in each gas. The 
deviation of the hydrogen curve from a 45-degree 
slope results from the presence of background 
gas contaminants which possess larger électron 
capture cross sections. The cross sections of 
helium and the contaminants seem to be about 
equal, and the helium curve has a 45-degree 
slope. For helium the measured lifetimes are 
80% of the calculated. For hydrogen at the higher 
pressures, where the effects of the contaminants 
are the smallest, the measured lifetimes are 60% 
of the calculated. The discrepancy is attributed 
to small errors in measurements and to the 
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FIG. 2. Amplitude of neutral-particle detector 
signal as a function of time. 


energy lost (approximately 10 kev) by the protons 
through ionization and excitation collisions with 
the gas. 

The residual gas curve of Fig. 1 was obtained 
by cooling the liner which surrounded the plasma 
region with liquid nitrogen and photographing the 
decay traces as the pressure decreased. Figure? 
shows a typical decay trace obtained at the lowest 
pressure. The semilog plot shows an exponential 
decay with a time constant of 3 seconds. That the 
decay traces are exponentials constitutes evidence 
for a single loss mechanism. 

Linear plots of the measured reciprocal life- 
times as functions of pressure show that the 
helium and residual gas data extrapolate to in- 
finite lifetime at zero pressure, as is expected 
if proton loss is due solely to charge exchange. 
As a result of the higher average electron cap- 
ture cross section of the background gas contam- 
inants, the hydrogen data extrapolate to a finite 
lifetime. 
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We believe that the present experiment is sig- 
nificant in that we have shown that a mirror 
machine is capable of containing a plasma for 
seconds at densities of 10’ cm™* with charge 
exchange as the only apparent mechanism for 
particle loss. Our future work will be under- 
taken to increase the density of stored ions, to 
study the extent of orbital disorganization, and 
to investigate the homogeneity of the ion distri- 
bution. 


*Operated by Union Carbide Corporation for the U. S. 
Atomic Energy Commission. 

‘C. F. Barnett et al., Proceedings of the Second 
United Nations International Conference on the Peace- 
ful Uses of Atomic Energy, Geneva, 1958 (United 
Nations, Geneva, 1958), Vol. 31, p. 298. 

2a, H. Snell, Fourth International Conference on 
Ionization Phenomena in Gases, Uppsala, Sweden, 
1959 (North-Holland Publishing Company, Amster- 
dam, 1960), Vol. I, pp. 997-1001. 

*J. L. Dunlap, IRE Trans. on Nuclear Sci. NS-7, 
No. 4, 19 (1960). 




















EVIDENCE ON THE LAMINAR NATURE OF THE EXOSPHERE OBTAINED BY MEANS OF GUIDED 
HIGH- FREQUENCY WAVE PROPAGATION 


R. M. Gallet and W. F. Utlaut 
National Bureau of Standards, Boulder, Colorado 
(Received May 1, 1961) 


The purpose of this Letter is to report the re- 
sults of a series of experiments which take ad- 
vantage of the laminar structure of the exosphere 
for testing the existence and properties of the 
sheets at large distances above the earth. In ad- 
dition, it is suggested that this new mode of prop- 
agation may be useful for observing, in a quasi- 
continuous manner from ground observations, 
the characteristics of the permanent ring cur- 
rent, recently observed in Explorer VI and 
Pioneer V measurements.' 

It has been known for some years? that the 
structure of the earth’s exosphere seems to be 
somewhat fibrous, in the sense that the elec- 
tron density distribution has slight increases 
distributed in thin sheets or columns along the 
magnetic field lines. This property is evidenced 
mainly by the observations made on “whistlers.”® 
These are a very specific type of natural radio 
signals of very low frequencies propagating a- 
long the curved magnetic lines of force and 
generated by lightning discharges. The obser- 
vations show that quite often the signals are 
guided along several, well-separated, discrete 
paths. Also, the existence of long trains of 
echoes, in which the signal is reflected back 
and forth many times between the two ends of 
the line of force with a remarkably small total 
attenuation, requires that the waves be guided 
very efficiently within such columns or sheets, 
in a way similar to the guided propagation of 
light in long glass fibers (“light pipes”) or radio 
waves in dielectric waveguides. In the whistler 






mode the frequency of the signal has to be small- 
er than the electron gyrofrequency, wy =eH /mc, 
and smaller than the angular plasma frequency, 
wp, given by wp?=4mNge?/m. Under such condi- 
tions, the refractive index is usually very large 
(from 7 to 30 within the exosphere) and there are 
strong refraction effects, complicated by a great 
anisotropy and a large dispersion. 

The energetic particles comprising the outer 
Van Allen belt have also been observed by Pio- 
neer IV* and by Explorer VII,° the latter under 
disturbed conditions, to be distributed in thin 
layers which are remarkably sharp. The re- 
lationship between the density distribution of 
the high-energy electrons of the Van Allen radi- 
ation, which have a low space density, and the 
more numerous thermal electrons forming the 
bulk of the exosphere is, however, not yet de- 
termined. 

The present series of experiments was stimu- 
lated by the suggestion made by Obayashi® of the 
possibility of high-frequency wave guidance 
through the exosphere. In this case, the fre- 
quency is well above the plasma frequency and 
gyrofrequency, and the refractive index is al- 
most equal to unity; consequently, the wave ve- 
locity is very nearly that for free space. We 
believe, however, that guidance along the line 
of force is afforded by a relatively small gradi- 
ent of electron density transverse to the field 
under circumstances where the gradient is suf- 
ficient to produce the slight bending equal to the 
curvature of the line of force. Also it is neces- 
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sary that the wave enter the exosphere in a direc- 
tion tangential to the magnetic field. 

It is interesting to note, historically, that 
Pedersen in a 1929 paper’ suggested the possi- 
bility of propagation along a band of electrons or 
ions which connects the northern and southern 
polar areas and he considered such bands to be 
formed by the magnetic field of the earth. This 
mode of propagation was suggested by Pedersen 
to provide an explanation of the long-delay ech- 
oes observed at Oslo, Norway, in 1927 by Hals® 
and later by Van der Pol and Stgrmer at a fre- 
quency of 9.56 Mc/sec.® A later search for 
such echoes, on frequencies of 13.455 Mc/sec 
(30 kw) and 20.675 Mc/sec (9 kw), was made in 
the period 1947-1949 by Budden and Yates,’° 
and during this period echoes considered to be 
of long delay were not observed. Observations, 
however, were made at several selected inter- 
vals of an hour or two in length one day each 
fortnight only and, if conditions permitting such 
echoes exist for short periods, a small percent- 
age of time, the relatively infrequent periods of 
observation may explain why long-range echoes 
were not detected. 

The experimental work of Obayashi was per- 
formed in Japan at a geomagnetic latitude of 
26°. The magnetic line originating there is a- 
bout 7000 km in length and, at most, reaches a 
height above the magnetic equator of 1500 km. 
Thus, most of the propagation path lay below 
the exosphere in the normal ionosphere. In 
order to test better the properties of the exo- 
sphere, experiments conducted at a much higher 
geomagnetic latitude are necessary so that the 
radio wave may traverse sufficiently far into the 
exosphere. Fortunately, a transmitter located 
near Washington, D. C. (geomagnetic latitude 
50.2°) was available for use in the present series 
of observations. The line of force beginning at 
this point, computed using spherical harmonic 
analysis of the geomagnetic field with 48 terms 
included (see Vestine and Sibley"), has a maxi- 
mum height above the earth’s surface of 10 460 
km (2.64 earth radii measured from the center 
of the earth) and a length of 33 650 km to the cor- 
responding magnetically conjugate point in the 
South Pacific. 

Attempts were made on various occasions 
during a two-month period from late April to 
June of 1960 to observe the signal propagating 
via the suggested mode and backscattering from 
the ground at the point in the other hemisphere 
which is the magnetic conjugate of the trans- 
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mitter. A high-frequency radar having a peak- 
pulse power output of approximately 100 kw was 
used. To aid in signal identification a special 
pulse group of two one-millisecond pulses sepa- 
rated in time by 8 milliseconds was transmitted 
at a group repetition rate of 2.5 cps. A super- 
heterodyne receiver, with a bandwidth of approxi- 
mately 1 kc/sec, was used in conjunction with a 
low-noise preamplifier to receive the signal of 
which range-time and range-amplitude presenta- 
tions were photographically recorded. The an- 
tenna used for both transmitting and receiving 
was an array of two vertically polarized 3-ele- 
ment Yagis which had been designed to operate 
on a frequency of 13.7 Mc/sec. This was tipped 
upward 71° to align it with the earth’s magnetic 
field. Observations were made only during the 
nighttime, because of the limited frequency 
range of operation imposed by the antenna and 
because of the marginal sensitivity of the radar, 
so that the benefits of low critical frequency and 
low absorption of the ionosphere could be utilized. 
Table I summarizes the periods of observations 
and gives the radar range and times during which 
a clearly recognized signal was detected. 
Backscatter echoes were found to exist for a 
length of time ranging from a few minutes up to 
a period in excess of 50 minutes. The received 
signals were weak, and it is estimated that the 
transmission loss for these signals exceeded 
200 db. Two echo characteristics of interest 
are their sharpness and range stability. Pulse 
broadening much greater than that caused by the 
receiver circuitry has not been observed. This 
is in contrast to the large amount of pulse elonga- 
tion found in normal ground backscatter echoes 
propagated via the ionosphere, where energy is 
returned from a large area of the earth, and 
suggests a selective mechanism of wave guidance 
such that only those waves whose normals lie 
very close to the direction of the magnetic field 
are propagated via this mode. The range of an 
echo usually remained constant throughout the 
period during which it was detectable but, as 
may be seen in Table I, echo ranges differed 
from one time of occurrence to another. Figure 
1 is a good example showing the stability of echo 
range as well as the sharpness of echo return. 
This record was made between the hours of 0200 
and 0300 U.T. on June 9, 1960, during which 
time an echo at a range of 27 600 km was found 
to exist for more than 45 minutes, gradually 
weakening in strength and finally disappearing. 
The total propagation path distance of this echo 
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Table I. Summary of one set of data. 














Observation period Echo return Approximate 
- U.T. U. F range (km) 
Date Start Stop Start Stop 
xi April 22 0100 0513 
22 0558 1309 1000 1040 25 000 
' 25 0440 1240 
27 0500 1244 
i= 28 0630 1149 
30 0930 1241 
May 2 0536 1205 
3-4 2235 1226 0140 0142 20 000 
0220 0222 20 000 
5 0456 1225 
1 6 0520 1156 0645 0648 21500 
0700 0704 22 000 
0951 0953 22100 
9 0444 1234 
10 0445 1104 
, 12 0529 1202 0530 0532 21500 
4 0630 0635 21500 
od, 13 0330 0950 0405 0420 24 100 
0800 0802 25 200 
' 16 0344 0907 
ch 18 0325 0846 
23 0346 1005 
25 0308 1026 
26 0254 0302 
28 0312 0906 0802 0815 21600 
31 0315 1006 
June 1 0325 1003 
7 0434 0903 
8 0244 0904 0605 0608 20600 
0637 0640 20600 
2 0650 0653 20 600 
0710 0729 20600 
a- 9 0200 1007 0204 0250 27600 
10 0116 1008 
11 0319 0945 0730 0732 24 500 
15 0307 0805 
i 16 0138 1003 
e 8 19 0409 1105 
& 20 0219 0900 
21 0249 0902 0300 0400 24 000 
s 30 000: 
% 24000) 
; UT 0200 
*) 
0 
FIG. 1. Exospheric backscat- 


ter, June 9, 1960, 0200-0300 


U.T. 





(55200 km) exceeded by at least 13 000 km the 
range of round-the-world echoes which have been 
found by various observations made in the high- 
frequency and the very-low-frequency bands to 
range from 41160 to 42000 km.’*>'% 

The characteristics of the echoes described 
above differ from those reported by Obayashi 
and shown in Fig. 8 of reference 6. His echoes, 
at ranges of about 7000 km, show continuous 
fluctuations in range with time and are not sharp 
but instead are spread out. Their characteristics 
and range are very similar to those of normal 
ionospherically propagated backscatter echoes 
shown by several investigators'*’** and the trans- 
equatorial echoes reported by Stanford observ- 
ers.'® 

All exospheric echoes received at Washington, 
D. C., were found at ranges which were less 
than the calculated field-line length and different 
ranges from echo to echo were observed. Short- 
ening of field lines may be explained in terms of 
the combined effect of the diamagnetism and 
drift currents, observed at altitudes beyond 6 
earth radii (Explorer VI and Pioneer V).' Time 
variations of the intensity and position of the cur- 
rents are, perhaps, sufficient to account quanti- 
tatively for the change in length of the line of 
force even though they probably do not cause a 
detectable change in the strength of the magnetic 
field observed at the surface of the earth. Thus, 
it is possible that this new method of exploring 
the exosphere may be able to provide information 
on the characteristics of the diamagnetism and 
drift currents from observations made on the 
ground. An additional quantity which is measur- 
able by this type of experiment, if our concept 
of the method of wave guidance is correct, is the 
gradient of electron density transverse to the 
magnetic field and in the plane containing the 
magnetic line. Under the conditions of our ex- 
periment (i.e., frequency and latitude) it has 
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been calculated that an electron density gradient 
equal to or greater than 8.1x10~? electron/cc cm 
transverse to the field was needed for wave guid- 
ance. 

Continuation of this experiment, at the same 
and at higher geomagnetic latitudes and trans- 
mitting nearly simultaneously on several differ- 
ent frequencies, is planned to permit further in- 
vestigation of the variations of transverse elec- 
tron density gradient and the effects of the ring 
current. Observations of one-way propagation, 
with transmitter and receiver located at mag- 
netic conjugate points such as New Zealand and 
southern Alaska, are also being considered. 

A more complete discussion of the conditions 
necessary for the high-frequency exospheric 
mode of propagation and the resulting informa- 
tion which can be deduced concerning the exos- 
phere is planned for publication in the future. 
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ELECTRICAL CONDUCTION AND BREAKDOWN IN HIGH-PRESSURE RARE GASES 


R. Forman 


Parma Research Laboratory, Union Carbide Corporation, Parma, Ohio 
(Received March 29, 1961) 


Some novel results have been obtained in the 
study of electrical conduction and breakdown 
effects in high-pressure inert gases (1-300 mm 
of Hg) using a high-temperature thermionic emit- 
ter as a cathode. 

Two types of diodes have been used in these 
measurements. The first is a cylindrical-type 
diode which consists of a 0.010-inch diameter 
filament about 3.5 cm long and a cylindrical 
anode 1 cm in diameter and 1.25 cm long. The 
second type consists of a horizontal filament 
about 1.25 cm long and 0.010 inch in diameter 
with a disk anode, 1.5 cm in diameter, placed 
0.5 cm above the filament. The filament was 
in general of tungsten, although similar results 
have been obtained on tantalum and rhenium. 
The anode material was either tungsten or tan- 
talum sheet. 

The diode assembly was mounted in a Pyrex 
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envelope and evacuated on an ultrahigh-vacuum 
system. The high-vacuum pump consisted of a 
forepump, three-stage water-cooled oil diffu- 
sion pump and a liquid nitrogen cold trap. In- 
cluded in the high-vacuum system was an Alpert- 
type metal valve, a molybdenum getter tube, a 
purified break-seal argon bottle, the tube under 
test, and a Bayard-Alpert ionization gauge. The 
system was baked out at 400°C for 8-10 hours 
and after degassing the pressure in the system 
was below 10~® mm of mercury. The filaments 
were degassed by heating to 2500°-3000°K and 
the anodes were degassed at about 1000°C by 
electron bombardment. 

After the tubes had been processed as de- 
scribed, anode voltage vs anode current vacu- 
um data were taken on these diodes. The fila- 
ments were heated by a dc voltage source. The 
points shown as dark squares in Fig. 1 illustrate 
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FIG. 1. Anode current-voltage characteristics of 
Tube 7B at vacuum and at an argon pressure of 107 
mm of mercury. This tube is a cylindrical diode with 
a 0.01-inch diameter tungsten filament and a 1.0-cm 
diameter tungsten cylindrical anode. The vacuum 
data are shown by the solid squares. The other data 
illustrate the diode characteristics in argon as a func- 
tion of filament temperature. 


these data taken on a cylindrical-type diode with 
a filament temperature of 2800°K. This curve is 
representative of space-charge data and above 
+1 volt is the same for all filament temperatures 
between 2400°-3000°K. However, if a high pres- 
sure of high-purity argon is the ambient, some 
very different anode current vs anode voltage 
curves are obtained. These are illustrated in 
Fig. 1 for the case where the tube was filled 
with argon to 107 mm of mercury. The space- 
charge-limited vacuum data, described previ- 
ously, are in the same figure for comparison. 
The high pressure—high temperature data of 
Fig. 1 have some unusual features. One obvious 
effect is that the current vs voltage curves at 
the high pressure do not follow a space-charge 
relation. They are extremely temperature de- 
pendent. In addition, at temperatures of 2400°K 
and above, the gas breaks down before the anode 
potential reaches the ionization potential of the 
gas. These measurements were made with dc 
on the filament where the negative end of the fila- 
ment was always the common terminal for the dc 
power supply in the anode circuit. Under these 
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FIG. 2. Anode current-voltage characteristics of 
Tube 8A at a given filament temperature with the pres- 
sure of argon as the other parameter. This tube has 
a disk tantalum anode and a 0. 01-inch diameter tung- 
sten filament. 


conditions the maximum accelerating voltage in 
the tube was the anode voltage and this type of 
arrangement eliminated any possibility of having 
the filament voltage contribute to early break- 
down. 

Figure 2 is a plot of the anode current vs anode 
voltage for a diode, with a filament temperature 
of 2800°K. The ambient pressure of argon is the 
variable parameter. This tube had a 0.01l-inch 
diameter tungsten filament and a tantalum disk 
anode. The vacuum line follows the usual V™ 
space-charge equation. One should note that, at 
pressures above 1 mm, there is no doubt that 
the gas breaks down before ionization potential 
is reached. At lower pressures (250 y» and be- 
low), breakdown occurs at ionization potential or 
above. 

Similar results have been obtained in xenon 
and krypton but not in helium or neon. 
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The highest temperature current-voltage curve 
of Fig. 1 resembles those obtained on a gas plas- 
ma with a Langmuir probe and can be analyzed 
on this basis to yield very reasonable values of 
the parameters. This plasma-like behavior 
strongly suggests that there are positive ions 
being generated in the gas below breakdown 
which neutralize the electron space charge. 

From the data in Fig. 1 it is also obvious that 

if this assumption is correct, there are a 
greater number of ions being generated at a 
higher temperature than at the lower ones, be- 


‘these ions is the thermal ionization of the inert 





cause the measured electron current rises with 
increasing temperature. 
A possible mechanism for the production of 






gas at or in the vicinity of the hot filament. In 
addition, it is probable that some of the high 
electrical conduction properties of the gas, 
shown in Figs. 1 and 2, can be partly attributed 
to the generally low collision cross section for 
electrons in the inert gases. 

A more complete report will be published in 
the Journal of Applied Physics. 





NEW ALLOTROPIC FORM OF He*t 


A. F. Schuch and R. L. Mills 


Los Alamos Scientific Laboratory, University of California, Los Alamos, New Mexico 
(Received May 8, 1961) 


The occurrence in solid He* of a change in 
phase from hexagonal closest packing to the 
face-centered cubic structure’? at a temperature 
and pressure above 14.9°K and 1100 atmospheres 
has prompted an investigation of the high-pressure 
and temperature structure of He’. 

This investigation employed the cryostat and 
x-ray diffraction arrangement used in the He* 
study.’ However, because higher pressures were 
required, the He® was solidified in a beryllium 
cell with a slightly thicker wall. The cell wasa 
cylinder with a 0.8-mm bore and a 0.5-mm wall. 
The pressure was measured with a Bourdon gauge 
which was calibrated with a free piston gauge. 
The temperature was measured with a hydrogen 
vapor pressure thermometer fastened to the 
sample cell. A four-stage mercury piston pump 
designed by Edeskuty*® was used to compress the 
gas from 4 psi to the required pressure. The 
diffraction patterns were obtained with copper 
and with molybdenum radiation and were photo- 
graphed on a 45-inch flat film placed 5 cm from 
the sample. During the exposure the cell was 
oscillated through an angle of thirty degrees. 

The He® was analyzed with a mass spectrometer. 
The impurities found were 0.04% He* and 0.04% 
of a mixture of hydrogen isotopes. The remaining 
99.92 % was He’. 

The sample was first prepared by cooling the 
cell to the neighborhood of 20°K and then com- 
pressing the helium into it at the desired pres- 
sure of the experiment. The gas was then solidi- 
fied by slowly cooling it. The intention was always 
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to have the solid at a temperature and pressure 
near the melting curve. 

It was found that He*® under high pressure also 
transforms to the face-centered cubic structure. 
The data from a diffraction photograph of He® at 
18.74°K and 1693 atmospheres are listed in 
Table I. In the table, d is the interplanar spacing 
found and a, is the length of the edge of the cubic 
unit cell derived from it. 

The average cell edge obtained from seven 
photographs of helium at an average temperature 
and pressure of 18.76°K and 1690 atmospheres is 
a, =4.242+0.016 A. This dimension gives a molar 
volume of 11.50 cc, which is to be compared to 
11.54 cc computed from the measurements of 
Grilly and Mills.* 

The hcp-fcc transition occurs at a higher pres- 
sure than it does for He*. Where the transition 
line joins the meiting curve has not been deter- 








Table I. Data from one diffraction photograph. 
hkl d (A) @, (A) 
111 2.425 4.200 
200 2.138 4.276 
220 1.495 4.228 
131 1. 283 4.255 
222 1.221 4.230 


Average = 4.238 
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mined in this study. It was observed, however, 
that He® at 15.98°K and 1341 atm still has the 
hexagonal close-packed structure. The cell di- 
mensions and molar volume at this condition are: 
C, = 4.986 A, a,=0.046 A, c,/a,=1.637, and 
V=12.07 cc. The computed* molar volume is 
12.09 cc. 

The new phase has been designated the y phase. 

We wish to acknowledge the many contributions 
of Donald F. Clinton and the members of his staff 
shop in the construction and repair of items of 


equipment used in this research. 





tWork performed under the auspices of the U. S. 
Atomic Energy Commission. 
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SIMPLE MODEL FOR THE MAGNETIC BEHAVIOR OF SUPERCONDUCTORS 
OF NEGATIVE SURFACE ENERGY 


B. B. Goodman 


Institut Fourier, Grenoble, France 
(Received April 12, 1961) 


It is well known that for certain superconductors 
the magnetic field H,* required for the complete 
suppression of the superconducting state consid- 
erably exceeds the critical field H,. Interest in 
this question has been revived by the discovery 
that, well below their transition temperatures, 
exceptionally high values of H,.* are found for 
Mo,Re* (= 15 kgauss) and Nb,Sn? (> 88 kgauss). 
While inhomogeneous strains and chemical in- 
homogeneities in a superconductor are known 
to influence its behavior, probably owing to a 
resulting spatial variation of the surface energy 
of a normal-superconducting boundary, it is the 
purpose of this Letter to propose a simplified 
model for the magnetic behavior of a homogeneous 
strain-free superconductor, in which a spatially 
independent negative surface energy could lead 
to large values of H,*/H-- 

Suppose that a superconductor of negligible 
demagnetizing coefficient divides itself, in the 
presence of a magnetic field H, into thin normal 
and superconducting laminas, of thickness d, 
andd,, respectively. Let us add to the resulting 
expression for g, the Gibbs free energy per 
unit volume,* an additional contribution arising 
from a configurational interphase surface energy, 
4’H,”/87. One then obtains 


cs q tanhp\] .. 
a= °/8 7) }-[1 “; -7 (1 a+ 


where H.”/81 is the difference in free energy 


» (1) 





per unit volume between the normal state and 
the superconducting state in zero field, h=H/H,, 
x =ds/(dy+ds) is the fraction of the specimen in 
the superconducting state, p =d,/2 A, and g=A’/,, 
A being the penetration depth of the magnetic field. 
Equation (1) supposes that the London equation 
is obeyed in the superconducting laminas, a hy- 
pothesis we shall comment on later. As London 
pointed out,® the condition for Eq. (1) to lead to 
the well-known magnetic behavior of a so-called 
ideal superconductor is that g>1. Hence A=A’-i, 
the total surface energy usually measured experi- 
mentally, is then positive. Pippard*’® has sug- 
gested that A’ should be of the order of —, the 
coherence distance in the superconducting state. 
For pure superconductors, é tends towards a 
limiting value £, characteristic of the metal; 
while for very impure superconductors (/<£,), 
— is of the order of J, the electronic mean free 
path.® 

For g<1, corresponding to a superconductor 
of negative total surface energy, Eq. (1) leads 
to the existence of two transition fields. Forh 
less than the lower transition field h’, given by 
h’=q, one finds x =1 and p=, so that here the 
specimen exhibits a Meissner effect, while for 
h>h’'’, where h”’ is given by 


g=h'’+(1 -h’”) tanh “(h’’), (2) 


the normal state is stable. For q«1, Eq. (2) 
becomes h’’=2/3q. For intermediate fields the 
laminar structure is stable, the thickness of the 
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FIG. 1. Reduced magnetiza- 
tion curves for the present 
model of a superconductor, for 
different values of q. 








superconducting laminas being determined by the 
equation gh~* =tanhp - psech’p, while that of the 
normal laminas remains negligible (x=1). The 
magnetization M, given by -47M =xH[1 - (tanhp)/p], 
leads to the curves given in Fig. 1 for different 
values of g. We may note that, since the surface 
energy does not contribute to g either for h =0 

or for h>h’’, the area under each magnetization 
curve is independent of q and satisfies the ther- 
modynamic condition, 


oo 
f MdH =-H ?/81. 


It is tempting to correlate the prediction of 
two transition fields for q<1 with the observation 
of Doidge’ and Davies® that in dilute tin alloys 
the magnetic and resistive transitions become 
distinct when /<A. However, the magnetic tran- 
sitions observed by these authors went practi- 
cally to completion at the lower transition field. 
The inadequacy of the present model, when 7 is 
not small compared with A, is hardly surprising, 
since here the London equation is known to be 
inexact,° and furthermore it is physically un- 
reasonable to assume, as we have done, that 
the normal laminas can be of negligible thickness. 
However, these limitations should be much 
less serious for very impure superconductors, 
for which / is small, and Fig. 2 shows magne- 
tization curves based on the present model which 
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have been fitted to the hysteresis-free results 
obtained by Calverley and Rose-Innes® for a 
single mixed crystal of composition Ta,,Nb,,. 

It has here been assumed that the fields H,.* 
quoted for the complete suppression of super- 
conductivity correspond to the predicted discon- 
tinuous disappearance of M ath’’. In fact, it is 
more probable that, owing to a field dependence 
of A which we have not considered, the observed 
magnetization went to zero continuously, as has 
been found for superconducting thin films and 
colloids.**"* The results of these authors sug- 
gest that the effect of allowing A to vary with H 
would probably be to distort the magnetization 
curves in Fig. 1 in such a way as (a) to yielda 
continuous disappearance of M, h’’ being in- 
creased by at most a factor of about V2, and 

(b) to decrease the fitted values of H, in Fig. 2 
by not more than 10%. The correction required 
by the field dependence of A has therefore been 
ignored. 

By extrapolation of the values of H, to 0°K one 
obtains H,=1600+100 gauss, so that H,/T,.=230 
gauss deg™, which is within 20% of the values 
found for this quantity for pure Ta and Nb. The 
value obtained for g, which is expected to be close 
to 1/A, may be discussed by using standard ex- 
pressions relating A, &,, 1/o (o being the conduc- 
tivity of the specimen =2.5x10™ emu), and T, 
(its transition temperature = 6.9°K) to the proper- 
ties of the Fermi surface.” These lead to 4=800A 
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and /=76(S/S,)"* A, where S is the free area of 
the Fermi surface in the alloy and S, the area it 
would have if the five conduction electrons per 
atom were perfectly free. Thus the result qg=0.16 
suggests that (S/S,) should be somewhat less than 
unity, which does not seem unreasonable in view 
of the values of this order of magnitude which 
have been found for several other metals.** 

The present model therefore seems to give a 
plausible account of the results of Calverley and 
Rose-Innes, and it confirms the impression that 
their magnetization curves represent a new type 
of magnetic behavior characteristic of homoge- 
neous strain-free superconductors of negative 
surface energy. The magnetic field needed for 
the complete suppression of superconductivity, 
h"'H,, appears as a true thermodynamic variable, 
rather than as a property associated with certain 
individual lattice flaws. According to this model, 
for large h’’, h’’=2/3q=22/3&, while using a 
different model, which takes into account the 
field dependence of A, Pippard® finds a value 
which is only slightly larger: h’’= (fv 6)Ay = 0/mé; 
here f depends on the details of the model and lies 


between 1 and V2. Since for 1A, we have ,«/~“” 


and §=/, we may expect that, other things being 
equal, h’’x<]~**, The persistence of supercon- 
ductivity in high magnetic fields, brought about 


H (gauss) 


by the existence of a negative surface energy, 
should therefore be favored by a short electronic 
mean free path, and thus a high resistivity, in 
the normal state. 
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EFFECT OF MAGNETIC FIELD ON THERMAL CONDUCTIVITY 
AND ENERGY GAP OF SUPERCONDUCTING FILMS* 


D. E. Morris and M. Tinkhamt 
Department of Physics, University of California, Berkeley, California 
(Received April 27, 1961) 


We have measured the change in the thermal 
conductivity of superconducting tin and indium 
films upon application of a magnetic field in the 
plane of the film. These experiments were under- 
taken to explore the dependence of the energy 
gap and the penetration depth upon magnetic 
field, and to determine the thermodynamic order 
of the field-induced superconducting transition 
in films. 

Sample films are evaporated onto glass sub- 
strates held at 77°K, after ion-bombardment 
cleaning of the substrate surface. During an ex- 
periment, a typical temperature drop along the 
film is 0.1°K with ~10~" watt total heat transport 
along the sample film and substrate. Because 
of its relatively large thickness (~0.01 cm), the 
substrate carries most of the heat despite its 
very low thermal conductivity. The change from 
superconducting to normal alters the temperature 
drop by ~0.001°K. Changes of ~2 microdegrees 
in the temperature drop along the film are detect- 
able. This extreme differential temperature 
sensitivity is achieved by use of a carbon re- 
sistor heater-thermometer at one end of the 
film, connected in a bridge circuit with an iden- 
tical reference thermometer located at the op- 
posite end of the film. A copper rod connects 
the latter end to a helium bath whose tempera- 
ture is controlled to ~0.001°K by a vapor pres- 
sure regulator. Balance of the thermometer 
bridge is measured by a dc microvoltmeter. 
Fields up to 1000 oersteds are produced by a 
water-cooled Helmholtz pair, and the angle of 
the field with respect to the surface of the film 
can be adjusted to within 0.02° by auxiliary coils. 
In practice, field tilts of +0.4° produce negligi- 
ble changes in the experimental results. To 
reach the high critical fields of very thin films, 
we use an iron-core electromagnet. The film 
thickness is estimated from the measured criti- 
cal field and published data’ relating the critical 
field to thickness. 

Experimental results for films of tin and indi- 
um of thickness ~ 700 A are shown in Fig. 1. 

For both metals the thermal conductivity in- 
creases nearly as H? up to the critical field. At 
the critical field, the thermal conductivity in the 
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FIG. 1. Change of thermal conductivity of thin super- 
conducting films with magnetic field. 


superconducting state joins smoothly onto the 

field-independent conductivity of the normal 

metal. This indicates that the superconducting 

transition in a magnetic field is second order 

for films of this thickness, in contrast to the 

first-order transition of bulk superconductors ; 

in a field. _ 
Orientation of the field parallel and perpen- : 

dicular to the direction of heat flow in the film ; 

produces the same effect upon the thermal con- 

ductivity to within the experimental accuracy, 

for films ranging between 700 and 2800 A. Since 

the induced diamagnetic currents in the super- 

conductor are perpendicular to the field, our 

data also indicate that the thermal conductivity 

does not depend upon the angle between the ther- 

mal gradient and the diamagnetic current. Thus, 

if the observed variation of the thermal conduc- 

tivity is explained in terms of a decrease of the 

energy gap of the BCS theory’ with increasing 

field, the modified gap remains essentially iso- 

tropic. This result is interesting in view of the 

theoretical prediction’ of a pr* Vgrig¢ term in the 

excitation spectrum of a superconductor carry- 

ing a uniform current. The absence of an ob- 

servable effect of this sort might be due to the 
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fact that in these films the two equal and oppo- 
site surface currents are only slightly separated 
in space. Alternately, the short-mean-free path 
in these “dirty” samples may upset the simple 
prediction. The lack of anisotropy is important 
in evaluating the results of Giaever and Megerle* 
concerning the magnetic field dependence of the 
tunnel effect between superconductors, since 
their method can only determine the energy gap 
normal to the surface of the film, which is per- 
pendicular to the directions of both the applied 
magnetic field and induced current. 

The results of Fig. 1 may be used to compute 
a dependence of the energy gap upon field in a 
simple way if the electronic term Kg dominates 
the thermal conductivity, and if Ke is primarily 
limited by electron scattering by lattice imper- 
fections. Under these conditions the ratio of the 
thermal conductivity in the superconducting state 
Kes to that in the normal state Kp, is given by 
Bardeen, Rickayzen, and Tewordt® as their Eq. 
(3.6). This may be written 


K °o 00 
f° = Gle,/T)= [ ery /es ae /f E(0f/8E)dE, 
en Eo 0 


where f(E/kT) is the Fermi function. If we as- 
sume that the effect of a field upon the supercon- 
ducting state can be adequately represented as 

a change of the €, of BCS, then the observed 
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FIG. 2. Magnetic field dependence of superconduc- 
ting energy gap computed from data of Fig. 1. 


quantity Ax plotted in Fig. 1 should be given by 
K, ,(#) -K,,© : Gle )(H)/kT} - G[e 9(0)/eT} 
1- Gle )(0)/k T| 
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This relation was used to invert the experimental 
data on indium III to yield €,(H), with the result 
shown in Fig. 2. Because G(e,/kT) initially drops 
only quadratically as €,/kT increases from zero, 
the gap.must drop more steeply to zero near H c 
than the approach of Kg,(H) to Ken there. The 
difficulty of distinguishing between a broadened 
transition because of film inhomogeneity and a 
truly rounded approach to Key, because of the 
form of G(e,/kT) results in an uncertainty in the 
choice of H-. The bars on the plotted experi- 
mental points in Fig. 2 indicate the size of this 
uncertainty. 

Our results for the dependence of the energy 
gap upon magnetic field in a thin film are quali- 
tatively similar to those from the more direct 
measurements of Giaever and Megerle,* although 
their data were taken on a thicker film at a much 
higher reduced temperature (7/T,) and in a dif- 
ferent metal (aluminum). From our data in Fig. 
2, €9(H)/€,(0) seems to approach 1 -(H/H,)? at 
low temperatures, but it is moderately well fit- 
ted at T/T, = 0.63 by €,(H)/e,(0) =[1 - (H/H,)?]”. 
The latter form is that given by Douglass® based 
on the Landau-Ginzburg-Gorkov theory, which 
is expected to hold near T,. 

With thicker films, new features appear as 
shown in Fig. 3. The thickest film, tin II with 
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FIG. 3. Change of thermal conductivity of supercon- 
ducting films of intermediate thickness with magnetic 
field. 
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d=2800 A, has Ax ~H? up io the vicinity of H,, 
where it reaches only ~25%; then K appears to 
rise discontinuously to K,, indicating that the 
transition is first order in films of this thickness. 
By our method of analysis, we find that €,(H)/ 
€,(0) drops abruptly from ~ 0.83 to 0. If we apply 
the theory of reference 6, the observed critical 
field ratio for this film, H,./H,(bulk) =1.53, im- 
plies that d/A =3.6 and €,(H_,)/€,(0)=0.78. The 
agreement with the above value of 0.83 is quite 
good considering the limited accuracy and the 

fact that the gap will not be completely uniform 

in a film as thick as this one. The film indium II, 
of intermediate thickness d= 1800 A, is still thin 
enough to display a second-order transition, but 
the increase of K near H, is much faster than -. 
This behavior may be explained in terms of an in- 
crease of the penetration depth A as the gap de- 
creases, as would be expected on the basis of sum- 
rule arguments.” Since the film thickness is small- 
er than a coherence length &,, we expect the energy 
gap to be nearly independent of position in the film. 
We may then apply the approach of Pippard,® con- 
sidering the balance between condensation energy 
and magnetic energy to determine the gap (or 
order parameter) as a function of field. Taking 
account of the dependence of film susceptibility 

on (d/,) and of the dependence of \ on €,, one is 
led to expect a change in the form of €,(H/H,) of 
the observed sense for films withd>a. A sim- 
ilar conclusion follows from the explicit theory 

of reference 6. Only in a film thin enough (d<,) 
so that the penetration is always nearly complete 
will there be a definite limiting relation between 
€,(H) /€,(0) and H/H,. 
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In indium II at the lowest temperature we see 
an apparent initial decrease of thermal conduc- 
tivity with field. We believe that this represents 
a decrease in the phonon conductivity K g? pro- 
duced by increased electronic scattering of pho- 
nons as the energy gap is decreased. This inter- 
pretation is reasonable, as the initial decrease 
is seen at the lowest temperatures and in a rela- 
tively thick film, both conditions favoring Ky 
over Kg, and it is observed in indium, having a 
lower Debye © than tin. The observed K, indi- 
cates an electronic mean free path limited by 
defect scattering to/~200 A. Then Kg, limited 
only by boundary scattering, should be compar- 
able in magnitude to Kg, as required for the 
validity of this interpretation. 
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GROUND-STATE ENERGY OF AN ELECTRON GAS IN A LATTICE OF POSITIVE POINT CHARGES 


A. Bellemans and M. De Leener* 


Faculty of Sciences, University of Brussels, Brussels, Belgium 
(Received February 27, 1961; revised manuscript received May 4, 1961) 


The ground-state energy of an electron gas ina 
rigid lattice of positive point charges has been 
studied using the quantum statistical formalism 
of Bloch and de Dominicis.’»? The grand partition 
function of the system was expanded in powers of 
the coupling parameter between (i) the electrons 
and the lattice and (ii) pairs of electrons; this 
parameter is of course the square of the elec- 
tronic charge. On account of the long-range 
forces, a certain number of divergences are 
met within this expansion; as usual one gets rid 
of them by (a) taking into account the electroneu- 
trality of the whole system, and (b) summing in- 
finite series of divergent terms. 

Finally, taking the limit of zero temperature, 
one obtains the ground-state energy which can 
be conveniently expressed in terms of the para- 
meter r,, related to the density p of the electron 
gas by 


a (3/42)“p-“’. (1) 
Using as units the Bohr radius for r, and the ryd- 
berg for the energy, we found the following ex- 


pression of the ground-state energy for the case 
of a simple cubic lattice: 


e( oe i -3CA"+ 
"s 5\4/) v2 3\a) \4 i 


ys 
-- 6) ~ +{0.0622 Inr_- 0.096} 
Ss 








1\4 
s 
1 } 1 (1 -x? l+u 
ea} ia "li-e +f 
+ O(r Inr |) ; (2) 


v is the ratio of the charge of a lattice point to the 
electron charge; the vector ti is equal to (1/3 v)“*i, 
where (n,,my, 2) are integral numbers; x is a 
constant defined as 


xo-ts L fh as, (3) 


where X is a lattice vector. 
The physical meaning of the various terms of 
(2) is the following: (a) The first term is the 


kinetic energy of a free electron gas; (b) the sec- 
ond term is the Madelung energy [see Eq. (3)] of 
a lattice of positive point charges in a uniform 
negative background; (c) the third term is the ex- 
change energy of the electron gas, as calculated 
by Wigner®; (d) the fourth term is the correlation 
energy of the electron gas, as calculated by Gell- 
Mann and Brueckner*; and (e) the fifth term is the 
polarization energy of the free electron gas by the 
lattice, in the linear approximation. 

The sum of terms (a), (c), and (d) in Eq. (2) is 
just the ground-state energy of an electron gas in 
a positive continuum; terms (b) and (e) arise from 
the lattice. 

Expression (2) is plotted on Fig. 1 (curve I) for 
the special case y=1; strictly speaking, this 
would correspond to “metallic hydrogen.” The 
importance of the “structural” terms (b) and (e) 
in €(75) is apparent when one compares curve I 
to curve II, the latter corresponding to an electron 
gas in a uniform positive background.* (For body- 
centered and face-centered cubic lattices, curves 
very similar to curve I are obtained.) 

The ground-state energy of “metallic hydrogen” 
has been investigated before by Wigner and Hunt- 
ington®; they first looked for approximate Bloch 
wave functions for independent electrons and next 
calculated exchange and correlation corrections. 
This method of approach is rather different from 
ours, which treats electron-lattice and electron- 
electron interactions in a symmetric way. It is 
then quite remarkable that the curve for €(7.) 
derived by Wigner and Huntington is very close 
to our curve I, as can be seen from Fig. 1. The 
abcissas and ordinates at the minimum are, re- 
spectively: r,=1.64, €(r,) =-1.05 and 7, =1.59, 
€(r,)=-0.99. We may further expect that the 
exact minimum value of €(7,) for “metallic hydro- 
gen” is near to (and probably somewhat below) -1 
rydberg. This indicates that expansions for €(r,) 
like (2) converge satisfactorily in the region 7, <2. 
This seems to be supported by the calculations of 
Dubois® for the term O(r, Inr,) of the correlation 
energy of an electron gas in a positive uniform 
background and also by our own rough estimations 
of the polarization energy O(r,) of the electrons 
by the lattice. 

A detailed account of the results reported above 
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will be published elsewhere. The authors wish to 
thank Professor I. Prigogine for his constant in- 
terest in this work. 





*Aspirant au F.N.R.S. 
'C, Bloch and C. de Dominicis, Nuclear Phys. 7, 


459 (1958). 
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SPLITTING OF THE EMISSION LINES OF RUBY BY AN EXTERNAL ELECTRIC FIELD 


W. Kaiser, S. Sugano, and D. L. Wood 
Bell Telephone Laboratories, Murray Hill, New Jersey 
(Received May 9, 1961) 


The fluorescence emission of ruby (ALO, +Cr**) 
has been investigated in great detail in recent 
years and it is now well understood.’ In highly 
diluted pink ruby with chromium concentrations 
below 0.05%, the emission consists essentially 
of two lines at R, = 14421 cm™ and R,= 14450 cm™ 
corresponding to transitions from the split com- 
ponents of the ¢,° E state to the ground state *A,, 
respectively. Since the ground state is split by 
0.38 cm™,? both emission lines are actually dou- 
blets, which can easily be observed at 77°K with 
a high-resolution grating instrument. 

We have superimposed an external electric 
field upon the internal crystal field and observed 
a new additional splitting of the emission lines, 
which is much larger than expected from a super- 
ficial consideration of the Stark effect. Ruby 
plates of 1-mm thickness with a chromium con- 
centration of approximately 0.01% were cut from 
a large single crystal parallel and perpendicular 
to the optic axis. The samples were immersed 
in liquid nitrogen and illuminated with a 100- 
watt Hg lamp. The emission of the crystal was 
focused on the entrance slit of a Bausch and 
Lomb grating spectrometer with a 50000-line/in. 
grating. With this system the ground-state split- 
ting of the R, and R, lines was clearly visible 
(see Fig. 1 top). Upon application of an external 
electric field parallel to the optic axis of the 
ruby, a surprisingly strong effect on the emis- 
sion pattern was observed. Both components of 
the R, and R, lines exhibit a distinct splitting 
which is shown in Fig. 1 for the R, line for sev- 
eral electric fields. On the top (no field) and the 
bottom (£,=1.58%10° v/cm) of this figure, actual 
microphotometer traces of the photographic plate 
are presented; at intermediate fields, only the 
positions of the individual components are indi- 
cated. The frequency shift is so large that the 
two center lines pass the midpoint between the 
two original R, components at a field of about 
65 kv/cm and cross over at still higher fields. 

It is interesting to note that the center of the 
split lines does not shift with applied electric 
field; furthermore, both components of the R, 
(see Fig. 1) and R, lines split by the same a- 
mount. It can be seen from Fig. 2 that the split- 
ting Av is directly proportional to the electric 
field. At a field of 1.7x105 v/cm, Av reaches 


a value of one wave number (3 x10'° cps). With 
increasing electric field, the emission lines 
slightly broaden. No detailed study was made 

of this effect. Polarizers in the light beam indi- 
cated that the polarization is unaltered by the 
electric field. In specimens with the electric 
fields perpendicular to the optic axis, no change 
in the emission pattern could be observed for 
fields as high as 1.7 10° v/cm. 

The observed splittings of the R, and R, fluo- 
rescence lines in the presence of an electric 
field are not normal Stark splittings. This is 
because both the initial and final states of the 
fluorescence are Kramers doublets whose de- 
generacies may not be removed in the presence 
of any electric field. In what follows, we shall 
show that the apparent splittings are due to the 
opposite Stark shifts of the energy levels of chro- 
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FIG. 1. Splitting of the R, line in an electric field 
parallel to the optical axis. 
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FIG. 2. Pseudo-Stark splitting 4v of ruby vs ap- 
plied electric field Ey. Data are obtained from the 
Ri (e) and R» (x) lines. 


mium ions at different lattice sites. 

It is well known that in ruby there are two 
kinds of chromium sites, A and B.* Transforma- 
tion of the A site into the B site is only possible 
by symmetry operations involving inversion 
around the chromium ion, while transformation 
of A to A or B to B is achieved solely by trans- 
lation or by both translation and rotation around 
the optic axis. The A and B sites are energeti- 
cally equivalent in the absence of an electric 
field. This circumstance may be simply ex- 
pressed by the fact that the odd-parity crystal 
field associated with the A site, V,,(A), is equal 
to minus the odd crystal field at the B site, 
V,,(B) = -V,(A). 

Let us consider the local electric field acting 
on the chromium ions. The local field may be 
written in the form, a+Eo, where E, is the ap- 
plied electric field and a is generally a tensor 
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determined by the dielectric properties of crys- 
tals. Polarization effects due to ionic displace- 
ments and deformation of electron clouds in crys- 
tals are taken into account in the tensor a. For 
a uniaxial crystal such as ruby, the term in the 
Hamiltonian for the interaction of the chromo- 
phoric electrons with the local electric field pro- 
duced by the external applied field is given by 


Vv. =U, [a 71X02 +a (x,E 5, +9, Eq.) , 


where x, y, andz are the electron coordinates 
and the z axis is taken along the optic axis. 

Since V, belongs to odd parity, the crystal 
field levels of chromium ions are perturbed by 
the combined action of V, and V;. This perturba- 
tion produces the Stark shift. Because of the op- 
posite sign in the odd-parity crystal fields at the 
A and B sites, the Stark shifts for these sites 
are expected in opposite directions. This gives 
the pseudo-Stark splittings observed in our ex- 
periment. Furthermore, since V,, corresponds 
to the z component of the equivalent irreducible 
tensor,‘ a simple symmetry consideration shows 
that the Stark shifts are vanishing when E, is 
perpendicular to the optic axis. 

The accurate calculation of this type of Stark 
shift is as difficult as that of the absolute intensi- 
ties of the crystal field spectrum. It is, however, 
possible to estimate the order of magnitude of 
the Stark shift Av by using the relation Av~ay 
xEQz (Vy) (z)/4, where (Vy) is the matrix element 
of Vy, between the state relevant to the fluores- 
cence and the high excited state of odd parity, 

{z) the similar matrix element of z, and A the 
energy difference between the two states appear- 
ing in the matrix elements. Putting a “es Be 

(Vy) ~ 1000 cm™, (z)~1 A, and A~50000 cm™, 
we obtain Av~0.1 cm™ at Ep, =10° volts/cm. 

The above-mentioned Stark shifts are expected 
in both the initial and final states of the fluores- 
cence. Therefore, the observed magnitude of 
the pseudo-Stark splitting should be compared 
with twice the difference of the Stark shifts in 
the initial and final states of the fluorescence. 
This would not change the estimated order of the 
magnitude, as the shifts are not expected to be 
the same in both states. 

The other characteristic features expected 
from our theory are as follows: 

(1) The center of the split components does not 
shift (no pseudo-Stark shift) because the crystal 
symmetry gives |V,,(A)| =|V,,(B)|. Here we neg- 
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lect contribution from the V,” term which is very 
small for the fields achieved in our experiments. 

(2) The magnitude of the splittings is expected 
to be the same for the R, and R, excited states 
because the Stark term is independent of spin. 
The same is true for the initially split ground 
levels having M, =+3/2 and M, =+1/2. 

(3) Optical anisotropies associated with the 
split components are identical with each other 
and also with the anisotropy of the original un- 
split line. All these points are in agreement with 
our experiments. 

Although we have not given any detailed dis- 
cussion on the local electric field, experiments 
of this type are very important because they allow 
a direct investigation of the local electric field. 
The quantitative discussion on the local field is 
only possible after the detailed knowledge of the 
odd-parity excited states is obtained. 

The effect discussed here is expected to be 





extremely fast (comparable to the optical fre- 
quencies of the lattice), and high-frequency 
modulation of the ruby emission can be antici- 
pated. The application of an ac electric field 
should allow frequency modulation of the output 
of a ruby optical maser.°® 
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ANISOTROPY OF THE MAGNETOELECTRIC EFFECT IN Cr,0, 


Vv. J. Folen, G. T. Rado, and E. W. Stalder 
U. S. Naval Research Laboratory, Washington, D. C. 
(Received May 8, 1961) 


On the basis of earlier considerations of mag- 
netic crystal symmetry effects by Landau and 
Lifshitz,! Dzyaloshinskii? showed that a magneto- 
electric effect (i.e., a magnetic moment produced 
by an externally applied electric field or an elec- 
tric moment produced by an externally applied 
magnetic field) may exist in antiferromagnetic 
Cr,0,. In his development of the thermodynamic 
potential which is compatible with the magnetic 
symmetry, Dzyaloshinskii obtained the relation’ 

4 %... = -a ye 2" - aE H. +8 BD, (1) 
where ,,¢ is the magnetoelectric part of the 
thermodynamic potential, and a and a, are the 
magnetoelectric parameters for the directions 
parallel and perpendicular to the crystallographic 
c axis, respectively. From this potential, one 
obtains (using 478 /eH = -B) the linear relations 

B, =a a”? B.=aE, and B, “af (2) 
As a result of measurements on a nonoriented 
Cr,0, single crystal, Astrov® obtained evidence 
that the magnetoelectric effect does exist in Cr,O,, 
and that the effect is a linear function of the elec- 




















tric field in accord with Eq. (2). 

In this Letter, we report measurements of the 
magnetoelectric effect (resulting from an applied 
electric field) in two x-ray-oriented single-crystal 
disks of Cr,O,, and show that a large anisotropy 
exists both in the magnitude and temperature de- 
pendence of this effect. One of the disks was fab- 
ricated with its plane surfaces parallel to the 
crystallographic c axis, and the second disk with 
its plane surfaces perpendicular to the c axis. In 
the former case, the surfaces were parallel to the 
(110) plane. The plane surfaces were coated with 
silver paint to provide uniform electrical contact 
with the specimen. The apparatus used to measure 
the magnetic moment was similar to that used by 
Astrov.* However, we found it necessary to make 
the measurements in vacuum in order to avoid 
electrical breakdown in the vicinity of the elec- 
trodes. The magnetic moment was obtained by 
applying to the sample a 1000-cps alternating 
electric field and using a narrow-band amplifier 
of approximately 0.15-yv sensitivity to measure 
the voltage induced in a 10000-turn pickup coil 
surrounding the sample. The largest values of a, 
and a correspond to induced rms voltages of 







607 








VoLuME 6, NUMBER 11 


PHYSICAL REVIEW LETTERS 


June 1, 1961 





8.35 uv and 124.7 pv, respectively. Electrostatic 
shielding was provided by inserting a grounded 
brass cylinder between the pickup coil and the 
assembly consisting of the sample and electrodes. 
A Lavite disk having approximately the same con- 
ductivity as the Cr,O, disks was used to show that 
the electrostatic shielding is adequate and that 
extraneous current effects are absent. Low tem- 
peratures were obtained by immersing a cylinder 
containing the sample and electrodes into a liquid 
nitrogen bath. High temperatures were obtained 
with a furnace or heated water bath. In all the 
measurements the temperature was determined 
by means of a copper -constantan thermocouple 

in thermal contact with one of the silver electrodes 
on the sample. The copper thermocouple wire 
also served as the ground connection for this elec- 
trode. Considerable care had to be taken to avoid 
electrical breakdown in the vicinity of the elec- 
trodes, particularly at the lower temperatures. 
Such breakdown even occurred when a helium at- 
mosphere was used in the apparatus, and conse- 
quently the measurements had to be made in 
vacuum so that no breakdown occurred. 

Figure 1 shows the results of our measurements 
on the temperature dependence of a, and a I" It is 
seen that there is a large anisotropy of the mag- 
netoelectric effect. The magnitude of a, and its 
temperature dependence is similar to that obtained 
by Astrov.* His measurements (which were made 
on a nonoriented crystal) extended down to 250°K 
only. The curve a y» on the other hand, shows 
that the magnetoelectric effect goes through a 
maximum and then becomes zero at 97°K. At 
slightly lower temperatures it appears to rise 
again. In both samples, the magnetoelectric ef- 
fect is zero at 306°K, and remains approximately 
zero at higher temperatures. Since antiferromag- 
netic susceptibility measurements by McGuire 
et al.* show the Néel temperature to be 307°K, it 
appears that the magnetoelectric effect does not 
exist in the paramagnetic region, in accord with 
the symmetry considerations of Landau and Lif- 
shitz. Measurements of the magnetoelectric ef- 
fect in the two disks as a function of applied elec- 
tric field showed that the effect is linear [in accord 
with Eq. (2)] for fields ranging from zero to 3200 
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FIG. 1. Temperature dependence of the magnetoelec- 
tric parameters @ , and a I The @’s are dimensionless 
in the Gaussian units used. 


rms volts/cm. The values of E, and E,, used in 
obtaining the results shown in Fig. 1 were 2160 and 
2050 rms volts/cm, respectively. 

A theoretical discussion of the temperature 
dependence of the @’s is given in reference 5. 
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MECHANISM OF THE MAGNETOELECTRIC EFFECT IN AN ANTIFERROMAGNET 


George T. Rado 
U. S. Naval Research Laboratory, Washington, D. C. 
(Received May 8, 1961) 


A crystal placed in an applied electric field 
may experience a magnetization which is pro- 
portional to the electric field. This magneto- 
electric effect was first observed experimentally 
by Astrov’ in a monocrystal of antiferromagnetic 
Cr,0,. Even earlier, Landau and Lifshitz? had 
pointed out the possibility that this effect (and its 
inverse) may exist in some spin-ordered materi- 
als, and Dzyaloshinskii® had predicted specifi- 
cally that it could exist in Cr,0,. Both of these 
theories are based on symmetry arguments and 
thermodynamics only. 

In the present Letter we propose a simple phe- 
nomenological theory of the temperature depend- 
ence of the magnetoelectric effect. This theory 
correctly predicts the experimental observations 
presented in the accompanying Letter‘ for both 
the parallel (||) and perpendicular (1) orienta- 
tions of the electric field with respect to the 
axis of spontaneous antiferromagnetism (z axis). 
In addition, we tentatively suggest a possible a- 
tomic mechanism which conforms to the restric- 
tions imposed by our phenomenological theory. 
Rough numerical estimates based on this mecha- 
nism explain the observed values of the magneto- 
electric effect to within an order of magnitude. 

We begin by introducing a fictitious magnetic 
field h which we define by the requirement that 
the magnetization M induced by h be identical 
with that induced by the applied electric field 
E*. Since the quantity a, /47, which may be 
called the “parallel magnetoelectric susceptibil- 
ity,” is defined® by a y = B, /E,°, our definition 
of h leads to 


a =40M/E,“ = 40x )h,/E,", (1) 
where the applied magnetic field H? is assumed 
to be zero (so that B, = 41M,) and the parallel 
antiferromagnetic volume susceptibility is de- 
noted by Xq- If we recall that experimentally 
M, is proportional to E,%, and that the magneto- 
electric effect vanishes unless the material un- 
der consideration possesses an ordered spin ar- 
rangement, it seems reasonable to assume that 
hz has the form 

h =a,E Ss), (2) 

z 2° 2’ av 


where (S_),y, which is proportional to the 
zero-field sublattice magnetization M,, rep- 
resents the thermal average of the expectation 
value of the z component of the ionic spin, and 
a, is an essentially temperature-independent 
constant of the material. By combining Eqs. (1) 
and (2) we obtain 
=4 

——r"e 
and on the basis of analogous arguments we ob- 
tain 


(S.) a.» (3) 


a, =4na,x,(S_) ae (4) 
The temperature dependence of a, anda, is 
seen to be determined by x Mo and xX,Mo, re- 
spectively. If we express these quantities in 
terms of the equations resulting from the Néel- 
Van Vleck molecular field theory® of a two-sub- 
lattice antiferromagnet, we can make the fol- 
lowing predictions: As the temperature T in- 
creases from 0°K to the Néel temperature, Ty, 
the quantity a, should decrease monotonically 
from some finite value to zero in the manner of 
the sublattice magnetization; the quantity a}, 
on the other hand, should increase from zero to 
some maximum value and then return to zero. 
Both of these qualitative predictions have been 
verified experimentally.‘ As shown in Fig. 1, 
even the quantitative agreement between the 
calculated and measured q’s is satisfactory. 
The existing discrepancies between theory and 
experiment are no worse in the case of the a 
measurements (Fig. 1) than in the case of the 
x measurements® (Fig. 2 of reference 6). In 
both cases these discrepancies may at least in 
part be attributed to the fact that a two-sublat- 
tice model is not strictly applicable to Cr,O,. 
According to the above considerations, elec- 
tric fields should have “tuning” (and other) ef- 
fects on antiferromagnetic resonance, and these 
effects should be calculable by simply using our 
fictitious h in place of (or in addition to) the 
actual H@ throughout the usual resonance equa- 
tions. We also suggest that by applying a strong 
magnetic field parallel to the antiferromagnetic 
axis and making use of the Néel “spin-flop” ef- 
fect,® it should be possible to “switch” from a 
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FIG. 1. Comparison of calculated and measured 


magnetoelectric susceptibilities. The theoretical 
curves are based on Eqs. (3) and (4) and the following 
numerical values: Ty =360°K, 6(=Curie-Weiss con- 


stant)=-0.38Ty, ay =38.7*10~%, and a, =1.32 107%. 


The experimental points are taken from Folen, Rado, 
and Stalder, reference 4. 


measurement of a i to a measurement of a,. 

In attempting to explain the physical origin of 
the magnetoelectric effect, we think of the elec- 
tron cloud of each magnetic ion as being dis- 
torted by E*. As a result of spin-orbit coupling, 
the spins are “aware” of this distortion and 
hence they may give rise to an induced magneti- 
zation. To construct a specific model, we as- 
sume that the unperturbed Hamiltonian consists 
of a part describing the free ion (say, Cr***), a 
part describing that portion of the crystalline e- 
lectric potential which is an even function of the 
coordinates, and a part describing the Zeeman 
energy arising from the exchange field. Ina 
material like Cr,O, the ground state is then non- 
degenerate. We further assume that the perturb- 
ing Hamiltonian is linear in the potential arising 
from E¢, linear in the potential arising from the 
constant part E° of the crystalline field, and 
quadratic in the spin-orbit energy AL-S. Thus 
the order of magnitude of the fourth-order en- 
ergy correction is given by 


W=(eE-r)(eE -r)(L-S)*/(a,4,4,), (5) 
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where A,, A,, 4,, are appropriate splittings 
calculable from the eigenvalues of the unper- 
turbed Hamiltonian. Alternatively, the quantity 


W=2 L gS-h (6) 


may be regarded as the Zeeman shift produced 
by h. From Eqs. (5) and (6), it follows that h 
is consistent with Eq. (2) provided we take 
eueE SP Lie pA 4549) (7) 
as an approximate measure of either of the (un- 
equal) questities ay anda,. The quantity r*® de- 
notes z* and xz in the persitel and poeennnns 
cases, respectively, and L? denotes Lz” and 
L,Lz in these two cases. It should be noted that 
the a’s resulting from the substitution of Eq. (7) 
into Eqs. (3) and (4) do fulfill the requirement 
of being invariant with respect to the product of 
the space inversion and time reversal transfor- 
mations, while being noninvariant with respect 
to either of these transformations separately. 
For a very rough numerical estimate applicable 
to Cr,0,, we use 4, ~A,=2x10* cm™, A,=1.6 
x 10° cm? (parity considerations), y=107* cm, 
4=10? cm", and eE,°z =1.7x10° cm“ (point 
charge calculations based on the structure of 
aFe,O,). The value’ a=0.6 107° resulting 
from Eq. (7) may be consistent with the values 
= 38.7x107% and a, = 1.32 x107° required to 
fit the experimental curves of a, anda, by 
means of Eqs. (3) and (4). More accurate es- 
timates based on actual calculations of the rele- 
vant matrix elements are under way. 
The writer wishes to thank V. J. Folen and 
R. D. Myers for helpful discussions. 
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DYNAMIC POLARIZATION ANOMALIES IN ORGANIC FREE RADICALS* 


Robert H. Webb 


Department of Physics, Stanford University, Stanford, California 
(Received May 12, 1961) 


As a part of a broader study of exchange effects 
in various paramagnetic substances, we have 
found an organic free radical in which the Over- 
hauser effect gives way to a “solid” effect at low 
temperature, and one in which the two effects 
are present simultaneously. 

The Overhauser effect! is produced by satura- 
tion of the electron spin resonance (ESR) transi- 
tion in a system in which the nuclei are coupled 
to the electrons by a hyperfine interaction. It 
appears as an enhancement of the nuclear polar- 
ization, and therefore of the nuclear magnetic 
resonance (NMR) signal. The “solid” effect,” 
on the other hand, depends on dipole-dipole 
coupling between the electrons and nuclei in 
a solid and requires the saturation of the for- 
bidden transitions at v= (ly, | +y,)H og, where v 
is the microwave frequency and y, and Y, are 
the electronic and nuclear gyromagnetic ratios. 
Ho is the external (dc) magnetic field. These 
forbidden transitions involve simultaneous flips 
of a nucleus and an electron. Saturation of the 
transition at (| Ye! +¥, Ho results in a negative 
enhancement of the NMR signal (the nuclear 
polarization is inverted), while that at (ly, -v,) 
xHg produces a positive enhancement. Typically, 
the forbidden transitions are more easily satu- 
rated than the allowed ones, so that much less 
saturating power is necessary to produce the 
solid effect than the Overhauser effect, if it 
is present. 

The Overhauser effect in the organic free rad- 
ical diphenyl picryl hydrazyl (DPPH) is relatively 
well known.* The enhancement obtained at 3400 


- oersteds is, however, small (about 20) at room 


temperature and at 77°K. It disappears in the 
liquid helium temperature range, being masked 
by a collapsing of the Knight shift and by micro- 
wave heating effects.* 

Wurster’s blue perchlorate (WB) has an ex- 
change-narrowed electron spin resonance about 
3 oersteds wide at room temperature.® At 300°K 
and at 77°K it gives only an Overhauser effect 
(Fig. 1). In the liquid helium range, it exhibits 
only a solid effect (Fig. 1). Here the ESR line 
is 14 oersteds wide and the proton NMR 14 oer- 
steds, giving rise to an effective broadening of 
the solid effect spectrum. Power dependence 
of the two effects is shown in Fig. 2, which is 


typical of Overhauser effects in all the free rad- 
icals, and of solid effects in general. 

Picryl amino carbazyl (PAC) resembles the 
other free radicals in its electron spectrum at 
room temperature.® It has one of the larger g- 
factor anisotropies among the organic free rad- 
icals. It exhibits both solid and Overhauser ef- 
fects at 300°K and 77°K, with a pure solid effect 
at 4.2°K and below (Fig. 3). 

These two free radicals differ from DPPH in 
the temperature dependence of their suscepti- 
bilities (see Fig. 2 of Rhodes et al.’). They also 
are known to show g=4 satellite’ lines below a 
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FIG. 1. Proton NMR signal in arbitrary units as a 
function of magnetic field in Wurster’s blue. The 
microwave frequency is held constant at 9.7 kMc/sec. 
ESR occurs at H)=0. 
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FIG, 2. Proton NMR signal in arbitrary units as a 
function of saturating power, at the peak of the en- 
hancements in WB. 


characteristic temperature T* which is higher 
in WB. The transition between the two types of 
dynamic polarization which we observe comes 
above T* in PAC, below it in WB. However, 
the existence of the g=4 lines is apparently 

due to the presence of dipolar coupling, which 
is in turn related to the decrease in suscepti- 
bility and its consequent effective magnetic 
dilution. It is this effective magnetic dilution 
which we suggest may provide the nuclear- 
electron dipolar coupling which gives rise to 
the solid effect. It does not follow that the tran- 
sition temperature for the nuclear-electron ef- 
fects need be the same as that for the purely 
electronic phenomenon. Support for this is 
provided by a fourth free radical, 1-3-bis- 
diphenylene-2-phenyl allyl (BDPA) which shows 
a g=4 line below 2°K,” but retains a reduced 
proton Knight shift® at this temperature—showing 
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FIG. 3. Proton NMR signal in arbitrary units as a 
function of magnetic field in PAC. The microwave 
frequency is held constant at 9.7 kMc/sec. ESR 
occurs at H)=0. 


that the proton-electron coupling is still pri- 
marily scalar (hyperfine). No solid effect has 
been observed in this substance at the lowest 
temperatures so far employed (about 1.8°K), 
and we have some evidence of‘a residual Over- 
hauser effect at 4.2°K. 

We are continuing our investigations of these 
materials in an effort to understand the nature 
of the susceptibility anomalies and the influence 
of exchange on the interactions among electrons 
and between electrons and nuclei in solids.°® 

The author wishes to thank Professor G. E. 
Pake for his help and interest in this work. Dis- 
cussions with R. S. Rhodes have also been inval- 
uable. 





*Supported by the Air Force Office of Scientific 
Research. 
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THERMAL CONDUCTIVITY OF NORMAL AND SUPERCONDUCTING LEAD ALLOYS* 


Peter Lindenfeld 
Rutgers University, New Brunswick, New Jersey 
(Received April 18, 1961) 


The thermal conductivity of lead specimens 
containing 3% indium, 6% indium, and 6% bis- 
muth has been measured in the normal and super- 
conducting states. The data are shown in Fig. 1. 
Their analysis shows that the addition of impurity 
leads to an increase in thermal resistance in both 
states, but that near 2°K the ratio of increase in 
the resistance to lattice heat transport is 40 
times greater in the normal state than in the 
superconducting state. Since in superconducting 
lead, at 2°K and below, the electronic conductiv- 
ity and the effect of the electrons on the lattice 
conductivity are expected to be negligible, it is 
seen that the large change in lattice thermal re- 
sistance with impurity in the normai state de- 
pends on the presence of the electrons. 

In the normal state the electronic thermal con- 
ductivity is equal to a7, where the coefficient 
a@ is related to the residual electrical resistivity 
Py by the Wiedemann-Franz law, a=L,/p, 

(L,= 2.445 107® volt?/deg?).! At sufficiently low 
temperatures the lattice conductivity is equal 

to 8T*. The total thermal conductivity K =aT 

+ 8T? is most easily separated into its compo- 
nents graphically, on a plot of K/T versus T 
which will then be a straight line. In Fig. 2 the 
experimental data for the normal state have 
been plotted on such a graph. For each speci- 
men the intercept a has been calculated from a 
measurement of the electrical resistance near 
4.2°K, anda straight line drawn through this 
point and through the low-temperature thermal 
conductivity points. The slope 6 of each line is 
shown in Table I. In addition Table I lists the 
values of p, for our specimens and of § and p, 
for a specimen of lead with 0.7% bismuth meas- 
ured by Montgomery.” The deviation of the data 
from the straight lines at higher temperatures 


is likely to be caused by the increased impor- 
tance of impurity scattering at these tempera- 
tures. 

The effect of adding impurity may be calcu- 
lated roughly by comparing the thermal resis- 
tances of specimens containing different amounts 
of impurity. The lattice resistances at 2°K are 
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FIG. 1. Thermal conductivity of lead alloys versus 
temperature. 
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FIG. 2. Normal state thermal conductivity divided 
by temperature versus temperature. 


listed in Table II. It is seen that the extra re- 
sistance of the 6% indium specimen compared 
to the 3% indium specimen is 190 cm-deg/watt 
in the normal state, and 4.9 cm-deg/watt in the 
superconducting state; i.e., the extra resistance 
is 40 times as large in the normal state as in 
the superconducting state. If the resistances of 
the 6% bismuth and 0.7% bismuth specimens are 
compared in the two states, the same ratio be- 
tween the extra resistances is obtained. Sup- 
porting evidence may be inferred from a com- 


Table I. Values of lattice conductivity coefficient 
B and residual resistivity pp. 


parison of Montgomery’s data with those of 
Olsen® for a specimen of lead with 10% bismuth. 

The results for the superconducting state are 
similar to those found in insulating crystals 
where boundary and point defect scattering pre- 
dominate. The difference between the effects of 
indium and bismuth may be at least partially ac- 
counted for by the difference in the mass of the 
impurity atoms. On the high-temperature side 
of the maximum, conduction and scattering by 
electrons will, of course, begin to affect the 
conductivity. 

That the lattice conductivity in normal metals 
decreased in an anomalously strong fashion with 
increasing impurity content has been known from 
work on other metals.’ The present measure- 
ments underscore our previous conclusion‘ that 
this change is not related to a large change in 
dislocation density since this would be expected 
to affect the conductivity in both states in the 
same way. The measurements show that the 
large increase in resistance with added impurity 
depends on the presence of the electrons, since 
it is greatly reduced in the superconducting state. 

This added resistance could, in principle, be 
caused by a change in the electronic density of 
states or the electron-phonon interaction. It 
seems more reasonable, however, that the 
presence of the electrons as scattering centers 
leads to an interference between impurity scat- 
tering and electron scattering similar to that 
discussed for a different case by Berman et al.° 
In that case the interference between isotope 
scattering and phonon scattering was such that 
the process of smaller magnitude turned out to 
dominate the thermal resistance. If a similar 
enhancement can occur for impurity scattering 
in metals, the anomalously large effect of im- 
purities on lattice conduction in lead as well as 
in other metals might be explained. 


Table II. Lattice resistance at 2°K in the normal 
state (1/Kgm) and in the superconducting state (1/Kgs). 











B Po 1/Kgn 1/Kgs 

(watt/cm-deg’) (ohm-cm) (cm-deg/watt) (cm-deg/watt) 
6% In 0.52x 1073 5.77 x10 6% In 480 13.7 
3% In 0. 86x 1073 2.33 x10~* 3% In 290 8.8 
6% Bi 0.77x 107% 6.09 x 1078 6 % Bi 320 4.4 
0.7% Bi® 1.3 x 107 0.750 10% 0.7% Bia 190 1.1 





a 
See reference 2. 
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a 
See reference 2. 
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The measurements were made with the help 
of J. S. Brown and W. B. Pennebaker. I would 
also like to acknowledge helpful conversations 
with Dr. P. G. Klemens and Dr. B. Serin. 
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DIRECT OBSERVATIONS OF ION DAMAGE IN CADMIUM 


P. B. Price 
General Electric Research Laboratory, Schenectady, New York 
(Received May 3, 1961) 


Pashley and Presland' have recently found that, 
during the normal operation of the Siemens elec- 
tron microscope, the filament gives off negative 
ions (of unknown origin) which are accelerated 
and strike the specimen. The measured flux 
density at the specimen ranged between 10’ and 
10° ions/cm? sec. They studied the effect of 
these ions on a number of metal films, partic- 
ularly gold,? and observed microstructural de- 
fects similar to those produced by fast neutron 
irradiation,* namely, tiny dot-like features which 
were sometimes resolvable as dislocation loops. 

This Letter reports evidence that ion damage 
in cadmium crystals irradiated in the electron 
microscope results in the nucleation and growth 
of vacancy-type dislocation loops. Dislocation- 
free platelets of cadmium, thin enough ($2000 A) 
to be transparent to 100-kv electrons, were 
grown from the vapor* and examined in trans- 
mission in the Siemens Elmiskop I, on a cold 
stage, at temperatures in the range +25° to 
-100°C. After a few minutes at high beam cur- 
rent, tiny, isolated dislocation loops appeared 
and began to grow in the region illuminated by 
the beam. They were first detectable as dots of 
about 40A diameter (the minimum size resolvable 
by the transmission technique). At low tempera- 
ture the rate of formation was high and the 
growth rate of individual loops was low, leading 
to a defect distribution like that in Fig. 1, which 
is a transmission electron micrograph, taken 
at -90°C about five minutes after the dots began 
to appear. At room temperature the rate of pro- 
duction was low, but those loops formed grew 
quite rapidly and eventually intersected one an- 





FIG. 1. Vacancy-type dislocation loops due to ion 
damage in cadmium at -90°C. Several of the smallest 
loops, clearly visible on the original micrograph, are 
indicated by arrows. 


other to form dislocation networks. Figure 2 
shows a typical distribution of loops (a) 5 min- 
utes and (b) 10 minutes after the beginning of 
irradiation of the same area at +25°C. 

The nature of the loops was determined by 
tilting a specimen and obtaining selected area 
diffraction patterns from individual loops in 
several different reflecting conditions. Virtually 
all of the loops were circular and oriented in the 
basal plane, which was parallel to the large 
platelet faces and perpendicular to the beam. 
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(a) (b) 


FIG, 2. (a) Vacancy-type dislocation loops due to 
ion damage at +25°C after 5 minutes of irradiation. 
(b) Same area after 10 minutes. Note the dislocation 
reactions which take place when expanding loops meet. 


The Burgers vector of all of the loops was found 
to be of the type $c +p, where p denotes a “half- 
dislocation” in the basal plane. It was inferred, 
from the observation that p could lie in any one 
of six equivalent (1010) directions, that each 
loop was formed by the collapse of a single layer 
of vacancies, together with a shear in the basal 
plane.® Each loop thus enclosed a single layer 
stacking fault which, when viewed in a 1010 re- 
flection, was visible because of the difference 
in contrast inside and outside of the loop. In 
Fig. 1 the faults are not visible. In Fig. 2, taken 
in a 1010 reflection, the loops have different 
degrees of contrast, indicating that they were 
formed at different distances from the surface 
of the crystal. The uniformity of the contrast 
inside each loop implies parallelism to the sur- 
face and therefore to the basal plane.® 

The following model is proposed to explain the 
nucleation and growth of the loops: 

The ions, which have an initial energy of 10° 
ev, behave in about the same way as the primary 
knock-ons of fast neutrons, losing their energy 
almost entirely by elastic collisions with cadmi- 
um atoms. Because of their small range, $200 
A,’ the damage is limited to a volume close to 
the surface of the specimen. It is estimated 
that about 200 Frenkel pairs are produced by 
each ion.* Of these, the interstitials have an 
extremely high mobility even at -100°C; the first 
few in the cascade are probably sputtered out of 
the crystal and the others diffuse rapidly to the 
surface, which in these thin crystals is an ex- 
cellent sink. 

This leaves about 200 vacancies, clustered in 
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the path of each ion, which may precipitate as 
a dislocation loop about 40 A in diameter.*»®,?° 
Because of its large curvature and line tension, 
such a small loop will be subjected to a very 
large radial stress tending to shrink it by the 
emission of vacancies, which will diffuse to the 
crystal surface. In the absence of any local 
vacancy supersaturation due to the presence of 
other clusters, it is calculated" that a 40A loop 
will disappear in ~10~° sec at +25°C, but only 
after ~10° sec at -75°C. A 100°C difference in 
irradiation temperature will therefore have a 
large influence on the kinetics of loop growth. 

At -75°C it is expected that a large proportion 
of the clusters will precipitate as stable loops, 
which will grow slowly as a result of the rather 
low mobility of vacancies arriving from nearby 
clusters. At room temperature fewer stable 
loops will be formed, due to the difficulty of 
maintaining a high local supersaturation, but 
these will grow rapidly because of the longer 
diffusion distance of vacancies. 

The transmission electron microscope tech- 
nique provides a means of studying ion damage 
directly and can therefore be used to study radia- 
tion hardening if stress-strain curves are taken 
of thin films before and after irradiation. It 
should be borne in mind by those who use the 
electron microscope to study dislocation distri- 
butions in thin specimens that point defects are 
created during observation and may, at suitable 
temperatures, create new dislocations (loops 
and/or tetrahedra’) and cause the existing dis- 
locations to climb. The latter effect has recently 
been observed in cadmium by the author,” in zinc 
by Fourdeux and Berghezan,’* and in magnesium 
by Lally.** 

Thanks are due to Dr. R. H. Pry, Dr. W. G. 
Johnston, Dr. R. M. Walker, and Dr. A. Howie 
for useful discussions, and to Dr. D. W. Pashley 
for communicating his results to me prior to 
publication. 
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DE HAAS-VAN ALPHEN EFFECT IN POTASSIUM* 


A. C. Thorsen and T. G. Berlincourt 
Atomics International, Division of North American Aviation, Canoga Park, California 
(Received April 28, 1961) 


A study of the alkali metals is of particular 
interest for an understanding of the electronic 
structure of metals since these metals should 
conform most closely to an idealized metallic 
model of free conduction electrons. Previous 
measurements of transport properties’ and spe- 
cific heats* of the alkali metals are in fair ac- 
cord with this model. In the free-electron pic- 
ture, a spherical Fermi surface in k space gives 
rise to de Haas—van Alphen oscillations whose 
period is given simply by P= (2e/hc)(3n?N)-*, 
where e, h, c have their usual meanings and N 
is the number of valence electrons per unit vol- 





FIG, 1. de Haas—van Alphen effect in potassium at 
1.77°K. The oscillating trace (~0.1-mv amplitude) 
shows the output from a pickup coil containing the sam- 
ple. The curved trace shows the field increasing from 
151.7 to 158.5 kilogauss during a sweep time of about 
1.0 millisecond (time increasing from right to left). 
The change in interval occupied by each successive 
period is due to the decrease of the time derivative of 
the magnetic field. 





ume. The period of the oscillations is related 
to the extremal cross-sectional area A of the 
Fermi surface perpendicular to the magnetic 
field by P=2zne/chA. For the case of potassium 
with 1.32 x10” electrons per cm‘ (at room tem- 
perature), one would expect a period of 5.69 

x 107° gauss~ corresponding to an extremal area 
of 1.68x10'® cm~*. Such a short period is of the 
same order of magnitude as those found in the 
noble metals,® and one would expect to observe 
such an effect utilizing the high-field pulse tech- 
nique.* 

In this paper we report on the observation of 
de Haas-—van Alphen oscillations in a single crys- 
tal of potassium in pulsed magnetic fields up to 
160 kilogauss. Figures 1 and 2 show the ob- 
served oscillations as taken directly from the 





FIG. 2. de Haas—van Alphen effect in potassium at 

1.07°K. The diagonal trace shows the field decreas- 
ing (from right to left) from 132.8 to 126.9 kilogauss 
during a sweep time of about 0.9 millisecond. 
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oscilloscope. The measurements yield a period 
of 5.75x107°+1% gauss~', indicating a cross- 
sectional area of the Fermi surface of 1.66 x 10** 
cm™~, in good agreement with that predicted 
from the free-electron model. It might be point- 
ed out that the theoretical value is only valid to 
approximately 2% since no correction for ther- 
mal contraction was made. Furthermore, be- 
cause the measurements were carried out for 
only one orientation nothing can be said as yet 
regarding deviations of the Fermi surface from 
sphericity. A plot of lIn(@/T) versus T, where 

a is the amplitude of the oscillations and T is 
the absolute temperature, gives a straight line 
with slope -277k/@H, where H = magnetic field 
and B=eh/m*c, from which we deduce the ef- 
fective mass of the electrons as (0.90+ 10%)m, 
(m,=free-electron mass). The specific heat 
data of Douglass et al.” give m* =1.3 my. 

A brief description of the pulsed-field appa- 
ratus used in the present experiment is given 
elsewhere.® The pickup coils used to detect the 
oscillations in magnetization of the specimen 
were wound on a 0.79-cm long form using No. 52 
copper wire. They consisted of 3050 turns wound 
series-opposed on top of 4800 turns, so that the 
effective sensitivity was proportional to the dif- 
ference in the number of turns. The pure metal 
was supplied by M. S. A. Research Corporation 
and had a ratio of room-temperature resistivity 
to helium-temperature resistivity of approximate- 


ly 1300. The crystal (0.019 cm diam by 0.81 cm 
long) was grown by zone-melting techniques in 
an open-ended Pyrex capillary and was protected 
from oxidation by a film of mineral oil. Attempts 
to determine the crystal orientation with x rays 
were unsuccessful due to the small size of the 
specimen and the small atomic scattering factor 
of potassium. 

The authors wish to acknowledge the sugges- 
tion of D. Shoenberg regarding the use of coaxial 
pickup coils and valuable discussions with J. Tri- 
visonno, Charles S. Smith, R. Bowers, and A. F. 
Kip concerning techniques for handling alkali 
metals. 
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OPTICAL SATURATION OF F-CENTER SPIN RESONANCE* 


B. R. McAvoy, D. W. Feldman, J. G. Castle, Jr., and R. W. Warren 


Westinghouse Research Laboratories, Pittsburgh, Pennsylvania 
(Received May 8, 1961) 


The influence of external light on the spin reso- 
nance of F centers in KCl has recently been ob- 
served,'*? with somewhat ambiguous results. 
This note reports results of preliminary experi- 
ments at 9.0 kMc/sec on very pure samples of 
KCl in which both the steady-state population 
difference and the spin relaxation time were 
measured during exposure to light. Irradiation 
by light in the F band was found to influence the 
populations of the ground-state sublevels as ex- 
pected by excitation to a bound state; lattice 
heating was found to be negligible. 

The ground state of F centers in KCl has been 
described in some detail by Gourary and Adrian.° 
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Relaxation‘ of the populations of the two levels 
of the ground state to the lattice temperature is 
accomplished by phonons. Absorption of light in 
the F band equally from both levels and the sub- 
sequent random return should drive the popula- 
tions toward equality and decrease the apparent 
relaxation time, t. The rate equations for n 
electrons in levels 1 and 2 may be written as 


dn ,/dt = -[2(w..+w2,)+2W+2W+AN]n, 
+ (wa, + W+W+AN,)n = -dn,/dt, (1) 


where®**® n=n,+n, with no storage in the optically 
excited states, w,, is the probability of transition 
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from 1 to 2 induced by the lattice, W is the proba- 
bility of transition induced by the microwave radi- 
ation, ‘“\ is the probability of transition from 1 to 
2 induced by the optical radiation, andA is the 
probability of transition by cross-relaxation to 

one of N impurity spins, of which N, spins are 
capable of flipping in pairs with F electrons in 
level 1. Since the N impurity spins are more 
strongly relaxed to the lattice than the F-center 
spins, N, is independent of mn, and the F -center 
population difference, An(t), is given by 


An(t) - An() = Bexp(-at), (2) 


where ) = 2(w,.+W2,)+2W+2W+AN=1/r7, An(~) 
=[w2,-Wy2+A(2N,-N)|/d, B is the initial pertur- 
bation, and 2(w,,+w,,)+AN=1/7, the reciprocal 
of the spin relaxation time observed with zero 
incident light. With low power input, we have 
found 2W« 1/1). Therefore, for low tempera- 
tures the expected relations between observed 
quantities are 


1/7 =(1/7.)+ 2, (3) 
An(~)/An(%),=T/To, (4) 


where An(~), is the steady-state population dif- 
ference without light. 

The samples were prepared from single crys- 
tals of pure KCl which had been grown in a hori- 
zontal furnace under an atmosphere of HCl. The 
method of purification will be described else- 
where.” The resultant purity was found by ionic 
conductivity to be about thirty parts per billion. 
The crystals were colored to saturation by brems- 
strahlung from 2-Mev electrons and kept in the 
dark until used. 

In the experiment, light from an Osram XBO 
2001 lamp is focused into a quartz prism mono- 
chromator. Light from the monochromator is 
directed through the multiple glass walls of a 
helium -nitrogen double Dewar into a microwave 
cavity of special design. Details of a similar 
cavity design are reported elsewhere.® The cav- 
ity is resonant in the TE,,, mode with the sample 
mounted at one end of the cavity near a window. 
The transmission of the window to light is ~ 0.8, 
and to microwave is ~10~*. Light transmission 
of the sample is monitored through a second 
window by a phototube. The power in the light 
beam incident on the Dewar is measured by a 
thermopile. The optical absorption spectrum 
of each sample is measured at 2.1°K. The abso- 
lute photon absorption rate is obtained from the 
measured absorbance of the sample with cor- 
rections for reflection losses at nine surfaces, 


transmission of the cavity window, and estimated 
scattering. The heat absorbed in the sample and 
cavity from the beam is found by the change in 
heater power required by the temperature regu- 
lator® to be in agreement with the thermopile 
values. The heat input (~20 microwatts into the 
sample) gives a negligible rise in temperature 

in the center of the sample whose outer surfaces 
are bathed in liquid helium II. A microwave sig- 
nal power of 10~® watt is employed together with 
a rapid field sweep of 100 psec duration twice 
each second. Resulting average microwave pow- 
er gives 2W much less than 1/7,. Resonance ab- 
sorption signals without saturation are observed 
in this way. Spin relaxation times are measured 
by the inversion-recovery technique®; for no 
light conditions, room light is excluded from the 
sample. 

Optical saturation is evident in Fig. 1 where 
the appearance of the resonance absorption sig- 
nals on the oscilloscope is given for the condi- 
tions with and without F-band light. Since the 
resonance line shape remains constant under 


t Light On 


SS 
= 





Light Off 


Detected i.f. Output, |Omv/cm 











Time, 100 psec/cm 


FIG. 1. Microwave resonance absorption by F 


centers in a crystal of KCl at 2.1°K and 3 kgauss 
with and without F-band light. 
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Table I. Effects of F light on ground-state properties of F centers in a y-irradiated sample of pure KCl at 








2.1°K. 
Photon absorption Observed spin Observed spin Calculated 
rate (10'* sec™') absorption relaxation time absorption 
signal (sec) signal 

(relative) [Eq. (4)] 

0 1.0 1350% 1.0 

1 0.5 540 0.4 

0.2 270 0.2 





8This rather long time appears to be the characteristic spin-lattice relaxation time and is under further study 


by Feldman, Warren, Castle, and Gourary. 


optical saturation, the signal height is sufficient 
for the comparison of Eq. (4). 

The observed resonance signals and relaxation 
times are given in Table I for one sample irradi- 
ated to a total dose of 5x 10° r and having 1 x 10*® 
F centers. Proportionality between the effective 
relaxation time and the steady-state population 
difference is found for two intensities of F light. 
This agreement confirms the model for optical 
saturation in which the excitation of isolated F 
centers occurs to these bound excited states and 
there is essentially no preferential pumping. 

In addition the data allow a comparison of ab- 
solute rates. For 1x10 photons absorbed per 
second and two photons absorbed for every op- 
tically induced spin flip, an effective relaxation 
time of 60 seconds is expected on the basis of 
Eq. (3) and 1x10" F centers. The number of 
F centers measured by spin resonance absorp- 
tion agrees with the number calculated from op- 
tical absorbance® within a factor of two. Uncer- 
tainties in the optical measurements probably 
account for the discrepancy from the observed 


tT of 270 seconds. 





*Work supported in part by a contract with the Air 
Force Research Division. 

1J, Lambe and J. Baker, in Quantum Electronics, 
edited by C. H. Townes (Columbia University Press, 
New York 1960), p. 93. 

"I, Wieder and J. S. Hyde, Bull. Am. Phys. Soc. 
5, 497 (1960). 

3B. S. Gourary and F. Adrian, Advances in Solid- 
State Physics, edited by F. Seitz and D. Turnbull 
(Academic Press, Inc., New York, 1960), Vol. 10, 
p. 127. 

‘Relaxation near room temperature for F centers 
in NaCl has been analyzed by W. E. Blumberg, Phys. 
Rev. 119, 1842 (1960). 

5N. Bloembergen, S. Shapiro, P. S. Pershan, and 
J. O. Artman, Phys. Rev. 114, 445 (1959). 

6J. G. Castle, P. F. Chester, and P. E. Wagner, 
Phys. Rev. 119, 953 (1960). 

"R. W. Warren (to be published). 

8B. R. McAvoy, Rev. Sci. Instr. (to be published). 

°J. H. Parker, Jr. (to be published); and R. H. 
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EXTREMELY LOW LOSS ACOUSTIC RESONANCE IN SINGLE-CRYSTAL GARNET SPHERES 


R. C. LeCraw, E. G. Spencer, and E. I. Gordon 
Bell Telephone Laboratories, Murray Hill, New Jersey 
(Received May 10, 1961) 


perature and frequency.’ 

The above observation was first made on a 
YIG sphere using magnetoacoustic resonance 
(MAR). With this technique a microwave pump 
signal is used to excite parametrically into 
acoustic oscillation a small spherical sample of 
YIG biased with a dc magnetic field near fer- 
romagnetic resonance.” At a critical rf mag- 


An acoustic resonance has been observed at 
room temperature and ~ 9.2 Mc/sec in single- 
crystal spheres of yttrium iron garnet (YIG) 
with losses nearly an order of magnitude lower 
than any other known material at the same tem- 
perature and frequency. Until now, quartz crys- 
tals vibrating in shear modes were recognized 
as having the lowest acoustic losses at this tem- 
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netic field, herjt, the sample oscillates in one 
of the acoustic modes of the sphere.* 

Using MAR the Q’s of a number of the acoustic 
modes of a YIG sphere have been obtained by ob- 
serving the decay time of the acoustic amplitude 
after the microwave pump has been turned off, 


@ being defined by the decay factor exp[-(w/2Q)é]. 


The modes S,, R,, and R,, shown in Fig. 1, have 
been studied in particular.‘ 

Although MAR was used for the original meas- 
urements, it was realized that a method permit- 
ting the selective excitation of a given acoustic 
mode was desirable. Thus the direct drive 
scheme shown in Fig. 2 was devised. Here the 
sphere is driven magnetostrictively at a fre- 
quency F/2, exciting a particular acoustic mode 
at a frequency F. The Q of the driven mode is 
then determined from the frequency difference 
AF between the points of one-half power trans- 
mission, using the relation Q=F/AF. It was 
established that the input and output coils con- 
tributed negligible loading. Acoustic oscillations 
have also been excited by driving at F. This is 
probably due to a “bias effect” in which the dc 
magnetization is not exactly parallel to the ap- 
plied dc magnetic field. The bias can be caused 
by several different kinds of slight sample im- 
perfections or bottom effects. In most cases 
the response when driving at F is greater than 
at F/2. At higher amplitudes, when driving at 
F, parametric excitation has been observed in 
which the measured @ depends strongly on the 
drive amplitude. 

The transmitted signal is monitored by a nar- 
row band superheterodyne receiver connected 





s - Ry 


FIG. 1. Schematic representations of three of the 
acoustic modes of a sphere. The letters S and R 
refer to spheroidal and radial modes, respectively. 
The subscripts refer to the number of zeros along a 
radius. The dashed lines indicate the extremes of the 
motions (omitted in R,). The solid curves show the 
relative amplitudes along the radius (for @=0 in S;). 
Both R; and R, are radially symmetric, R, having a 

breathing” type motion, and R, having an inner sphere 
and an outer shell with opposite radial displacements. 








to the output coil. The response is square law 
when driving at F/2 and linear when driving at 
F. The input frequency must be highly stable, 
since the linewidths are of the order of one 
part in 10’ of the resonant frequency. It may be 
noted that the output coupling is primarily 
through near-field effects, since mode R,, for 
example, has no oscillating dipole moment. 

The results obtained using direct drive and 
measuring transmission vs frequency are in 
good agreement with those obtained using MAR 
and observing the exponential decay. These re- 
sults are tabulated in Table I for the sample 
both in air and in vacuum and are compared 
with the highest published Q’s for quartz! and 
silicon.® 

At present the physical implications of these 
data are somewhat obscured by the fact that con- 
tact losses at the bottom of the quartz tube have 
not yet been eliminated. A levitation scheme is 
being developed to remedy this. We should like, 
however, to make several suggestions concerning 
these yuite unexpected results: Most YIG crystals 
grown to date are known to contain relatively large 
numbers of imperfections such as dislocations,° 
small inclusions, tube-like voids,’ and chemical 
impurities. Hence by the usual criteria for lattice 
perfection, the crystals used in these experiments 
are quite poor. Thus, it is not surprising that 






EVACUATION HOSE 


QUARTZ TUBE 


OUTPUT COIL, FREQ. F 
( xy PLANE ) 








SPHERE F 
INPUT COIL ,FREQ. > 


( xz PLANE ) 


FIG. 2. Schematic diagram of the direct drive 
method. The sphere rests loosely on the bottom of 
the quartz tube and is oriented by the dc magnetic 
field along the [111] crystal axis. The principal pur- 
pose of the dc field is to saturate the sample and elimi- 
nate the domains. (Below saturation the observed Q’s 
are very low.) The input and output coils are adjusted 
to give minimum direct coupling when the sphere is not 
acoustically vibrating. 
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Table I. Comparison of the Q’s for several modes 
in a single-crystal YIG sphere 0.060 inch in diameter 
and the highest published Q’s for quartz and silicon. 
The frequencies are tabulated here to only three signifi- 
cant figures, and are quite close to those given for an 
isotropic sphere by Love.* The velocities used were 
measured on polycrystalline YIG by H. J. McSkimin 
and are vz=3.87 10° cm/sec and vj = 7.17% 10° cm/sec. 
Contact losses due to the quartz tube are contained in 
these data. Essentially the same results have been 
observed on more than one sample. 








Freq. 
(Mc/sec) Q (air) Q (vac) 
YIG (S;) 2.10 0.07x10° 1.2 x 10° 
YIG (R;) 3.97 0.26x10° 0.45 10° 
YIG (R,) 9.22 2.5 x10° 9.5 x10° 
Quartz 9.22 1.5 x 10° 
Silicon 9.22 0. 42x 10° 





4See reference 3. 


mode S,, which has both shear and compression- 
al components and is quite analogous to similar 
modes observed in quartz, has a lower @ than 
quartz at the same frequency.’ 


The radial modes R, and R,, however, are not 
analogous to any motion thus far observed in 
quartz or silicon. In these modes the crystal 
unit cell experiences no shearing stress; and 
the motion, to first order, is purely “hydro- 
static” in form. It is possible that this has a 
twofold effect: Damping due to several types 
of dislocations vanishes, or is greatly reduced, 
when no shear is present.* Also, even in a per- 
fect crystal, certain phonon-phonon relaxation 
terms may vanish or be considerably reduced 
when the motion is such that the unit cell remains 
cubic.* With respect to the latter point, it is 
worth noting that YIG has an unusually isotropic 
sound velocity, there being only about one percent 
variation in the longitudinal velocity along vari- 
ous crystal axes.’° This fact is probably quite 
significant, particularly in the radial modes 
where dislocation effects are less. Quartz, for 
example, due to its large elastic anisotropy would 
not remain nearly as radially symmetric as 
YIG under motion such as in the R modes. In- 
formation about the change in velocity with hy- 
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drostatic pressure in YIG would be very help- 
ful but is unavailable at present. 

We cannot explain why the Q for R,, the first- 
order radial mode, is so much lower than for 
R,, especially since its frequency is even lower 
than R, (the normal frequency dependence for Q 
above a few Mc/sec is Q~f~'). This large dif- 
ference has been observed on more than one 
sample, and must await levitation experiments 
before being resolved. Although contact losses 
are undoubtedly greater for R, than for R,, cal- 
culations show that the ratio of the surface ve- 
locity squared to the stored energy is only ap- 
proximately twice as great in R, as in R,. 

It is important to determine if the unusually 
large Q observed for mode R, is due primarily 
to the material or to the radial symmetry; and 
if it is the material, whether it is the whole 
class of garnets or YIG in particular. Detailed 
studies are in progress to achieve a greater 
understanding of this new acoustic material and 
little-investigated type of motion. 





'H, E. Bommel, W. P. Mason, and A. W. Warner, 
Jr., Phys. Rev. 99, 1894 (1955). 

*E. G. Spencer and R. C. LeCraw, Phys. Rev. Let- 
ters 1, 241 (1958); J. Appl. Phys. 30, 149S (1959). 

3A, E. H. Love, The Mathematical Theory of Elas- 
ticity (Dover Publications, New York, 1944), 4th ed., 
Chap. 12. 

‘In the notation of P. M. Morse and H. Feshbach, 
Methods of Theoretical Physics (McGraw-Hill Book 
Company, New York, 1953), p. 1873, the S and R 
mode amplitudes can be written 5 = L. min*+ AN min» 
where m <I and n is added to indicate the number of 
zeros along the radius. The function L is derivable 
from a scalar potential and N from a vector potential. 
Mode Sy corresponds to m=0 andl=1. Mode Ry cor- 
responds tol=0, for which Noon =0. 

5J, Lamb, M. Redwood, and Z. Shteinshleifer, Phys. 
Rev. Letters 3, 28 (1959). The value here is extrap- 
olated and is for longitudinal waves. These authors 
report in private communication a somewhat higher Q 
for shear waves. 

°T, H. Ramsey, Jr., J. Am. Ceram. Soc. 42, 645 
(1959). 

"J. Nielsen and J. F. Dillon, Jr. (private communi- 
cation). 

‘Ww. P. Mason, Physical Acoustics and the Properties 
of Solids (D. Van Nostrand Company, Inc., Princeton, 
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OBSERVATION OF FORBIDDEN RESONANCES IN OPTICALLY DRIVEN SPIN SYSTEMS 


William E. Bell and Arnold L. Bloom 


Instrument Division, Varian Associates, Palo Alto, California 
(Received May 1, 1961) 


In a recent Letter! we described experiments 
in which a macroscopic moment precessing at 
the Larmor frequency was induced in a spin sys- 
tem by the synchronous modulation of an incident 
beam of optical pumping radiation. It was pointed 
out there that this resonance effect does not in- 
volve a perturbation directly coupling the mag- 
netic substates; instead it may be regarded as 
an indirect effect employing the coupling of each 
of the magnetic sublevels to a common optically 
excited state. This statement carries the im- 
portant corollary that such a resonance might be 
observed even if a direct coupling between the 
sublevels does not exist, i.e., even if there is 
no possible matrix element for the direct tran- 
sition. We report here the observation of such 
a forbidden resonance in the Am =2 (g=4, 5.6 Mc 
sec gauss) transition between the sublevels 
m=1 and m=-1 in the 2°S metastable state of 
helium. It is clear that the method is generally 
applicable to other optically pumped systems 
having more than two sublevels. 

In order to describe the experiment, we bear 
in mind that the transformation properties of 
triplet state wave functions under rotations are 
the same regardless of whether the angular mo- 
mentum is spin or orbital. It is therefore con- 
venient to take advantage of the spatial repre- 
sentation afforded by the orbital functions, whose 
angular parts can be written as follows: 


m=1, 4,=(x+iy) exp(-iw,f), 
m=0, y=V2z, 
m=-1, p_,=(x -iy) exp(éw,f), 


the functions being suitably normalized on a 
sphere of radius ry. w, is the Larmor angular 
frequency in a field H, parallel to the z axis 

for the usual magnetic transition Am =1. If an 
ensemble of triplet helium metastables is ir- 
radiated by resonance radiation, then, as shown 
by Colegrove and Franken,’ optical pumping re- 
sults in an alignment process tending to populate 
the state described by V2z, the effect being a 
maximum when the axis of the light beam is in 
the z direction. (The axis of the light is either 
the direction of propagation of unpolarized light 
or the direction of the electric vector of plane- 
“ized light. For our purposes it is imma- 


terial whether the actual alignment results in an 
excess or a deficit in the specified state.) 
Consider now what happens if the light axis is 
rotated so that it is perpendicular to H,, specif- 
ically if it is in the x direction. Instantaneously 
the light pumps a small alignment in the state 
described by V2x. This alignment immediately 
proceeds through a time development given by 


V2 (x coswot +y sinw,f) = (v, + p.)/2, 


whose expectation value for population excess in 
the x direction varies as cos’w,¢ and is therefore 
periodic at twice the Larmor frequency. If the 
light is left on continuously, then the successive 
smearing-out of phases still results in a net 
alignment, but much weaker than the maximum.’ 
However, let the light be modulated at frequency 
2w,; then the result is a synchronous reinforce- 
ment of the time-dependent effect resulting ina 
relatively large precessing alignment. Bya 
treatment exactly analogous to that used in refer- 
ence 1, it can be shown that the observable effect 
should be an increase in intensity of the transmit- 
ted light when the light modulation frequency w is 
within about one linewidth of 2 w,. 

The possibility of obtaining double-frequency 
optical effects upon driving a spin-1 system at 
its Larmor frequency has been pointed out else- 
where.*® However, the unique feature of the 
present experiment is the fact that the time- 
dependent alignment defined by the state (y, +-,)/ 
v2 does not correspond to any kind of directly 
observable macroscopic moment in the case of 
helium. Furthermore, note that the m =0 level 
does not enter in the description of the resonant 
state. Thus the resonance should occur as de- 
scribed even if the m =0 level is displaced from 
its median position between the m =+1 levels for 
any reason. The resonance linewidth is expected 
to be about the same, in frequency units, as that 
of the Larmor resonance and therefore only half 
as large in terms of magnetic field changes. 

The experiment was performed with substan- 
tially the same apparatus described in reference 
1 but without the circular polarizer. Since align- 
ment effects in He are much smaller than the 
orientation obtained with a circular polarizer,* 
it was necessary to assure virtually 100% modu- 
lation of the light beam in order to make the sig- 
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nal detectable. This was accomplished most 
simply by operating in a weak field H,=0.01 
gauss and an observed Larmor frequency w,/27 
of 24 kc/sec. The g=4 resonance at 48 kc/sec 
was clearly observable on the oscilloscope with 
or without a linear polarizer in the light beam. 


Rotation of the linear polarizer verified, within 
signal-to-noise limitations, that the signal am- 
plitude obeyed the sin?@ dependence to be expec- 
ted by generalizing the above arguments (@ is the 
angle between light axis and H,). 

We are indebted to N. Steiger for assistance 


in performing the experiment and to R. G. Ald- 
rich, of the U. S. Naval Postgraduate School, 
for pioneering in the production of large signals 
by modulated light. 





‘w. E. Bell and A. L. Bloom, Phys. Rev. Letters 6, 
280 (1961). 

*F. D. Colegrove and P. A. Franken, Phys. Rev. 119, 
680 (1960). 

3A. Kastler, Compt. rend. (to be published). 

‘L.. D. Schearer, Second International Conference 
on Quantum Electronics (to be published). 





PION-PION INTERACTION IN PION PRODUCTION BY 1*-p COLLISIONS* 


D. Stonehill, C. Baltay, H. Courant, W. Fickinger, E. C. Fowler, 
H. Kraybill, J. Sandweiss, J. Sanford,f and H. Taft 


Yale University, New Haven, Connecticut and Brookhaven National Laboratory, Upton, New York 
(Received May 12, 1961) 


Since the first conjectures! that rise in the total 
7™~-p cross section between 300 and 600 Mev might 
be caused by a pion-pion interaction, this subject 
has received considerable attention. Theoretical 
analysis’ of high-energy electron scattering on 
protons and neutrons has predicted a resonance 
in the pion-pion interaction at a total di-pion en- 
ergy (w) of 4 to 5 pion masses, with isotopic spin 
and angular momentum both equal to one. Several 
analyses of ~-p experiments’ in the 1-Bev energy 
range have tended to confirm this prediction, and 
application* of the Chew-Low method has indicated 
a steep rise in the pion-pion cross section above 
w=4. Recent work® with 1.9-Bev 1~-p collisions 


shows a peak in the pion-pion interaction at w~5.5. 


We report here evidence of pion-pion interaction 
in 7*-p collisions at three scparate energies, 
which show striking effects attributable to a pion- 
pion resonance with w of about 5.5 pion masses. 
The data presented are results of a systematic 
study of pion-proton reactions at kinetic energies 
of 910 Mev, 1090 Mev, and 1260 Mev, which is 
still in progress. Photographs taken at the Cos- 
motron in the Brookhaven 20-inch hydrogen bub- 


ble chamber have been scanned for all interactions. 


All two-pronged collisions have been measured on 
_ a projection microscope of the Franckenstein type 
and have been processed by the Yale spatial re- 
construction and kinematical fitting programs on 
the IBM 704 of the New York University Computing 
Center. Each possible identification assigned by 
the computation has been compared by a physicist 
with other information (including ionization densi- 
ties) available from the photographs, to establish 
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the final identification. Cross sections for the 
various reactions, based upon the first compila- 
tion of these events, are shown in Table I. 

The influence of pion-pion interaction will ap- 
pear most readily in the single pion production 
processes: 1*+p—+-p+a*+m° and a*+p2-n+n* +n". 
Accordingly, the @ value of the two outgoing pions 
(that is, the kinetic energy of their relative motion 
in their mutual center of momentum) has been com- 
puted for each individual occurrence of single pion 
production. The distribution of Q values at each 
energy is shown in Fig. 1. This figure includes 
all identified events, without additional selection. 

In the reaction 1*+p+p+7*+7°, a definite peak 
appears at each energy in the Q-value region of 
400-500 Mev, extending well above the number of 
events to be expected from a uniform momentum- 
space distribution of secondary particles. The 
peaks also extend well above the distribution to 


Table I. x*-p cross sections (in mb) at 910-Mev, 
1090-Mev, and 1260-Mev kinetic energy. 





910 Mev 1090 Mev 1260 Mev 





40.3+2.2 
16.5+1.4 
11.9+1.2 


24.5241.3 
10.3+0.9 
10.4+0.9 

2.6+0.4 


30.1+1.6 

12.6+1.1 

10.8+1.0 

2.50.5 4.6+0.7 

1.3+0.3 3.940.6 6.9+0.9 
+0.018 


0.034 9 o12 0.25+0.02 0.42+0.07 


o 
total 


0 
elastic 
o(pa*n*) 


o(nn*n*) 


° multiple m prod, 
o(=-K) 
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FIG. 1. Distribution of pion-pion Q values (that is, 
kinetic energy of the two outgoing pions in their mutual 
center-of-momentum system) for the reactions 7*+p 
—p+nt+n° and 1'+p—-n+7'+n° at 910-Mev, 1090-Mev, 
and 1260-Mev laboratory kinetic energy of the incident 
pion. The curved lines are the Q distribution resulting 
from uniform distribution of the secondary particles in 
momentum space. The straight lines give the Q distri- 
bution resulting from isotropic decay of a pion-proton 
isobar of unique mass 1230 Mev. 


be expected from isotropic decay of the 3-3 pion- 
nucleon isobar, which is known to influence this 
reaction.® 

In the pion-pion Q distributions for 1*+p-n+7* 
+m*, on the other hand, no peaks appear in the re- 
gion 400 to 500 Mev. (There may be a peak for the 
1260-Mev pion energy in the vicinity of 350 Mev, 
but this is not statistically certain.) 

Now, two 7* mesons can only have a total isotop- 
ic spin of 2, whereas the 1*+-1° combination con- 
tains isotopic spin values of 1 and 2. Therefore, 

a resonance in the J=1 state can only appear in 





the pr*n® reaction. On the other hand, if charge 
independence is assumed, a resonance in the /=2 
state should appear 80% as nm*tn* reactions and 
20% as pr*n®. We conclude that the interaction 
responsible for the observed peaks is overwhelm- 
ing in the state J=1, and possibly not at all in the 
state J=2. This is consistent both with theoretical 
predictions!»? and with the conclusions of Walker 

At present, only an approximate estimate of the 
pion-pion resonance energy can be made from 
these data. At 910 Mev, the peak is attenuated 
and shifted to lower energies by the lack of avail- 
able momentum space at higher Q values. The 
peak at 1090 Mev also seems to be about 50 Mev 
lower than at 1260 Mev. However, the data are 
not inconsistent with a single resonance with w 
between 5.0 and 5.5, corresponding to a total di- 
pion mass of 700 to 770 Mev. The full width at 
half maximum appears to be ~90 Mev. (The rms 
error in the determination of Q values is less than 
10 Mev.) 

We are pleased to acknowledge the considerable 
help we have received from the regular Brookha- 
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EMISSION OF CARBON GROUP HEAVY NUCLEI FROM A 3+ SOLAR FLARE 


Herman Yagoda 
Air Force Cambridge Research Laboratories, Bedford, Massachusetts 


and 


Robert Filz and Katsura Fukui 
Emmanuel College, Boston, Massachusetts 
(Received May 8, 1961) 


Discoverer satellite XVII was launched into a 
polar orbit 2042 U.T. on November 12, 1960, a- 
bout 37 minutes after ground-level neutron moni- 
tors indicated their maximum cosmic-ray in- 
crease following a 3+ solar flare which com- 
menced at 1325 U.T. During the 50 hours be- 
tween launching and successful recovery, it is 
estimated that the outside of the capsule received 
a bombardment of ~2x 10° protons of energy E 
230 Mev.’ The instrument capsule, which car- 
ried a 772-gram block of Ilford G5 emulsion, 
had an apogee of 993+ 4 km and a perigee of 
188+3 km. The average latitude of apogee was 
20°S and the perigee occurred at 18°N. Of the 
31 orbits of 96.44-min duration, a total of 0.766 
day was spent above geomagnetic latitude 55°N 
and S. 

The 10x15 cm face of the emulsion block was 
placed adjacent to one wall of the aluminum cap- 
sule and radiation arriving perpendicular to this 
plane penetrated 2.2 g cm™ of condensed matter 
composed of light nuclei. This admitted carbon 
nuclei with energies in excess of 88 Mev per 
nucleon over a solid angle of 1.8 steradians. 
Owing to other instrumentation in the capsule, 
the rear of the block was shielded by approxi- 
mately 21 g cm~ of condensed matter. 

After several unsuccessful attempts at full 
development, which caused almost every grain 
of silver bromide to be reduced to metallic 
Silver, the bulk of the emulsion stack could be 
processed to a point where the preparations 
transmitted light (Fig. 1) by employing a very 
dilute amidol developer at near 0°C tempera- 
ture and adjusted by means of a sulfite-bisulfite 
buffer to a pH of 6.5. Details of the processing, 
of the flare exposure, and of further measure- 
ments on long-lived radioactive nuclides pro- 
duced in the metallic components are available 
in a preliminary report on the Discoverer XVII 
flight.?»* 

The final ability to discern the tracks of slow, 
highly charged nuclei above the dense, random, 
single-grain background is not due in its entirety 
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to the extremely weak development. Owing to the 
rapid diminution of the flare radiation following 
maximum, we estimate that about 50% of the tot- 
al dose was delivered during the first hour in or- 
bit. This intense radiation appears to have de- 
stroyed the basic sensitivity of the G5 emulsion, 
It is very likely that the hydrogen peroxide pro- 
duced by the radiochemical decomposition of the 
water retained by the gelatin played a significant 
role. It is known that hydrogen peroxide* has a 
powerful destructive action on the latent image 
produced by ionizing radiations. By diffusion 
into the silver bromide grains all surface latent 
images would be destroyed, but in the case of the 
largest or most sensitive grains a residual inter- 
nal latent image might persist if the nuclear par- 
ticle crossing the grain had a very high rate of 
energy loss. 

Prior to development, the flown emulsion and 
a control piece from the same batch were ex- 
posed to polonium alpha particles. The tracks 
of the 5-Mev rays were readily visible on the 





FIG. 1. 


Photomicrographs of three terminating 

flat tracks recorded in the Discoverer XVII emulsion 
block. The black bars define 10-micron intervals. 

The integral grain densities in the terminal 300 microns 
are as follows. A: 192 grains (Z =4 or 5); B: 308 
grains (Z =6 or 7); and C: 871 grains (Z ~20). 
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FIG. 2. Histogram of rela- 
tive abundance of particles 
based on a scale of grains in 
terminal 300 microns of track. 
The insert is not a part of the 
systematic survey and is in- 
cluded to show the gambit of 
the charge spectrum. 
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control sheet, but could not be discerned on the 
flare-irradiated recording medium. Examina- 
tion of the desensitized emulsion at 2000x magni- 
fication shows star-like structures of 3 to 6 
tracks originating from a common center. The 
grain density and range of these associated 
tracks clusters suggests that they are the spal- 
lation products of Z 2 3 commonly produced when 
Br, Ag, or I emulsion target nuclei are evapo- 
rated. 

Systematic examination of the first three 600- 
micron emulsion sheets facing the thin window 
at 600x magnification shows a group of tracks 
which enter the stack and terminate their range 
by ionization. For particles making small dip 
angles with the emulsion plane the grain density 
can be measured after suitable correction for 
coincident background grains. The histogram 
of the grain density in the terminal 300 microns 
of range (Fig. 2) shows a pronounced peak at 
275+ 25 grains. Most all of the light tracks have 


600 800 1000 


GRAINS IN TERMINAL 300 MICRONS OF TRACK 


the same sense of direction and make a spatial 
angle with respect to the longitudinal axis of the 
rocket of < 50°. 

If we plot the logarithm of the integral grain 
count N as a function of the logarithm of the re- 
sidual range R, tracks of varying grain density 
yield linear functions of essentially the same 
slope. This indicates that the strongly irradiated 
emulsion is still serving as a quantitative tool 
for measuring ionizaiion. It can be shown that 
a line of constant grain density intersects these 
functions at ranges R roughly inversely propor- 
tional to the charge Z. If we assign Z =4 to the 
tracks with an average count of 183+ 17 grains, 
then the peak at 275+ 25 grains yields a charge 
of 5.9. The charge assignments for heavier nu- 
clei, following this model, are indicated in Fig. 2. 

Table I is a summary of the frequency of heavy 
primary nuclei terminations as observed on Dis- 
coverers XVII, XVIII, and at balloon elevations. 
The Discoverer XVIII block was exposed during 








Table I. Comparison of heavy primary terminations in balloon and satellite exposures. 
Flight Date Altitude Z => 6 terminations 
(km) (apogee) (cc! day~') 
Bemidji 532 July 29, 1960 42.5 7.78 +1.15 
Bemidji 533 August 5, 1960 42 8.60 +1.04 : 
Discoverer XVIII December, 1960 638 8.65 +0.93 
Discoverer XVII November, 1960 993 900 +300 
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a relatively quiescent period of solar activity and 
could be developed normally. It exhibits essenti- 
ally the same frequency of terminating Z > 6 par- 
ticles as the balloons launched at 55°N, and the 
tracks enter isotropically. A comparison of the 
two satellite exposures indicates that during the 
3+ flare the flux of Z 2 6 particles capable of 
terminating in the detectors increased about 100- 
fold. The greater part of the 900+ 300 cc™ day 
. termination frequency is due to the Z = 6+ 2 group. 
If we assume that our low-magnification esti- 
mate® of 13+ 3 cc~! day~! applies to the more 
conspicuous tracks of Z > 10, then the M/H ratio 
during a flare is ~ 70 as compared with a value 
of 3 for the galactic heavy primaries. The large 
uncertainty in the Z 2 6 population is not statisti- 
cal in character, but originates from the sam - 
pling position, which reaches a maximum in the 
first emulsion sheet facing the thin window. The 
over-all evidence indicates that a beam of parti- 
cles of Z =6+2 was emitted from the sun during 
the 3+ flare. 

This identification is consistent with the ob- 
servation by Fichtel and Guss* of a flux of car- 
bon, nitrogen, and oxygen nuclei associated with 
a flare of magnitude 3, detected by nuclear emul- 
sions flown to 130 km at 1480 U.T. on Septem- 


ber 3, 1960. On the other hand, Kurnosova and 
co-workers’ report an anomalously large flux 

of nuclei with Z 2 15 associated with a Class 17 
flare which occurred at 1127 U.T. on September 
12, 1959. Their observations are based on two 
Cerenkov counters carried on the second U.S.S.R. 
cosmic rocket which responded to heavy particles 
with energies exceeding 1.5 Bev/nucleon. During 
the 17-minute period when the flux of Z > 15 in- 
creased about twelvefold, the fluxes of particles 
with Z 2 2 and Z 2 5 remained essentially normal. 
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EVIDENCE FOR A 1-1 RESONANCE IN THE /=1, J=1 STATE* 


A. R. Erwin, R. March, W. D. Walker, and E. West 
Brookhaven National Laboratory, Upton, New York and University of Wisconsin, Madison, Wisconsin 
(Received May 11, 1961) 


Since the earliest data became available on pion 
production by pions, certain features have been 
quite clear. The main feature which is strongly 


exhibited above energies of 1 Bev is that collisions 
are preferred in which there is a small momentum 


transfer to the nucleon.’ This is shown by the nu- 
cleon angular distributions which are sharply 
peaked in the backward direction. These results 
suggest that large-impact-parameter collisions 
are important in such processes. The simplest 
process that could give rise to such collisions is 

a pion-pion collision with the target pion furnished 
in a virtual state by the nucleon. The quantitative 
aspects of such collisions have been discussed by 
a number of authors. Goebel, Chew and Low, and 
Salzman and Salzman? discussed means of extract- 
ing from the data the m-m cross section. 
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Holladay and Frazer and Fulco* deduced from 
electromagnetic data that indeed there must be a 
strong pion-pion interaction. In particular, Frazer 
and Fulco deduced that there probably was a reso- 
nance in the J=1, J=1 state. A qualitative set of 
1-p phase shifts in the 400-600 Mev‘ region were 
used by Bowcock et al.° to deduce an energy of 
about 660 Mev in the 7-n system for the resonance. 
The work of Pickup et al.° showed an indication of 
a peak in the 7-7 spectrum at an energy of about 
600 Mev. 

The present experiment was designed to explore 
the 1-m system up to an energy of about 1 Bev. 
The 1~ beam was produced by the external proton 
beam No. 1 at the Cosmotron. A suitable set of 
quadrupole and bending magnets focussed the pion 
beam on a Hevimet slit about 10 ft from the Adair- 
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FIG. 1. The angular distribution of the nucleons from 
the processes 1~+p—> 1~+n°+p and 1-+p—>a-+nt+n. 


Leipuner 14-in. H, bubble chamber. The pions 
were guided into the chamber by another bending 
magnet. The measured momentum was 1.89+ 0.07 
Bev/c. 

Events selected for measurement were taken in 
a fiducial volume of the chamber. The forward- 
going track was required to be at least 10 cm 
long. Measurements were made on a digitized 
system and the output was analyzed by use of an 
IBM-704. The events were analyzed by means 
of a program based on the “Guts” routine written 


by members of the Alvarez bubble chamber group. 


Figure 1 shows the combined angular distribu- 
tion for the nucleons from the two processes, 7~ 
+p+1~+7°+p and 1~+p+7~+2*+n, which appear 
to be identical within statistics. The results in- 
dicate a large number of events with small mo- 
mentum transfer to the nucleons. 

We concentrate our interest on those events 
with small momentum transfer since these events 
satisfy the qualitative criterion of being examples 
of r-n collisions. Somewhat arbitrarily, we cen- 
ter our attention on cases in which the momentum 
transfer to the nucleon is less than 400 Mev/c. 
Table I gives the ratios of the three possible final 
states m-1m, n-2°p, and 1°nm, assuming the 1-7 
scattering to be dominated, respectively, by the 
I=0, 1, 2 scattering states of the 7-1 system. 

The experimental results in the last column in- 
dicate a strong domination by /=1 state. For the 
I=1 state the basic 1-2 scattering cross sections 
O(n “1% 2~n°) and o(n-1* +2~2*) are equal. 


Table I. Ratios of final states. 








Experiment 
1=0 I=1 I=2 (4 <= 400 Mev/c) 
ann 2 2 2/9 1.7 +0.3 
n np 0 1 1 1 
r°'n'n 1 0 4/9 <0.25+0.25 





The nucleon four -momentum transfer spectrum 
seems to be in qualitative agreement with the 
theory for the process in which a 7 is knocked 
out of the cloud. Figure 2 shows ideograms for 
the mass spectrum of the di-pions for cases with 
4 <400 Mev/c and A>400 Mev/c, where A is the 
four-momentum transfer to the nucleon. The curve 
for 4 < 400 Mev/c clearly shows a peak at 765 
Mev/c. In the ideogram for A>400 Mev/c the 
peak is still present but seems to be smeared to 
higher values of the di-pion mass, m*. One wor- 
ries that diagrams other than the one involving 
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FIG. 2. The combined mass spectrum for the 17° 
and s~*+ system. The smooth curve is phase space 
as modified for the included momentum transfer and 
normalized to the number of events plotted. Events 
used in the upper distribution are not contained in the 
lower distribution. 
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one-pion exchange might be contributing to the 
observed peak in this m* spectrum. In particular, 
an important contribution at lower energies comes 
from a diagram in which one of the m’s rescatters 
off the nucleon and ends up in the 3-3 state with 
respect to the nucleon. If one restricts the data 
to cases with A < 400 Mev/c this diagram does 
not seem to be very important, but if one takes 
cases with A>400 Mev/c many cases consistent 
with rescattering are found. 

In order to deduce values of the 7-7 cross sec- 
tion, we use the formula? 


Ko , 
1-7 


do a 3f? A? m* 2 
dm*d&* a (A7+1)" qr 

where 6,_, is the mean of o(1-1°+2~7°) and 

o(n“1*++a-a*). In the above formula all momenta 


and energies are measured in units of pion masses. 


q,; =momentum of the incoming pion measured in 
units of the pion mass. K =momentum of the pions 
in the di-pion center-of-mass system. Then 


3 2 2 Amax 2 2 
0-2 (Va 5m* B._ 
t\q,,) 0 1-m A. ( +) (A +1) 
. min” 
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FIG. 3. The 7-12 cross section as deduced from cases 
with the four-momentum transfer less than 400 Mev/c. 
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The results of this calculation using the experi- 
mentally determined 50’s are shown in Fig. 3. 
The results indicate a peak in the neighborhood 
of 750 Mev with a width of 150-200 Mev, which 
is about 3/4 of what it would be (12 74’) for a 
resonance in the J=1, J=1 state. Since this 
cross section was determined off the energy 
shell, it is difficult to estimate the effect of the 
interference of other diagrams and also the effect 
of line broadening.’ Whether or not the other 
peak and the S-wave scattering indicated in Fig. 3 
are real will have to await better statistics for 
verification. 
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INELASTIC ELECTRON-DEUTERON SCATTERING 
AND THE ELECTROMAGNETIC STRUCTURE OF THE NEUTRON* 


Loyal Durand, III 


Brookhaven National Laboratory, Upton, New York 
(Received April 25, 1961; revised manuscript received May 11, 1961) 


Recent experiments on elastic electron-proton 
scattering performed at Stanford’ and Cornell? 
have provided detailed information on the be- 
havior of the Dirac (charge) and anomalous 
magnetic moment form factors Fy and Fo 
of the proton for electron 4-momentum trans- 
fers up to q?=36 f~*. Corresponding informa- 
tion on the form factors Fy, and Fo, of the neu- 
tron has been obtained from studies of the cross 
section d°a/(dQ,dE, ’) for the inelastic scattering 
of electrons from the deuteron, e+d~-e+n+p.*° 
Analysis of the recent Stanford data on this proc- 
ess, using the theoretical cross section calcula- 
ted by Jankus® with the modifications introduced 
by Hofstadter and others,’ leads to two interesting 
conclusions. First, Fo, and Fo,, which were 
thought on the basis of earlier work to be ap- 
proximately equal for g* < 10 f~,® were found 
to differ markedly for qg* >5 f~?,* reflecting a 
difference in the spatial distribution of the anom- 
alous magnetic moment in the two nucleons. Sec- 
ond, the Dirac form factor of the neutron F 
was found to differ from zero for q?25f~, this 
result providing the first evidence for an extended 
distribution of charge in the neutron.* The values 
of the form factors F in and Fo, obtained from 
Fig. 1 of Hofstadter and Herman‘ are given in 
Table I for representative values of gq’. The 
values of Fy, are also shown in Fig. 1. 

We wish to report in the present Letter on the 
marked changes in the neutron form factors which 
result from a reanalysis of the Stanford data on 
inelastic electron-deuteron scattering* using a 
theoretical cross section obtained from a com- 
pletely covariant theory.* This work, which ex- 
tends considerably the semirelativistic theory 
given in a previous paper,® is based on the tech- 
niques of dispersion relations. Calculation of 
the cross section d*0/(dQ,dE,') is especially 
simple in the region of the large peak which cor- 
responds to the quasi-elastic scattering of the 
electron by a single nucleon. The dominant terms 
in the transition amplitude then arise from single- 
nucleon pole terms which can be calculated exact- 
ly. Retention of these pole terms alone would 
correspond to the replacement of the complete 
deuteron wave function in the nonrelativistic 


theory by its asymptotic form. The largest re- 
maining terms arise at the quasi-elastic peak 
from anomalous regions in the dispersion rela- 
tions. The spectral functions in these regions 
are closely related to the short-range parts of 
the nonrelativistic deuteron wave function. The 
proximity of the single-nucleon poles to the phys- 
ical region for the transition results in a relative- 
ly weak dependence of the peak cross section on 
the high-mass regions of the spectral functions, 
or alternatively, on the short-range structure of 
the deuteron wave function. It is consequently a 
good approximation to use in the dispersion rela- 
tions spectral functions calculated directly from a 
nonrelativistic model for the wave function. The 
denominators in the dispersion relations may 
nevertheless be retained in covariant form; the 
usual ambiguities’»® with respect to the kinemat- 
ics and the relativistic treatment of the deuteron 
wave function are thereby eliminated. We may 
also remark that the nucleon form factors appear 


Table I. Changes in the neutron form factors F,,, and 
F»,, which result from changes in the theoretical peak 
value of d?0/ (dQ, dE e) « The form factors with the super- 
script J result from an analysis of the experimental cross 
sections based on the modified Jankus theory.» The 
form factors without superscripts are obtained when the 
data are reanalyzed using the present theory. The cal- 
culations were based on the method of intersecting ellip- 
ses.° The tabulated values of the form factors denoted 
by asterisks correspond to the point of closest approac 
in those cases in which the ellipses failed to intersect. 











q° (f-) +” Fn hs Fy 
5.1 0.16 0.09 0.76 0.75 
10.0 0.12 -0.09* 0.49 0.38" 
11.5 0.03 -0.11* 0.46 0.34* 
15.0 0.18 -0.08* 0.40 0.28" 
18.0 0.25 -0.02 0.40 0.28 
21.0 0.17 0.00 0.33 0.22 

pSee footnote 13. 


See reference 7. 
See reference 4. 
ae footnote 14. 
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FIG. 1. Values of the Dirac form factor F,, for the 
neutron obtained by analyzing the data of Hofstadter 
etal. (reference 4) by the method of intersecting el- 
lipses using the modified Jankus result for d*o/(dQ,dE ,') 
(open circles) and the improved result of the present 
paper (solid circles). For the points indicated by open 
and solid squares, “experimental” cross sections were 
constructed using the smoothed values of the nucleon 
form factors given in Table I and footnote 13. The 
points indicated by solid triangles represent cases for 
which no intersection was obtained unless allowance 
was made for possible experimental errors. The un- 
certainty in these points is very large. It is expected 
that the tabulated values of Fin for these points are 
too large in magnitude, F,,, ~ 0 being possible within 
the accuracy of the data. 


naturally in the present theory, and are evalu- 
ated on the mass shell in the dominant pole terms. 
Off-mass-shell corrections have been estimated 
in the remaining expressions, and are thought 
to be small. 

The cross section d’q/(dQ dE ,') is well repre- 
sented in the neighborhood of the quasi-elastic 
peak by the approximate expression: 


2 Ce 
d o/(dQ dE | )=o 


Mott™™P™(m/EIM(b, aN\G, +6,), 


(1) 





where 
Onsott =1a°E sin~*(3 6) cos?(3 @), (2) 
and 


G, =F ii +(¢°/4 m*)k "Fo 
+(q?/2m*)(F, + KF, tan?($@). (3) 


Here E, is the energy of the incident electron, 
@ is the electron scattering angle in the labora- 
tory system, m is the nucleon mass, Ky and k 
are the anomalous parts of the nucleon magnetic 
moments, and p is the momentum and E the en- 
ergy of either of the outgoing nucleons in their 
center-of-mass system. The effects on the 
scattering of the structure of the deuteron are 
contained entirely in the function M (>, q). 

The cross section given in Eq. (1) neglects the 
following: (1) the effects of interactions between 
the outgoing nucleons; (2) scattering involving 
the D-state component of the deuteron wave func- 
tion; (3) small corrections involving interference 
between scattering by the neutron and by the pro- 
ton (<0.5% error); (4) the effects on the functions 
G,, and G, of the initial Fermi momentum of the 
nucleons in the deuteron (<0.5% error). The ef- 
fects of final-state interactions have been calcu- 
lated® by the methods of reference 9 using ap- 
proximate final-state wave functions matched 
to the neutron-proton scattering phase shifts of 
Bryan, Signell, and Marshak,’° and are found to 
decrease the peak cross section by roughly 1.5- 
2.5%, depending on gq and @. We have adopted as 
a reasonable approximation to these effects a 
uniform 2% decrease of the peak cross section. 
The D-state component of the deuteron wave 
function scatters less efficiently than the S-state 
component, and the peak cross section is de- 
creased by an additional 3% for a large range 
of q when this effect is included. The function 
M(p, q) may be shown to reduce for q21 to the 
form (constant) /p? at the quasi-elastic peak. 
Evaluating the constant for a repulsive core 
wave function,® and incorporating the foregoing 
corrections, one then obtains for the peak value 
of d*o/(dQ,dE,') only 


[d*0/(dQ dE (4.57 x1075)(1 + 0.05) 


e ) eak e "Mott 


x (m?/pE)(G +G,). (4) 
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The factor (1+ 0.05) represents our estimate of 
the uncertainties in the theoretical cross section 
associated with the foregoing effects, with the 
contributions to the scattering of meson currents 
in the deuteron, and of off-mass-shell correc- 
tions to the form factors. This result is essen- 
tially equivalent in form to that given by Gold- 
berg!! as an approximation to his impulse-type 
calculations, but includes some additional cor- 
rections, and is furthermore free of the ambigu- 
ities attendant on the impulse approximation. 

The results obtained from Eq. (4) differ from 
the peak value of the modified Jankus cross sec- 
tion’ in three respects: (1) The value of the over- 
all numerical factor 4.57 107° in Eq. (4) includes 
the corrections discussed above; (2) Eq. (4) 
contains a factor (m/E) originating in the three- 
body phase space which is absent in the Jankus 
cross section, and decreases the peak value 
significantly for large values of p; (3) the func- 
tions G; are replaced in the modified Jankus 
cross section by functions G,’, 


G,'=G,-(q/2m*)F 2 + Of(q?/4m?)*}, 


The terms of order (q’/4m*)? in G;’ cannot be 
retained consistently in the nonrelativistic theory,° 
and should be omitted from the modified Jankus 
cross section. The appearance of the extra term 
-(¢°/2m*)F};° in G;’ may be traced to an incorrect 
modification of the Jankus theory intended to make 
the peak cross section reduce for large ¢’* to a 
sum of Rosenbluth cross sections’ for the scat- 
tering of electrons from free nucleons. The re- 
sulting errors in the theoretical value of d*0/ 
(dQ,dE,') are quite large, especially at small 
scattering angles. 

The effect of the foregoing results on the values 
of the neutron form factors Fj, and Fo, deter- 
mined from the Stanford electron scattering experi- 
ments, has been investigated by reanalyzing the 
experimental data of Hofstadter et al.* for rep- 
resentative values of g*. The calculations were 
based on the theoretical expression for the peak 
value of the inelastic deuteron cross section d*o/ 
(dQ,dE ,’) given in Eq. (4). The function G, which 
appears there was evaluated when possible using 
the Rosenbluth cross section for elastic electron- 
proton scattering,” 


(do/dQ ) 


2 p= “Mottll + 2E,/m) sin(40)"G,, (5) 


in conjunction with experimental values of the 
elastic scattering cross sections.’ This result 


was combined with the experimental electron- 
deuteron scattering cross sections of Hofstadter 
et al.*»’® using Eq. (4) to obtain the Rosenbluth 
factor G,, for the neutron. The values of Fy, and 
Fo, were determined by combining values of G,, 
corresponding to a fixed q’, but different energies 
and scattering angles, using the method of inter- 
secting ellipses.” For g*=10 f~? and 15 f~*, values 
of d*o/(dQ,dE,’) and (do/dQ¢), at properly matched 
energies and scattering angles were not available 
in the published data. In these cases we have con- 
structed “experimental” cross sections using the 
Rosenbluth and the modified Jankus theories with 
smoothed values of the nucleon form factors taken 
from references 1 and 4.* Corrected values of the 
neutron form factors were then obtained by reana- 
lyzing these “data” using the corrected theoretical 
cross section of Eq. (4). We note finally that the 
ellipses failed in three cases [q?=10, 11.5, and 
15 f~*] to intersect unless allowance was made for 
possible errors in the cross sections.“ In these 
cases, we have determined rough values of F in and 
Fo, by choosing as the “intersection” the midpoint 
of the region of closest approach. It is expected 
from an examination of the effects of possible 
experimental errors that the value of Fo, so 
determined is slightly too small, and that of Fiy 
considerably too negative. Fy =0 would be con- 
sistent with the data within the possible errors. 
The results of this analysis are summarized 
in Table I. The corrected values of Fj, are in 
addition plotted in Fig. 1 along with the values 
obtained by Hofstadter et al.* using the modified 
Jankus theory. The changes in F,, are quite 
striking. The results of Hofstadter et al.* sug- 
gested that this form factor was a positive in- 
creasing function for g’ > 5 f~*, with a value 
Fin®™ 0.2 for g=20f~*. The present results favor 
small or zero values of Fj, over the entire range 
of q’ which was considered. This result, if sub- 
stantiated by the analysis of newer experimental 
data,”® would indicate that the outer structure of 
the neutron is more nearly neutral than has been 
supposed.*»* The changes in Fo, are less strik- 
ing, but range up to 30% for large values of q’, 
and bring this form factor closer in magnitude 
to Fo ’ A detailed analysis of the relevant data 
will be necessary before the variation of the neu- 
tron form factors with increasing g° can be spec- 
ified with any assurance. We wish in particular 
to emphasize that the corrected values of Fj, and 
Fo, given in Table I are quite sensitive to changes 
in the cross sections d*a/(dQ,dE,’) and (do/dQ,) p? 
and may be altered when the results of recent ex- 
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the same scattering angles and energies given by Hof- 
stadter, Bumiller, and Croissiaux.’ For q’=5.1, the 
values of d’o/(dQ,dE e ) interpolated from the data of 
Sobottka® by Hofstadter, deVries, and Herman! were 
combined with proton cross sections calculated from 
the exponential proton model of E. E. Chambers and 
R. Hofstadter [Phys. Rev. 103, 1454 (1956)]. This 
model was used by Sobottka to normalize his deuteron 
cross sections.’ The cross sections for g?=10 f~? and 
15 f-? were calculated as described in the text, using 
smoothed values of the nucleon form factors taken from 
Fig. 1 of Hofstadter and Herman.‘ The values assumed 
for the neutron form factors F J and F/ resulting 
from an analysis using the J s theory are given in 
Table I. The values assumed for the proton form fac- 
tors are as follows: q’=10 f?, Fip =0.48, Fo = 0.39; 
q@’=15 f*, Fy, =0.43, Fo, =0.21. 

“The failure of the ellipses to intersect indicates 
that the input data are not consistent with the theo- 
retical form of d’0/(dQ,dE e’) given in Eq. (4). How- 
ever, intersections could be obtained by changing 
either the large-angle or the small-angle deuteron 
cross section d*0/(d2,dE,') by 10% or less. An error 
in the cross sections of this magnitude is consistent 
with the estimated accuracy of the experiments .'» 
Furthermore, such errors are undoubtedly present, 
as indicated by the deviation of the experimental val- 
ues of G, from the values calculated for the same q’ 
and @ using neutron form factors taken from the smooth 
curve drawn through the points in reference 4. The 
deviations are of both signs, and range up to about 
50% of G,, but are generally smaller. Since scatter- 
ing by the proton contributes on the order of 1-3 of 
the peak deuteron cross section, depending on q* and 
@, the corresponding errors in d’o/ (dQ2,dE,’) are re- 
duced by factors of 2-3, and lie within the expected 
range. Conversely, it is clear that values obtained 
for G,, hence, for the neutron form factors, can be 
changed significantly by rather small changes in the 
experimental values of d’0/(dQ,dE,’) and (do/dQ,),. 

‘’The Stanford group has recently extended its 
measurements of the inelastic electron-deuteron, and 
the elastic electron-proton scattering cross sections 
to a number of energies and scattering angles not pre- 
viously considered, with the largest value of i now 
25 f?. Preliminary analyses at several of the new 
points suggest that F,,, may indeed be greater than 
zero, but somewhat smaller in magnitude than indi- 
cated in reference 4. [R. Hofstadter (private commu- 
nication). The author is greatly indebted to Professor 
Hofstadter for a number of communications concerning 
the Stanford experiments.] The Cornell group has also 
extended considerably its measurements at larger 
values of g? [R. R. Wilson (private communication)] . 
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We have investigated the applicability of the of the energy-independent boundary condition em- 
modified boundary condition model’? to the ex- ployed to represent the shorter range interaction. 
tensive proton-proton scattering data inthe energy The precision fit we have obtained requires only 
range extending up to 350 Mev. The success ob- nine parameters plus the acceptance of the usual 


tained with the model using a field-theoretical po- fourth-order static “perturbation” theory results 
tential tail indicates at once the validity of the one- for the potential tail. 


and two-pion exchange potentials (TPEP) in their For r>r, (7, is the same in all states), we as- 
appropriate range, and the physical significance sume a potential of the form 
V=V,+V,, 





7 s 3Y\e’ 
Va(r)=ftuc? * - 2 [ae +5.( = oa ur’ 





2 2M -2 —— 
vain =ftue?|-(-a 2H oak, (2 ur) re 6 as lal -R,(ur) +#-FRA(UP) +5,.Ry(u7)] 
(1) 


where 


R ioe) = FFF [te ark, (ax) +(4 + 3) Kol2x) 
21 = 1 4 4 s.23.3 - 
+&(3 27.79] (3 +a+ 5) K,6) (2 +a+ Ke " > 


x x 
>. 2.2 1 
(2+3+3) K,«)+(2 + 3) Kt) e" . 


——s 


1 > > 
R,(x) = \f (5 + oy) K,(2x) + 4 K,22)| + 9&(3 - 27.7?) 


x 


2 
riin)=2\l(S+ =), (2x) + 3K, @x)|- -44-27-79[(4+ 5. x) K, (x) (2+3) K,(e)]e" *. 


f and f, are the dimensionless rationalized pseudovector coupling constant, and are approximately 
(0.08)” when evaluated from meson scattering. We distinguish the numerical values of f and f, in 
order to test separately the significance of V, and V,. y and M are the pion and nucleon masses, 
respectively. A variable “pair suppression” is provided by the parameter \. We shall call é the 
“ladder” parameter. The other operators and functions are standard. h=c=1. 
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For f,?=f* and §=1, V is the TMO potential.® 
For §=0 we have the BW potential.* Thus the 
terms proportional to — represent the contribu- 
tion of the iterated (ladder) second-order diagram 
to the fourth-order potential which has been the 
subject of so much controversy.° It should be 
noted that pair suppression (A<1) is a nonpertur- 
bative result (or at best comes from high-order 
terms) but is usually included in the perturbative 
potentials. No spin-orbit term has been included, 
as both the phenomenological evidence® and the 
theory’ indicate a very short range behavior. On 
the basis of earlier results’? the boundary radius 
r, is expected to include the spin-orbit region. 
The phase shifts and the corresponding experi- 
mental uncertainties, as listed in Table I, were 
taken from Breit et al.° and from effective-range 
parameters. Only 0, 95, 210, and 310 Mev were 
fitted in the expectation that the results for inter- 
mediate energies would interpolate well. 


Starting with Coulomb functions at r=8 f, the 
wave functions at each energy were numerically 
integrated inwards (on the IBM-709 at M.I.T.) 
to r,, for a set of “potential parameters” r,, /,’, 
A, and —. For some runs f? was fixed at 0.08; 
for others we set f? = "Ay and they were varied 
together. At 7,, the logarithmic derivative of 
the wave function F 7), = 7,|(1/~77,)d¥,yj,/ar|, -_ , 
was evaluated. For coupled states,’ F is a 2x2 
Hermitian matrix, 


c 
Fyg-1,1 "% 
c* 


“s "zens 


At zero energy the scattering length, a,, was 
used to determine the asymptotic 'S, wave func- 
tion, which was then integrated inward to obtain 
F 499(0 Mev). The singlet effective range Yeff (calc) 


Table I. Model parameters at minimum. 





f2=f,2=0.0837; r= 0.7011 f; A=1.00; —=0.50; M=11.10; (My=14) 





Energy 6 " 46 ” 
dls (Mev) (radians) (radians) F AF F 
000 0 b 1.02 0.06 
95 0.419 0.026 1.15 0.07 
210 0.087 0.026 1.05 0.08 
310 -0.175 0.044 1.12 0.14 1.08 
220 95 0. 068 0.012 1.63 4.43 
210 0.140 0.017 3.66 2.39 
310 0.192 0.017 3.25 0.94 3.14 
011 95 0. 209 0.044 -7.33 2.42 
210 -0.035 0.061 -7.14 2.06 
310 -0.209 0.035 -8.65 1.55 -7.82 
111 95 -0. 209 0.009 10.11 6.37 
210 -0. 364 0.017 17.34 12.36 
310 -0.480 0.026 15.21 9.70 13.37 
211 95 0.168 0.009 0.11 0.46 
210 0.311 0.017 -0.78 0.35 
J=2 310 0.291 0.017 -0.19 0.12 -0.26 
c 
coupling 95 -0.279 0.070 -0.14 0.17 
210 0.131 0.035 0.13 0.16 
310 -0.105 0.035 0.03 0.14 0.01 
231 95 -0. 003 0.008 -4.04 0.05 
210 0.003 0.026 -3.98 0.09 
310 0.017 0.017 -4.43 0.27 -4.06 
331 95 -0. 026 0.052 -3.04 115.74 
210 -0. 039 0.017 -3.34 0.73 
310 -0. 062 0.017 4.31 55. 29 -3.24 





,Nuclear Blatt-Biedenharn phase shifts. 


@y=-7.7 +0.1 f; regexp) = 2.67 +0. 03 f; ro¢e(calc) = 2.68 f. 


“For 6 read the nuclear Blatt-Biedenharn coupling parameter €,. 
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was then calculated directly from F,,.,. and the 
potential parameters. 

The F 5 8° obtained were in general energy 
dependent and ro¢¢(calc) did not agree with 
reff (exp). The potential parameters were then 
varied to make the /’s as energy-independent as 
possible and to bring 7vo¢¢(calc) into agreement 
witn regf(exp). To do this in a statistically sig- 
nificant way, utilizing the experimental uncer- 
tainties in the phase shifts 4677, and in the ef- 
fective-range parameters Aa, and Are¢¢(exp), the 
function 


ok 2 = 2 
F ng*) -F yng i 7 ogg CA1C) - 7 seep) 
AF Jis” AP ¢¢\xP) 





M= 
JIsE 


was calculated and then minimized (using a gra- 
dient search) with respect to the potential param - 
eters. 


F 7p y,©)/4F ©) Lgl AF 7, (E)] 


is the weighted average of the F 7), at all different 
energies. AF 315(E) = (@F 715/85, 715) 45J15; where 
the 4577, are the experimental uncertainties from 
Table I. Note that (aF /86)~ measures the sensi- 
tivity of 5to F. Those 6’s in the sum in (2) which 
are insensitive to F have a reduced relative im- 
portance. 

The most probable value of M, corresponding 
to 50% correlation probability, M,, is given by 
the number of phase shifts less the number of 
F jj, less the number of potential parameters. 
Counting the relevant states and energies, as 
seen in Table I, M,=14. M,<14 implies that if 
one puts F =F at all energies, the corresponding 
phase shifts will usually lie within their experi- 
mental limits. The word “usually” is meant in 
its appropriate statistical sense. 

Equation (2) neglects the effect of correlations 
between the phase-shift errors because they are 
difficult to take into account numerically. The 
effect of including the correlations in M is usu- 
ally to decrease its value. Thus a “good fit” re- 
mains a good fit while a “bad fit” may become a 
“good fit” if the correlations are very strong. 
Thus our criterion suffices for an acceptable 
potential but can only weakly reject one. However, 
if M >M, an unusual set of correlations would be 
required to make the potential acceptable. We 
intend to inquire further into the effect of corre- 
lations. 


With f*=0.08, a minimum of M = 10.73 was ob- 
tained for f,? = 0.0863, A=1.000, —=0.500, and 
v,=0.7011 f. With f?=f,?, a minimum of M =11.10 
was obtained for f* =/,? = 0.0837 and no change in 
the other parameters. The values of the F 77, 

F Jig» ASF zg, and regg(calc) corresponding to the 
latter case are given in Table I. The values for 
yr, and F,.. are very similar to those obtained 
previously for special forms of the modified 
boundary condition model.’ 

It is highly significant that the best fit is ob- 
tained for the coupling constant determined by 
pion-nucleon data,® f?=f,?=0.082. A variation 
of f,2=f? of about 0.002 would increase M by 
unity. Previous OPEP analysis’® had determined 
that f?~ 0.08 was required for the higher angular 
momenta. We see at once a similar result for the 
TPEP region and for lower angular momenta, but 
not independently of the model for the internal 
region. However, if this model was not a sig- 
nificant description of the scattering, it would 
be expected that M>M,. The parameters A and 
€ fall within the expected range. The sensitivity 
to A and £ is small however. Varying & between 
0 and 1 only increases M to 12. This is presum- 
ably due in part to the negligible difference be- 
tween TMO and BW in the singlet even state. 
Putting A = 0 while the other potential parameters 
remain unchanged makes M =58. If we then min- 
imize holding A=0, we obtain M = 24.8, but re- 
quire f?=f,7=0.093. This complete pair suppres- 
sion is not consistent with our model, especially 
if we require f? =f,” =0.08. 

There is then a strong statistical indication of 
the quantitative correctness of the fourth-order 
p-p static potential within limits of the order of 
the size of the “ladder” contribution. Conversely 
the energy-independent boundary condition is a 
statistically satisfactory representation of the 
inner region. If the fourth-order perturbation 
static potential is accepted on theoretical grounds 
for y>7,, then the parameters given by the F’s 
and 7, are sufficient to account for all the scat- 
tering data. This is a much smaller number of 
parameters than those inherent in previous pro- 
posed phenomenological potentials.® The fact 
that r,=0.7 f indicates that the TPEP approxi- 
mation breaks down at this distance and implies 
the onset of nonlocal or energy-dependent effects. 
In more conventional representations the short- 
range spin-orbit potential has the same implica- 
tion. 

For />3 the phase shifts are independent of 
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F jig except for a very narrow range which gives 
rise to a sharp resonance.“ Thus these phase 
shifts can in general be considered to depend only 
on TPEP. For example, a test of TPEP inde- 
pendent of the boundary condition can be made by 
examining the AF’s for high angular momenta at 
low energy.’* The extremely large value of AF,,, 
x (95 Mev) is such a verification of TPEP, and 
the decrease of AF... with increasing energy is 
also significant. 

The following calculations are now in progress: 
(1) The phase shifts for higher / are being evalu- 
ated with our potential to further test TPEP. 

(2) The phase shifts at interpolating energies 
are being calculated. (3) The nonperturbative 
potentials KMO"’ and Klein-McCormick™ are 
also being tried in place of Eq. (1). (4) The neu- 
tron-proton phase-shift data’® are under analysis 
employing the modified boundary condition model 
with the potential forms used in the p-p analysis. 
Here the triplet even potential depends strongly 
on the ladder ambiguity so that one may hope to 
obtain a better value of &. 

The above extensions of the analysis will be 
published elsewhere together with a more de- 
tailed discussion of the method and results of 
this paper. We conclude from the present anal- 
ysis that the field-theoretic fourth-order static 
potential combined with an energy-independent 
boundary condition is a simple and accurate, 
physically significant representation of the data. 

We are grateful to the IBM-709 staff and in 
particular to Miss M. Wirt. 
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CONSEQUENCES OF A SCALAR ZAn COUPLING 


Jeremy Bernstein*? and Reinhard Oehmet 


Institute for Theoretical Physics, University of Vienna, Vienna, Austria 
(Received May 2, 1961) 





It is the purpose of this note to point out some effect of this branch line is to increase the value 
new experimental and theoretical possibilities (3) of the coupling constant. Aside from the 12 
which follow from the assumption that the rela- and KN thresholds for s>(m A +m,)*, we have 
tive DA parity is odd." Using a crude dispersion- _ajso left out the circular branch line |s| =m ,? 
theoretic estimate of the scalar 2Azm-coupling -m,,”, which is due to the intermediate pion 
constant, we discuss the rates for the decay states of the reaction7+7—A+A. This neglect 
processes £* +e*+v+A° and Z°~y+y+A°, and seems to be reasonable, provided there are no 
we remark briefly upon some other consequences low-mass two-pion resonances with T =0. 
of a possibly strong scalar 2 Az coupling. In what follows we make the working assump- 

Let us begin with the evaluation of the constant tion that the scalar coupling constant is of the 
8S Aq" Because the ZA parity is odd, the = can order one, say Bsaq/4t =1.5. Asa first ap- 
be considered as a composite An system in an plication, we estimate the rate for the decay 
sstate.2 The s-wave amplitude for 7A scatter- 57 ~A°+e" +5 on the basis of several assump- 
ing, which is defined by f,(s) =q~*(s) sin6,(s) tions. 
xexpido(s), bare 7 (s)=(4s)~[s - (my +m,)°] We assume that the strangeness-conserving 
x[s-(m, -m,)*], has a pole at s=my*. The vector current is conserved. Therefore we may 
residue of this pole may be expressed in terms neglect the vector contribution to the decay rate 
Cag owe since the DA transition has AT#0. The vector 

Bong Mytm,-m* 1 part vanishes in the allowed approximation in 

fols) = 4n 2m m.--Ss accordance with the selection rule AT=0. If 
= = the vector current is not conserved, one would 
+other (nonpole) terms. (1) expect that its over-all contribution amounts to 


an increase of the total decay rate. 


In order to make a simple first approximation, The matrix element of the axial vector current 


let us assume that f,[s(q)] is a regular function may be written in the form (noting the assump- 
of the c.m. momentum q, except for the = pole tion of odd ZA parity), 

at q=q(s=m>,")=ik=i(2m,,qB)”. From the (AIP. 127) 

unitarity condition we find a 


-k +ig(s) = aly a(-P*) - ip, b(-p") +0, .b,0(-P*) Ju» 


f(s) = Fisye (2) 


If we compute the residue’ of the pole at s=ms” 
from Eq. (2) and compare the result with Eq. (1), 


(5) 
where p =P -Py» and a(0)=G,’ is the correspond- 





we obtain‘ ing renormalized axial coupling constant. The 
Exae’ 16m,*x coefficients in Eq. (5) are real analytic functions 
ms 3 : : : x71 1.8. which are regular in the complex p* plane except 
4s Lm, +m.) —. Im," - (m,"-m_") J for the branch lines® -p? > (3 m,)’ and, in the case 
(3) of b(-p?), the pion pole at -p*=m,”. We assume 
that in the neighborhood of the interval 0< -p? 
How good is this “zero-range” approximation? < (my -m,)° the function (m5 -m ,)a(- ?) - p?b(-p?) 
Presumably the most important correction is is dominated by this pion pole. By well-known 
due to the = pole in the crossed reaction, the methods we find the relation’ 
— of which is again determined by gy,,- 
or the s-wave amplitude f,(s), the relevant 
contribution from this a ‘appears in the form (1N2)64'= 85 Flom my -m 4); (6) 
5 
of a short logarithmic branch cut, where F(m,’) is determined by the pion lifetime. 
(m A - m_*)*/m,*< s<2m Oy +2 m_* - m,”. (4) We may compare Eq. (6) with the Goldberger - 
Treiman formula for nuclear 8 decay. Denoting 
Preliminary estimates seem to indicate that the the usual axial vector coupling constant by G A? 
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we can write 
G4'/G 4 = (85 qqN2 Byy [2 my/ im, -m 4)]- 


According to our zero-range calculation we may 

take 8sAq /ENNa * 1/10 and find G,’/G,=6. We 
assume that in the region of interest the function 

a(-p*) is slowly varying so that a(-p*)=G,’; then 

we obtain 


b(-p*) =G,'(m,, -m )/(m_* +"). (8) 


In the computation of the decay rate for ~ + A° 
+e~+D, we can neglect the 6 term in Eq. (5) be- 
cause its contribution becomes proportional to the 
electron mass. The c term in Eq. (5) is difficult 
to estimate; it is unlikely to make a very large 
contribution, and we neglect it. With these as- 
sumptions, the rate is given in terms of our con- 
stant G,’ by the “allowed” B-decay formula, 
Po pose 47 Og /00m ims -m,). (9) 
We compare this decay rate with the experimental 
rate for the strangeness-changing process 2° -n 
+e@~+v, which is taken to be about 1/10 of the rate 
obtained from a straightforward calculation using 
the universal constant. With our value G4’=s6G,4 
weflei l.-  .0,..°*.e/t oo - -#1/10, a 
number which 4 about 30 fo 40 times larger than 
one would predict on the basis of phase-space con- 
siderations. This larger rate is a reflection of 
the renormalization effects due to the rather large 
scalar =An coupling. We believe that, as far as 
our approximations are concerned, we have con- 
sistently underestimated this effect. We note that 
the states |=~A°) and |Z*A°) are related by G con- 
jugation. If the current P,, is such that it trans- 
forms like #y,y5p under G conjugation, then we 
have G, ‘(=~ ~M+e7+7) = -G4'(Z ~A°+e +v). 
Consequently, within our approximations, the rate 
for the decay £* + A°+e*+v should be the same as 
that for the 2° decay. 

We hope that there will soon be enough data on 
8 decays of = hyperons to make possible a serious 
test of our considerations.® 

It has been mentioned by Okun and Rudik® that 
for odd DA parity an intermediate 7° state may 
play an important role in the calculation of the 
process 5°+A°+y+y. Since we propose a rather 
large value for the scalar ZAn coupling, it seems 
reasonable to estimate the rate on the basis of this 
contribution. We have 


= 2 
Sea evey” (4/945)(1/m)(g,, */4a)E oo 


2y 


+ magnetic terms, 
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and with a 7°-decay rate’°of 0.3 x 10** sec™ we see 

that, even for 25 a /4n #10, the rate for the decay 
r°+A°+y+y is negligible compared to that for the 

decay £°+A°+y, which may be of the order of 10° 
sec™’. The “magnetic terms” have been discussed 
in reference 9; they do not change the situation. 

A strong, scalar ZAm coupling gives rise toa 
corresponding spin-independent two-pion exchange 
force between A and N. Spin-dependent forces 
must have a somewhat shorter range because they 
require the exchange of three pions or a K meson. 

Finally we would like to remark that a large 
value of g>,,_, may well give rise to a J=1/2 
resonance in the 7A system." 

It is a pleasure to thank Professor R. Karplus 
and Professor W. E. Thirring for stimulating 
conversations. 





*Permanent address: Brookhaven National Laboratory, 
Upton, New York. 

tNational Science Foundation Post-Doctoral Fellow. 

tPermanent address: Enrico Fermi Institute for 
Nuclear Studies, University of Chicago, Chicago, 
Illinois. 

‘For a discussion of arguments in favor of odd EA 
parity, as well as for references to the relevant papers, 
see Y. Nambu and J. J. Sakurai, Phys. Rev. Letters 
6, 377 and 506(E) (1961). See also S. Barshay and 
M. Schwartz, Phys. Rev. Letters 4, 618 (1960). 

*§, Barshay and M. Schwartz (reference 1). 

3M. L. Goldberger, Y. Nambu, and R. Oehme, Ann. 
Phys. 2, 226 (1957); this paper contains an analogous 
determination of the residue of the deuteron pole in 
np scattering. 

‘A similar determination of gy, has been made by 
Nambu and Sakurai, reference 1. Essentially, these 
authors evaluate the pole term in Eq. (1) at threshold 
and equate it to -x~', It is easily seen from Eq. (2) that 
the extrapolation to q=0 gives rise to a factor } in the 
expression for gy, ,”/4 7. 

5It should be noted that the branch line (4) corresponds 
to a force with maximal range 7) = (2m,)-'. This is 
because, for systems of two particles with unequal mass, 
the “force range’ should be determined from the position 
of the left-hand singularities in the g? plane. The branch 
line (4) lies inside the circle |s|=m,?-m,?; it is mapped 
into the second sheet of the q? plane and determines the 
weight function along the cut -m,*<q?<-m,? [see R. 
Oehme, Phys. Rev. Letters 4, 246 (1960)]. The upper 
limit of this branch line corresponds to a range 7 
=(2m,)"'. 

‘We have (0|3,P,|27)=0 if we assume that 99Pq 
transforms like the pion field under G conjugation. 

Note that this requires that P, contains the ZA terms 
in the combination Kygt* -=Iy,A. 

"Equation (6) and related relations have been written 
down by many authors. See, for instance, L. Okun, 

J. Exptl. Theoret. Phys. (U.S.S.R.) 39, 214 (1960) 
(translation: Soviet Phys. -JETP 12, 154 (1961)]; 
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note to use Eq. (6) in order to produce a number for OR. F. Blackie etal., Phys. Rev. Letters 5, 384 
the decay rate. (1960). 
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POSSIBLE MECHANISM FOR THE PION-NUCLEON SECOND RESONANCE 


Ronald F. Peierls* 


Department of Mathematical Physics, The University, Birmingham, England 
(Received February 17, 1961) 


Several phenomenological models! have been angle q, @: 
proposed to discuss pion-nucleon interactions s=(E+w)?+iA(1+w/E), 
near the higher resonances. There seems to be 
considerable evidence from these analyses that S=(E +w)? -24°(1+cos@) +iA(1 -w/E), 
final states described by the 33 isobar model? t = -2.q7(1 -cosé), 
play an important role, but so far there has been 
little quantitative discussion or explanation of the $+5+t=2M*+2 ue +2iA, (4) 


resonant peaks. In the present note it is shown 
that a simple model can give a good quantitative 
description of the second resonance, a rough 
description of the third, and a plausible explana- 
tion of the peak above 1.0 Bev. The model con- 
sists in taking seriously the concept of the 33 
isobar as an unstable particle (N*) of spin and 
isospin 3 and having a complex mass. This means 
that we assume the low-energy pion-nucleon scat- 
tering to be dominated by the one N* direct pole. 
The position and width of the 33 resonance then 
enable us to find the mass of the isobar M(N*): 








[M(N*) P =M? +i, 
M=8.94u , |All =8.05p ?, (1) 
7 T 


where Le is the pion mass. 
We consider next the reaction 


1+N—7+N*. (2) 
This is strongly coupled to the process 
1+N* -7+N*. (3) 


Let us try to apply the usual formalism to reac- 
tion (3), pion-isobar scattering. For the moment 
we ignore the spin and isospin and assume that the 





reaction is described by some amplitude A(s, 5, t) FIG. 1. (a) Crossed one-nucleon diagram for N* 
where s, S, ¢ are defined in terms of the pion scattering. (b) Schematic coupling between 1N* scat- 
energy w, the real part E of the isobar energy, tering and sN* production in 7N collisions. (c) Sche- 


and the center-of-mass momentum and scattering matic coupling between sN* production and 7N scattering. 
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The amplitude A should have a pole whenever s, 
S, or tis equal to the mass of a possible inter- 
mediate particle. The term represented by Fig. 
1(a) has the unusual feature that if A=0, it be- 
comes infinite in the physical region for g,@. The 
explicit contribution of this term may be taken 
empirically to be 


A,=8°/(5-m?’), (5) 





where m is the nucleon mass and g is the “coupling 


constant” for the NN*z vertex® (determined by the 
n-N cross section at the 33 resonance). We can 
rewrite this as 





- & 
A,(W, @)= (w? - Ww.) +iA(1 -w/E) -2 ql -cos@)’ 
W,?=2M?+2 uw? =m’. (6) 


The resulting total transition probability is 


+1 

t(W) -{ |A,|?d cosé 
= g tan™ ( 4q°4, 

2q°A, A,°(W* - W,°)(W? - W,” -4q"))’ 





A, =1A1(1-w/E). (7) 


- 40 


- 20 


-10 








hw 
1400 sa 


This is a sharply energy-dependent effect and 
behaves like a resonance. This interaction is 
coupled to the mN inelastic scattering [reaction 
(2)] as indicated in Fig. 1(b) and hence, by uni- 
tarity, to the elastic scattering, Fig. 1(c). If the 
effect is sufficiently strong, it should produce an 
appreciable contribution to the total 7N cross sec- 
tion at pion energies T - corresponding to the 
center-of-mass energies W where f(W) is large. 
Figure 2 shows f plotted as a function of te and 
compared with the experimental results on 1N 
scattering.* The agreement with the position and 
width of the second resonance (N**) in 1-N scat- 
tering seems very good, in view of the fact that 
there are no free parameters in the model. 

Let us consider the effects of spin and isospin. 
In the energy range in question, 2q*<|Al; hence 
we may, as a rough approximation, ignore the @ 
dependence of A,(W, @) in (6). In other words, the 
amplitude should contribute mainly to the s state. 
The only effect of taking account of the spin is to 
introduce several amplitudes A ; each having some 
angular factor y;(9) in the numerator. But if the 
8 dependence of the denominator is not important, 
this results in multiplying f(W) by some weakly 
energy-dependent factor, and does not change the 
main features of the result. Thus, we can say 
that this effect, looking like a “resonance” in the 


FIG. 2. The function 
f(W) plotted as a func- 
tion of the pion lab en- 
ergy T, required to 
produce a total energy 
W in the center-of- 
mass system. The ap- 
proximate positions of 
the sN** and 7N*** 
resonances are also 
shown and these theo- 
retical predictions are 
compared with the ex- 
perimental results for 
the total cross sections 
for ™N scattering. 
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s-wave 1-N* scattering, should produce a strong 
effect in the (3 -) state of the 7-nucleon scatter - 

ing via the reaction of Fig. 1(b).° The isotopic 

spin dependence is not so clear; however, the 
contribution of A, to the T= 4 state and T= 

state have opposite signs and the interference of 
these two contributions with the (appreciable) back- 
ground may be constructive and destructive, re- 
spectively. This might even lead to a dip, rather 
than a peak, in the 7*p cross section. 

The physical interpretation of the mechanism 
discussed above is easily understood. The isobar, 
being unstable, can decay into a pion and a nucleon 
which subsequently absorbs another pion. This 
[Fig. 1(a)] represents a real process (apart from 
the finite isobar lifetime) and it is not surprising 
that it should show up in the physical region. The 
same phenomenon should be possible whenever an 
unstable isobar interacts with one of its decay 
products and may be the explanation of the “reso- 
nance” in KN scattering due to the possible ine- 
lastic excitation of a 1Y* state (Y* being the 7A 
isobar whose mass lies below the KN threshold).® 

The mechanism described corresponds to an s- 
wave 7N* interaction. The p wave can be studied 
in the static limit by the conventional effective- 
range theory which leads to predicted resonances 
in the (§ 3+) and (3 $+) states. These are the 
states suggested by Wong and Ross,’ who con- 
sider essentially the iteration of Fig. 1(a), but 
carry out the calculation in the static-nucleon, 
rather than static-isobar limit. (The mechanism 
discussed here does not occur in either static 
limit.) The third 7-N “resonance” (N***) near 
900 Mev can be interpreted as the manifestation 
of these p-wave resonances. 

Lastly, the arguments leading to the prediction 
of N** can be applied to possible 1N** and 1N*** 
interactions, with intermediate single-nucleon or 
N* states, giving rise to four expected peaks in r- 
N interactions in the energy range indicated in 
Fig. 2. The situation is much more complex in 


this energy region, but this may be a partial ex- 
planation of the ohserved broad peak in the 7N 
cross section above 1 Bev. 

The model presented here is only approximate; 
it would be desirable to study it more quantitative- 
ly, in particular taking proper account of the 
spins and of the coupling between the various 
channels. Further work along these lines is being 
carried out and will be reported later. 

Iam very grateful to J. E. Paton for many val- 
uable discussions during the course of this work 
and wish to acknowledge many stimulating com- 
ments and suggestions from L. Castillejo, D. J. 
Thouless, Professor R. E. Peierls, and Profes- 
sor S. Mandelstam. 
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SCATTERING MATRIX APPROACH TO THE A-7 RESONANCE* 


K. C. Wali, T. Fulton, and G. Feldman 
The Johns Hopkins University, Baltimore, Maryland 
(Received May 1, 1961) 


Recently Dalitz' interpreted the A-7 resonance? 
as a bound state of the K-N system. His method 
made use of concepts familiar in low-energy nu- 
clear physics. He was able to obtain the position 
and width of the resonance independent of the rel- 
ative 2-A parity. 

We attempt, here, to analyze the data involving 
the three channels KN, =z, Az, in the J=1 state 
in the total energy range m,j +m, to my +my by 
means of an effective-range theory derived from 
a dispersion approach (analogous to the Chew-Low 
approximation in the one-channel case). If we as- 
sume that the ©-A relative parity is odd, and A-K 
odd, we find that there can exist a resonance in 
the J=1, J=(4-) state and we can fit all of the 
known data to a reasonable set of renormalized 
coupling constants. A similar approach has been 
used by Amati et al.* for the case of even Z-A 
parity. They include only the 7-A and 7-2 chan- 
nels and establish the possibility of a resonance 
in the J=(3+) state. Since the actual position of 
the resonance is closer (but below) the KN thresh- 
old than it is to the Yn threshold, one might ex- 
pect (as Dalitz has observed) that the KN channel 
is important. 

We next present a brief outline of our method 
of calculation. 

It was shown by Matthews and Salam® and by 
Feldman, Matthews, and Salam* that by confining 
oneself only to the two-particle channels of the 
S matrix, one can express the unitarity condition 
nonrelativistically as 


-1 oe 
Im{T | dig ko 


3 


if’ (1) 
where T is the inverse of the open-channel T 
matrix defined in MS,’ and k, is the c.m. momen- 
tum in the 7th channel. The analytic continuation 
of T _" below threshold of one of the channels is 
easily accomplished.® Dispersion relations for 
the T ey matrix elements were postulated in FMS. 
We use here a modified version of these relations 
in order to take into account properly the symme- 
try of T J 

For the case of J=(4-), odd Z-A parity, and 
neglecting terms of order w/M, A/M, where w is the 
c.m. energy of the 7 meson and A is the Z2-A mass 
difference, the dispersion relations reduce to the 
approximate form® (T,,. _~!=T~) 
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dw’. 





“| Cw’) mT, “1(w!) 


lo 
m. ~_ 
i 


(2) 


B~'(w) is the inverse of the Born approximation 
matrix in terms of the renormalized coupling 
constants Enq 8dq° Faq SEK? and BEAK" All 
of the Born matrix elements have a single or 
double pole at w=0. C;(w) is a diagonal matrix 
defined in such a way that 


T(w) - B(w) as w-0. (3) 


The constant m i is the physical threshold value 
for w in each channel. We estimate that the er- 
rors introduced by our approximations are of the 
order of 10-30% in the energy region of interest. 

We now examine the matrix elements obtained 
from Eq. (2) to see if they give rise to resonance 
phenomena. In order to do this we make use of 
certain general properties® which are satisfied 
at resonance. The open-channel K matrix is de- 
fined as 


(T y ) -ik.6 (4) 


- Ng ig? 


if 


where of course these matrices (T", Kz) are 
the inverses of the open-channel parts of the T J 
and K 7 matrices, respectively. The condition 
for a resonance at w=w, is 


det K~*(w _)=0. (5) 
Yr 


If a resonance exists at w= Wy for the (Az, =7) sys- 
tem, then 


kK 


, A AA 
Ts" SM, kK . (6) 
fs ole FA, ~W)K 4 4]; m 
and 
r= . + T. (8) 
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where I, is the partial An width and T is the full 
width at half-maximum. Kaq and Kyy> are the 
elastic An and =7 K-matrix elements. 

The conditions (5) through (8), when applied to 
the elements obtained from Eq. (2), lead to im- 
plicit equations involving the four unknown coup- 
ling constants gay, 85q» SAK» 20d g>K (g). 
Another relation involving the coupling constants 
is obtained by setting 


T,,-=A 


KK (9) 


at the KN threshold, where A is the complex scat- 
tering length defined by Dalitz and Tuan.° We have 
searched for reasonable values of the coupling 
constants (g/g_,,<1.5) which would lead to the 
following experimental results!>?: 


w=2 wu, A=(-1.06+0.27) fermis, 


ay 


A 57 15, 


I'<50 Mev. 

The best known data are for w, and A. Fixing 
w, at 2u,’° we can use Eq. (5) to determine one 
of the unknown coupling constants in terms of the 
other three. For A within its quoted experimental 
error,' we find that the constant &sxK /4 n is almost 
uniquely determined to be 1.6 for a wide range of 
the others: 


<g 7/4n<5.1; 5.1<g 7/40 <32; 
0 ‘ /4n<5.1; 5.1 Bs /4m 
2 
1.3<g, /4n <32. 


The demand that we fit w, and A also leads to a 
narrow width, namely ['~11 Mev. If we allowA 
to have a somewhat larger real part, we obtain 
the following values for the coupling constants: 


2 ~" 2 = 
BoA /47=5.1, 85 /4n=32.0, 
2 a 2 an 
&aK /4n=11.5, ByK /4n=2.1, 


and the following experimental fits: 


A = -1.50 + 0.243, Pr ,/T,=3-4, T=38 Mev. 


It is gratifying that the scalar coupling constants 
&jq /40 and gy,"/47n are a good deal smaller than 
the pseudoscalar constant &N7'/4 nm. They are 
somewhat larger than those obtained by Ferrari 
and Fonda,” using criteria based on hypernuclear 
physics data, and by Nambu and Sakurai,’* who 
treat the 2 as a bound 7A system. The fact that 
the 2K coupling constant is almost uniquely de- 
termined seems to support Dalitz’s contention 


that the KN channel is most important in deter - 
mining the properties of this resonance. In fact, 
if we neglect the KN channels altogether, we find 
it impossible to fit the resonance data using our 
method. 

We have also briefly examined the J=(4-) case 
for even 2-A parity. In this case the Born approx- 
imation matrix elements have no poles near w=0. 
This is analogous to the S-wave 7-N scattering 
where we know that the Chew-Low approximation 
does not explain the results. 

From our results, we conclude that a resonance 
in the J=(4-) state can be explained in terms of 
odd =-A parity, whereas it would be difficult to 
understand in terms of even £-A parity. 

In a future publication we intend to present our 
calculations in greater detail, report on a more 
refined search program, and present other pre- 
dictions near thresholds. 

We would like to thank Mrs. Doris Ellis and 
Mr. M. Nussbaum for assisting with the compu- 
tations. 
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dation and the Office of Scientific Research. 
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NEW SOLID PHASE IN He*. J. H. Vignos and 
H. A. Fairbank [Phys. Rev. Letters 6, 265 
(1961) ]. 

The last factor in Eq. (1) should read (AVq,y/ 
AV,, 1) instead of (AV, 1/AV a, y)- 


LOSS OF C"' FROM PLASTIC FOILS AND ITS 
EFFECT ON CROSS-SECTION MEASUREMENTS. 
J. B. Cumming, A. M. Poskanzer, and J. Hudis 
[Phys. Rev. Letters 6, 484 (1961)]. 


The authors’ attention has been called to a paper 
on the loss of activity from betatron-irradiated 
samples by diffusion, by H. Fuchs and K. H. Lind- 
enberger, Nuclear Instruments and Methods 7, 
219 (1960), which previously reported the same 
effect for C'' produced by a (y,m) reaction. Their 
results and ours obtained for C™ produced by 
(p, pn) and (2, 2n) reactions are in good agree- 
ment which indicates that the effect is not very 
sensitive to the means of producing C™. 
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ABSTRACTS 


In this section are printed the abstracts of Articles 
that have been forwarded to The American Institute of 
Physics for publication in THE PHYSICAL REVIEW. 
In quoting information obtained from this section be- 
fore the appearance of the corresponding Article, ref- 
erence should be made to “Physical Review (to be 
published)” rather than to this Journal. 











ELECTRIC-FIELD-INDUCED ANISOTROPIES IN 
AN INHOMOGENEOUS PLASMA. Lewis Wetzel, 
Brown University, Providence, Rhode Island (Re- 
ceived February 21, 1961; revised manuscript 
received May 1, 1961). 


It is shown on the basis of the Boltzmann equa- 
tion that in the absence of a dc magnetic field, an 
electric field alone can produce anisotropies in 
electron diffusion and conductivity in a weak in- 
homogeneous plasma. These anisotropies have 
their origin in the ac component of the electron 
density resulting from the interaction between 
the electric field and electron density gradients. 
Electron diffusion is reduced in the direction of 
the field, while the induced ac current acquires 
an additional component in the direction of the 
gradient of the ac electron density. This compo- 
nent will in general lie in a direction different 
from that of the exciting field. Although these 
effects are usually small, they may become sig- 
nificant under suitable circumstances. 


ION MOTION IN SUPERFLUID LIQUID HELIUM 
UNDER PRESSURE. Lothar Meyer, Institute for 
the Study of Metals, University of Chicago, Chi- 
cago, Illinois, and F. Reif, Department of Phys- 
ics, University of California, Berkeley, Califor- 
nia (Received March 17, 1961). 


Recent investigations of superfluidity by a study 
of the mobilities of ions in liquid He II have been 
extended to the liquid under pressure. At a fixed 
temperature the positive-ion mobility decreases 
appreciably as the pressure is increased, partic- 
ularly at low temperatures. At a fixed pressure 
the mobility increases less rapidly with decreas- 
ing temperature at higher pressures. The nega- 
tive-ion mobility, smaller than that of the posi- 
tive ion at zero pressure, becomes equal to that 
of the latter above 7 atm. In high electric fields 
and at high pressures the drift velocity of the 


negative ions approaches a limiting value roughly 
equal to the Landau critical velocity for a body 
moving through the superfluid. The theory, which 
discusses the mobility in terms of ion scattering 
by rotons and phonons, is reviewed. It is pointed 
out that previously neglected effects concerned 
with the importance of small-angle scattering of 
the ion ought to be taken into account; some ear- 
lier estimates of scattering cross sections are 
revised accordingly. It is then shown that this 
theory, making use only of the known change of 
the roton dispersion relation with pressure, can 
account quantitatively for the observed pressure 
dependence of the positive-ion mobilities. 


RADIATION IN A PLASMA. I. CERENKOV EF- 
FECT. Marshall H. Cohen,* Observatoire de 
Meudon, Seine et Oise, France (Received March 
27, 1961). 


We start with the linearized plasma equations 
containing an isotropic pressure term, plus ex- 
tra source terms J° and p* in the Maxwell equa- 
tions. The fields of (J*, p*) can be decomposed 
into two modes. The electromagnetic (EM) mode 
has all the magnetic field and no charge accumu- 
lation; it is the ordinary EM field of (J°, p*) in 
a dispersive medium of relative dielectric con- 
stant €,=1-(w»/w)’. The plasma (P) mode has 
all the charge accumulation and no magnetic 
field; at great distances from the source it be- 
comes a longitudinal (radial) plasma wave with 
the usual dispersion relation for plane plasma 
waves. Various potentials for the EM and P 
modes are given by the inhomogeneous Klein- 
Gordon equation. 

The fields of a uniformly moving charged par- 
ticle are found by a Lorentz transformation. 
When (u/v,)< 1 (u=particle velocity, v,=rms 
thermal velocity) the EM and P fields are ex- 
ponentially screened outside oblate spheroids 
foreshortened in the direction of motion. When 
(u/v,)> 1 the field exists only within the Mach 
(Cerenkov) cone trailing the particle. The fre- 
quency and angular spectra of the Cerenkov radi- 
ation are found, and the total radiated energy is 
found by assuming an arbitrary high-frequency 
cutoff due to Landau damping. The expression 
for total radiated energy agrees with that given 
by Pines and Bohm, except for the logarithmic 
terms. 

*On leave from Cornell University, Ithaca, New York. 
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SELF-CONSISTENT PAIR INTERACTION FOR 
MANY- FERMION SYSTEM. Katuro Sawada* and 
Toshio Soda, University of California, La Jolla, 
California (Received August 15, 1960). 


The equations which determine the one-particle 
energy and effective two-body interaction in an 
interacting Fermi gas are constructed within the 
approximation which sums up all pair creation- 
annihilation processes. The equation corresponds 
to the familiar equation for the K matrix which 
represents the interaction between particles (or 
holes) and sums up the particle-particle (or hole- 
hole) scattering processes. The method of the 
equation of motion is used in this paper. Our re- 
sult for the one-particle energy is shown to lead 
to the result previously obtained by Quinn and 
Ferrell, and by Rockmore for the case of the elec- 
tron gas with Coulomb interactions, when we re- 
place screened potentials by bare potential in the 
self-consistent energy equation. For nuclear 
matter, it is shown that the presence of an attrac- 
tive interaction in the equation of motion for num- 
ber density causes an “enhancement” of exchange 
forces, whereas in the electron gas repulsive 
Coulomb interactions lead to “screening” of the 
exchange force. The strength of the isospin den- 
sity interaction pseudopotential is enhanced by 
a factor of two when one solves the self-consis- 
tent equation; and a simple estimate shows that 
the Goldhaber-Teller mode lies about 15% higher 
than the value bra/m, previously estimated by 
Glassgold et al. (g: momentum of the oscillation, 
pr: Fermi momentum). 


*On leave of absence from Tokyo University of Educa- 
tion, Tokyo, Japan. 


LATTICE ANHARMONICITY AND OPTICAL AB- 
SORPTION IN POLAR CRYSTALS. Il. CLASSICAL 
TREATMENT IN THE LINEAR APPROXIMATION. 
A. A. Maradudin, Westinghouse Research Labora- 
tories, Pittsburgh, Pennsylvania, and R. F. Wallis, 
U. S. Naval Research Laboratory, Washington, 

D. C. (Received March 24, 1961). 


An expression for the elements of the dielectric 
susceptibility tensor for an ionic crystal has been 
derived in a manner analogous to that employed 
by Kubo in his treatment of magnetic susceptibil- 
ity. In the high-temperature (classical) limit, 
this expression reduces to the Laplace transform 
of the autocorrelation function of the single nor- 
mal coordinate which interacts directly with the 
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external radiation field. We have applied this 
formalism to the calculation of the high-temper- 
ature linear optical absorption coefficient of an 
anharmonic ionic crystal for which we have re- 
tained only cubic anharmonic terms in the lattice 
Hamiltonian. We have solved directly for the 
Laplace transform of the autocorrelation function 
to lowest order in the anharmonic coupling con- 
stant by solving the equations of motion for the 
normal coordinates after they have been linear- 
ized. The linear absorption coefficient obtained 
from the susceptibility tensor is of Lorentzian 
form with a frequency-dependent damping constant 
which varies linearly with temperature. The ab- 
sorption coefficient for a diatomic linear chain 
with nearest neighbor interactions has been eval- 
uated. 


RECOIL-FREE RESONANT ABSORPTION IN 
Au’, D. A. Shirley and M. Kaplan, Lawrence 
Radiation Laboratory and Department of Chem- 
istry, University of California, Berkeley, Cali- 
fornia, and P. Axel, Physics Department, Univer- 
sity of Illinois, Urbana, Illinois (Received March 
27, 1961). 


The Mdssbauer absorption in Au has been meas- 
ured at 4°K for the 77-kev gamma ray emitted by 
Au’™ nuclei embedded in gold, platinum, stain- 
less steel, iron, cobalt, and nickel. In each case, 
a Doppler -shift curve was measured to find the 
effective width and the chemical shift. The recoil- 
free fractions, f, are obtained with the aid of a 
straightforward analysis, which incidentally shows 
the errors that can be made in f if the chemical 
shift and effective width are not taken into account. 
The observed recoil-free emission fractions were 
found to be approximately 0.06(Au), 0.34 and 
0.14(Pt), 0.24(steel), 0.32(Fe), 0.27(Co), and 
0.35(Ni). The relative f values are correct; the 
correct absolute f values might require a multi- 
plicative correction factor that could be as small 
as 0.53. 

Relatively large f values were obtained when 
the Au’™ radioactive nuclei were in high Debye- 
temperature lattices composed of light nuclei. 
Particularly low f values were found when the 
largest radiation-damage effects were expected. 

The observed chemical shifts were (in units of 
107° ev) <0.13(Au), 0.26(Pt), 1.3(steel), 1.4(Fe), 
1.3(Co), and 1.1(Ni). These chemical shifts give 
information more directly interpretable than, but 
related to, the optical-isotope and isomer shifts, 
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and the Knight shift; the magnitudes of the shifts 
measure directly either the depletion or the en- 
hancement of electron density at the radioactive 
nuclei in the different lattices. 

Zeeman splittings of the nuclear energy levels 
caused by local magnetic fields at Au had relative 
magnitudes of 1.0:0.43:<0.10 for Fe, Co, and Ni 
lattices. If the magnetic moment of the 77-kev 
excited state in Au’®”’ is 1.6 nm, the local mag- 
netic fields in Fe and Co were 282 and 122 koe, 
respectively. 


LATTICE DISTORTION DUE TO ISOTOPES IN 
SOLID HELIUM. P. G. Klemens and A. A. 
Maradudin, Westinghouse Research Laboratories, 
Pittsburgh, Pennsylvania (Received March 17, 
1961). 


Observations of the thermal resistivity of solid 
helium containing He* isotopes have indicated 
that the phonon scattering cross section of the 
isotope is three times larger than the theoretical 
value obtained by one considering only the mass 
difference. Using as a model a helium atom in 
the potential well formed by its neighbors, and 
taking account of the zero-point energy of that 
atom and of the strain energy in the surrounding 
solid when the atomic cell is expanded, it is 
shown that around a He® isotope the lattice suf- 
fers an outward distortion of about 1-2%. This 
distortion can account for the discrepancy in the 
phonon scattering cross section. 


MAGNETIC PROPERTIES OF HCrO, AND 
DCrO,. R. G. Meisenheimer and J. D. Swalen, 
Shell Development Company, Emeryville, Cali- 
fornia (Received March 13, 1961). 


Magnetic susceptibility measurements and elec- 
tron paramagnetic resonance measurements have 
been made on powdered chromous acid (HCrO,) 
and deuterated chromous acid (DCrO,). The 
Curie constants, derived from the susceptibility 
results, are consistent with three unpaired elec- 
trons associated with each chromium ion. In 
addition, a large Curie-Weiss 6 temperature, 
indicating a strong exchange interaction, is neces- 
Sary to account for the results. The most unusual 
feature of this particular exchange is that it is 
strongly influenced by deuteration. 

The electron paramagnetic resonance spectra 
also show strongly exchange-narrowed lines 





with a variation between HCrO, and DCrO,. From 
a comparison between the calculated and observed 
linewidths and second moments, a value of the 
zero-field splitting is estimated. 

Although no transition has been observed to an 
antiferromagnetic state, the Curie-Weiss 6 tem- 
peratures indicate that the ground state is un- 
doubtedly antiferromagnetic for both HCrO, and 
DCrO,. Superexchange through the intervening 
proton or deuteron is probably the main source 
of any antiferromagnetic exchange coupling. 

Some of the various mechanisms proposed for 
superexchange are discussed in relation to chro- 
mous acid. 


ABSORPTION SPECTRA OF TRANSITION IONS 
IN CdS CRYSTALS. R. Pappalardo and R. E. 
Dietz, Bell Telephone Laboratories, Murray 
Hill, New Jersey (Received March 21, 1961). 


The low-temperature optical absorption spec- 
tra of single crystals of wurtzite CdS containing 
impurity ions of the first transition series and 
one rare earth ion, ytterbium, were studied in 
the range of 4.5 to 0.4 u. The main features of 
the spectra were found to be in good agreement 
with a cubic crystal field model, while structure 
found in the spectral bands of nickel and cobalt 
impurities could be described by a simple, first- 
order treatment using the free ion spin-orbit 
coupling constants. The relevant crystal field 
parameters, the site symmetry, and the formal 
charges of the impurity ions were determined 
where possible. 


RECOVERY OF ELECTRON RADIATION DAM- 
AGE IN 2-TYPE InSb. F. H. Eisen, Atomics 
International, a division of North American 
Aviation, Incorporated, Canoga Park, California 
(Received March 22, 1961). 


The production and recovery of electron radia- 
tion damage in n-type InSb has been studied by 
means of Hall coefficient and electrical conduc- 
tivity measurements. Irradiations were per- 
formed mainly at 80°K, since no recovery was 
observed between 4°K and 80°K. The damage re- 
covered in five well-defined stages with the re- 
covery nearly complete at 320°K. Isochronal 
and isothermal recovery was monitored in each 
of the stages, allowing a determination of the 
activation energies for recovery and a study of 
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the recovery kinetics. None of the recovery 
kinetics fit any simple models. There is evidence 
that the two lowest temperature recovery stages 
involve the annihilation of close interstitial- 
vacancy pairs and that interactions of primary 
defects with impurities do not occur. However, 
the first-order kinetics expected for close-pair 
recovery is not explicitly observed. A possible 
explanation for the observed kinetics, involving 
the independent annihilation of two types of close- 
pair configurations in the same stage with an 
electrostatic interaction between the interstitial 
and vacancy, is proposed. 


LATTICE ABSORPTION BANDS IN SiC. Lyle 
Patrick and W. J. Choyke, Westinghouse Research 
Laboratories, Pittsburgh, Pennsylvania (Received 
March 29, 1961). 


Using four phonon energies, TA =0.045, LA 
=0.067, TO=0.0955, and LO =0.1055 ev, it is 
possible to explain, as summation bands, the 
ten absorption bands lying in the energy range 
0.130 to 0.300 ev. Phonon energies close to three 
of the four given here have been obtained by an 
analysis of the indirect interband absorption in 
SiC. 


REFLECTION OF NOBLE GAS IONS AT SOLID 
SURFACES. Homer D. Hagstrum, Bell Telephone 
Laboratories, Murray Hill, New Jersey (Re- 
ceived March 22, 1961). 


A method has been developed for distinguishing 
the reflection of ions at solid surfaces as ions 
or as metastable atoms. Results are given for 
Het, Ne*, and Ar* ions incident on clean W, Mo, 
and Si(100) and on contaminated W, Hf, and 
Ge(111) surfaces. The reflection coefficient of 
ions to ions (Rjj) is found to be small (0.0004 to 
0.002) and essentially independent of incident ion 
energy. The reflection coefficient of ions to 
metastable atoms (R;,,) is found to increase with 
ion energy from values comparable to Rj; at 10 
ev to values as high as 0.04 at 1000 ev. Dis- 
cussion of these results is given in terms of the 
known resonance and Auger transitions which 
can occur near a solid surface for ions of suffi- 
ciently large ionization energy. It is shown that 
the results can be accounted for only if ions are 
transformed to metastable atoms very close to 
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the surface and a possible mechanism for this 
process is proposed. 


EFFECT OF HYDROSTATIC PRESSURE ON IONIC 
CONDUCTIVITY IN DOPED SINGLE CRYSTALS 
OF SODIUM CHLORIDE, POTASSIUM CHLORIDE, 
AND RUBIDIUM CHLORIDE. C. Ballard Pierce,* 
Department of Physics, University of Illinois, 
Urbana, Illinois (Received January 30, 1961; re- 
vised manuscript recelved May 8, 1961). 


The effect of hydrostatic pressure up to 9000 
kg/cm? on ionic conductivity in NaCl, KCl, and 
RbCl single crystals doped with divalent impuri- 
ties has been studied over the temperature range 
200°C to 500°C. The conductivity in this temper- 
ature range is due almost entirely to the motion 
of extrinsic cation vacancies. The activation vol- 
ume, AV,,, for motion of the cation vacancies is 
7.7+0.5 cc/mole in NaCl doped with CaCl, and 
7.0+ 0.5 cc/mole in KCl doped with SrCl,. The 
results are in fair agreement with values pre- 
dicted on the basis of Keyes’s empirical expres- 
sion relating activation volume to activation ener- 
gy and isothermal compressibility. Sample ma- 
terials were chosen with the view of testing for a 
correlation between activation energy and shear 
modulus. The small shear modulus in KCl and 
RbCl decreases with increasing pressure, while 
the reverse is true for NaCl. However, the data 
are not adequate to draw definite conclusions a- 
bout such a correlation. The conductivity of RbCl 
doped with BaCL, increases by an order of mag- 
nitude at the phase transition from the NaCl to 
the CsCl structure. At 300°C, the transition oc- 
curs at 6100 kg/cm’. 


*Now at Sandia Corporation, Albuquerque, New Mexi- 
co. 


PHOTOSENSITIVE-ULTRASONIC PROPERTIES 
OF CADMIUM SULFIDE. Harmon D. Nine, Gen- 
eral Motors Research Laboratories, Warren, 
Michigan, and Rohn Truell, Brown University, 
Providence, Rhode Island (Received January 9, 
1961). 


Ultrasonic attenuation in single crystals of CdS 
has been observed to be a function of light irradi- 
ation. Two distinct types of behavior have been 
observed. Some crystals (type A) show a de- 
crease of attenuation with white light application, 
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and others (type B) show an increase of atten- 
uation with white light. The ultrasonic attenua- 
tion and the conductance of the crystals were 
measured as a function of temperature. An ex- 
cellent correlation between the conductance and 
the ultrasonic attenuation is evident for type B 
cyrstals. A strong correlation between the di- 
rectional piezoelectric properties of CdS and 

the ultrasonic attenuation in both types of CdS is 
also shown. An interaction between conduction 
electrons and stress waves coupled by the piezo- 
electric property of CdS is isolated as the domi- 
nant mechanism for the photosensitive ultrasonic 
attenuation in CdS of type B. Two or more com- 
peting effects appear to be responsible for the 
type A photosensitive attenuation behavior. 


LATTICE DYNAMICS OF ALPHA URANIUM. 

D. O. Van Ostenburg, Argonne National Labora- 
tory, Argonne, Illinois (Received September 14, 
1960). 


The method developed by Begbie and Born has 
been applied to alpha uranium, where equations 
are developed which give the macroscopic elas- 
tic constants in terms of the microscopic force 
constants. Interactions of an atom with its first 
through fourth nearest neighbors are considered; 
these involve twelve atoms. Through symmetry 
considerations, nineteen atomic force constants 
enter into this force system. An independent 
determination of the force constants is required 
before a valid verification of the solutions can 
be made. However, using measured values of 
the nine elastic constants, two sets of force 
constants are evaluated, one based upon central 
forces and the other upon neglect of fourth 
nearest neighbors. 


Fe’ MOSSBAUER EFFECT IN Cu-Ni ALLOYS. 
G. K. Wertheim and J. H. Wernick, Bell Tele- 
phone Laboratories, Murray Hill, New Jersey 
(Received March 30, 1961). 


The Moéssbauer effect of Fe*’ has been used to 
Study the properties of iron impurity atoms in 
the complete range of composition of Cu-Ni al- 
loys. The isomer shift indicates a small de- 
crease in the total electronic density at the Fe®*”’ 
nucleus in going from pure Ni to pure Cu. The 
magnetic field at the iron nucleus at 0°K de- 


creases by 9% in the range from 0 to 40% Cu. 
Both results show that only minor changes take 
place in the atomic configuration of the iron. 
The linewidth, measured in the paramagnetic 
alloys, is smallest in pure copper and largest 
near the middle of the composition range. The 
rapid change in linewidth with small nickel ad- 
mixture can result from quadrupole splitting 
due to the field gradients arising from the spatial 
charge fluctuations around an impurity atom. 
Some of the broadening can also arise from an 
inhomogeneous isomer shift due to a range of 
surroundings of an iron atom in the alloy. 


EFFECTIVE X-RAY AND CALORIMETRIC 
DEBYE TEMPERATURE FOR COPPER. P. A. 
Flinn, G. M. McManus, and J. A. Rayne, West- 
inghouse Research Laboratories, Pittsburgh, 
Pennsylvania (Received March 27, 1961). 


The Debye-Waller factor for copper was deter- 
mined from x-ray intensity measurements on a 
single crystal over the temperature range 4.2°K- 
500°K. From a machine calculation of the vibra- 
tional spectrum of copper, the values of the 
specific heat and Debye-Waller factor were ob- 
tained and compared with those found by experi- 
ment. The agreement indicates that for copper 
the central force model with nearest and second 
neighbor interactions is adequate for the inter- 
pretation of effects depending on simple aver- 
ages over the frequency spectrum. 


IMAGE OF THE FERMI SURFACE IN SPIN-WAVE 
SPECTRA OF RARE EARTH METALS. Edwin J. 
Woll, Jr., and Stephen J. Nettel, Department of 
Physics, University of California, San Diego, 
California (Received March 20, 1961). 


Calculations of spin-wave spectra in rare earth 
metals were carried out to find whether images 
of the electronic Fermi surface might be observ- 
able. In the space of spin-wave vectors q there 
should occur surfaces on which the frequencies 
have an infinite gradient with respect tog, the 
location of such abrupt changes, “kinks” in the 
dispersion curves, being determined by the shape 
of the Fermi surface. The spin-wave spectrum 
is found by assuming that the coupling between 
ionic spins takes place primarily through exchange 
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scattering of conduction electrons, paralleling 
the calculation on the coupling of nuclear spins 
by Ruderman and Kittel. Spin-wave dispersion 
curves in two directions of high symmetry are 
computed numerically. The sought-for kinks in 
the dispersion curves are found to amount to 
about two percent of the maximum excitation 
frequency. The development is for ferromag- 
nets, but extension to spiral antiferromagnets 
is taken up briefly. 


ANGULAR DISTRIBUTIONS OF SPUTTERED PO- 
TASSIUM ATOMS. R. P. Stein* and F. C. Hurl- 
but, Missiles and Space Division, Lockheed Air- 
craft Corporation, Sunnyvale, California (Re- 
ceived January 19, 1961; revised manuscript 
received April 24, 1961). 


The angular distributions of potassium parti- 
cles issuing from a potassium surface under 
bombardment by noble gas ions were observed 
under moderately good vacuum conditions. Sput- 
tered potassium atoms were detected for inci- 
dent ion energies above approximately 15 elec- 
tron volts and useful observations of angular dis- 
tributions were obtained for incident ion energies 
in the range 50 to 450 electron volts for all avail- 
able values of the incident angle. A means was 
discovered for the discrimination between the 
total sputtered flux and that fraction of it pos- 
sessing particle energies above a certain thresh- 
old. The apparatus and experimental procedures 
are described and the observed distribution pat- 
terns and a two-collision sputtering mechanism 
are discussed, along with related observations. 


*R. P. Stein, Lockheed Missiles and Space Division, 
Palo Alto, California. 

TF. C. Hurlbut, University of California, Berkeley, 
California. 


FLUORESCENCE AND OPTICAL MASER EF- 
FECTS IN CaF,:Sm**. C. G. B. Garrett, W. 
Kaiser, and D. L. Wood, Bell Telephone Labo- 
ratories, Murray Hill, New Jersey (Received 
March 17, 1961). 


Measurements are reported of absorption, 
emission, and activation spectra in CaF,:Sm*t, 
and also of fluorescence lifetime. A revised 
level scheme is proposed. Observations of opti- 
cal maser effects were made at liquid hydrogen 
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and liquid helium temperatures over a wide range 
of pumping intensities. For illumination in the 
6400A band, the threshold intensity of illumina- 
tion was about 20 watts cm~*. Observations are 
reported of the dependence of the intensity of the 
maser beam on the pumping intensity. Photo- 
graphic observations, made both with and without 
a Fabry-Perot etalon, show that on one particular 
sample five distinct frequencies were present in 
the maser signal, and that the number of modes 
excited was of the order of one thousand. 


MAGNETIC MOMENT OF LANTHANUM MAG- 
NETOPLUMBITE FERRITE. A. Aharoni and 
M. Schieber, Department of Electronics, The 
Weizmann Institute of Science, Rehovoth, Israel 
(Received March 31, 1961). 


The composition La*+Fe’tFe,,**+O,, was pre- 
pared by sintering the oxides, and its magnetic 
moment was measured from liquid air tempera- 
ture to the Curie point. At the lowest tempera- 
ture a value of 17.5 Bohr magnetons per formula 
weight was obtained, as compared with 19 Bohr 
magnetons of BaFe,,°*O,,. The difference is 
qualitatively understood if the Fe** ion is a near- 
est neighbor to the La‘* in the crystal, i.e., in 
the [12 K] site, thus decreasing the total magnetic 
moment. The Curie temperature is 695°K, com- 
pared to 725°K of BaFe,,**O,,. The lower Curie 
point can again be qualitatively explained by less 
exchange interaction with the divalent iron of 
LaFe’*Fe,,°*O,,, or possibly by the change in 
unit-cell dimension which modifies the short- 
range exchange interaction between the iron ions. 


REMANENT STATE IN ONE-DIMENSIONAL MI- 
CROMAGNETICS. Amikam Aharoni, Department 
of Electronics, The Weizmann Institute of Sci- 
ence, Rehovoth, Israel (Received March 1, 1961; 
revised manuscript received April 20, 1961). 


A first integral is found for Brown’s nonlinear 
equations in one dimension. When the external 
field is zero, another first integral can be found, 
which enables complete integration of the equa- 
tions. For a unidirectional anisotropy with an 
easy direction perpendicular to the plane of the 
film, one of the integration constants is not de- 
termined by the boundary conditions, as if in- 
dicating a possibility of continuum of different 
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remanence values for different histories. How- 
ever, when a small field in the plane of the film 
is introduced as a perturbation, this degeneracy 
is removed, and the magnetization can change 
only in the plane defined by the field and the di- 
rection of anisotropy. There are still many 
discrete possible values for the remanence, 
each of which determines uniquely the suscepti- 
bility at that state. 


DISPERSION OF GYROMAGNETIC RATIOS IN 
COMPLEX SPECTRA. Norbert Rosenzweig, 
Argonne National Laboratory, Argonne, Illinois, 
and C. E. Porter, Brookhaven National Labora- 
tory, Upton, New York (Received March 20, 
1961). 


It is pointed out that the random-matrix hypoth- 
esis, which was previously used to explain the 
“repulsion” of energy levels in complex atomic 
and nuclear spectra, leads to a dispersion in the 
values of the gyromagnetic ratios which depends 
on the relative strengths of central and spin- 
dependent interactions and the number of states 
which interact with one another in the LS coupling 
scheme. The dispersion based on experimentally 
known atomic g factors is computed in three re- 
gions of the periodic table, and the empirical 
trends are found to be in qualitative agreement 
with the theoretical analysis. 


DECAY OF Tm’”. R. G. Helmer and S. B. Bur- 
son, Argonne National Laboratory, Argonne, 
Illinois (Received March 16, 1961). 


The radioactive nuclide ,,Tm’” was produced 
by successive capture of two neutrons in erbium 
oxide enriched in Er’. The irradiations were 
made in the Materials Testing Reactor at Arco, 
Idaho. In addition to three thulium activities, 
these samples contained six active contaminants. 
Pure thulium sources were obtained by use of an 
ion-exchange column. Studies were conducted 
with a 256-channel coincidence scintillation spec- 
trometer. These measurements indicate the pres- 
ence of at least 17 gamma-ray and 5 beta-ray 
transitions. The beta-ray spectrum was studied 
with a 180° magnetic beta-ray spectrometer. 

This spectrum was analyzed by use of a com- 
puter program compiled by the authors in col- 
laboration with members of the Argonne Applied 


Mathematics Division. The level scheme pro- 
posed for Yb’” has states with energies, spins, 
and parities of 0.0(0*), 0.079(2*), 0.260(4*), 
1.17(3), 1.46(2), 1.54(3), 1.60(1), 1.64(?7), and 
1.73(3) Mev. The total decay energy is found to 
be 1.88 Mev. The experimental data are con- 
sistent with the previously proposed interpre- 
tation that the first two excited states are mem- 
bers of a K =0 rotational band based on the ground 
state. The states at 1.46 and 1.54 Mev are inter- 
preted as members of a rotational band with 
K=2. The states at 1.60 and 1.73 Mev are tenta- 
tively interpreted as members of a rotational 
band with K =0 and negative parity. It is sug- 
gested that the state at 1.17 Mev has K=3. From 
the analysis of the beta spectrum it is concluded 
that the ground state of thulium has J/=K =2 and 
negative parity. 


ENERGY LEVELS IN atu™ FROM THE DE- 
CAY OF ,,Er’”. R. G. Helmer and S. B. Burson, 
Argonne National Laboratory, Argonne, Illinois 
(Received March 27, 1961). 


The radioactive nuclide ,,.Er’” was produced 
in the Materials Testing Reactor, Arco, Idaho, 
by the successive capture of two neutrons in 
erbium oxide enriched in Er’”°. In addition to 
three erbium activities and the radioactive Tm’” 
daughter, these samples contained six active 
contaminants from which the erbium was sepa- 
rated by use of an ion-exchange column. Scin- 
tillation studies, conducted with a 256-channel 
coincidence scintillation spectrometer, indicate 
the presence of at least eight gamma-ray transi- 
tions. Two of these transitions are either highly 
K converted or their transition energies are 
approximately 50 kev. The energies of the other 
six transitions are about 610, 450, 408, 200, 
160, and 125 kev. Beta-gamma coincidence ex- 
periments indicate the presence of two beta-ray 
components at approximately 310 kev (in coinci- 
dence with the 610-kev gamma ray) and 370 kev 
(in coincidence with the 408-kev gamma ray). 
The level scheme deduced for Tm'” has excited 
states at 408, 450 (or 160), 470+ 15, 530, and 
610 kev. From the beta-ray intensities, the 
states at 530 and 610 kev are assigned spins of 
either 0 or 1 with positive parity. The spin and 
parity (2~) of the thulium ground state have pre- 
viously been assigned from the properties of the 
decay of Tm?”. 
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INTERPRETATION OF EXPERIMENTAL (n, 2n) 
EXCITATION FUNCTIONS. D. W. Barr, C. I. 

Browne, and J. S. Gilmore, Los Alamos Scien- 
tific Laboratory, University of California, Los 
Alamos, New Mexico (Received March 6, 1961). 


Radiochemically determined (m, 2n) excitation 
functions for Sc**, Ti*®, Ni®, Cu®®, Ge”, As”, 
sr*™*, Rb*®, Rb*”, y*®, Zr”, Sn! , Cal, Sb'?3, 
and U** have been interpreted in terms of the 
statistical model of nuclear reactions. Values 
of the level density parameter a are obtained and 
correlated with mass number. A procedure is 
outlined for predicting the magnitude of any 
(n, 2n) cross section from the nuclear content of 
the target material. Two level density formula- 
tions are studied, and approximations customar- 
ily made in calculations of this sort are examined 
quantitatively. 


MEASUREMENT OF THE DEUTERON BINDING 
ENERGY USING A BENT-CRYSTAL SPECTRO- 
GRAPH. A. Halim Kazi,* Norman C. Rasmussen, 
and Hans Mark, Department of Nuclear Engineer - 
ing and Department of Physics, Massachusetts 
Institute of Technology, Cambridge, Massachusetts 
(Received March 20, 1961). 


The deuteron binding energy has been determined 
by measuring the neutron-proton capture gamma- 
ray energy. This energy has been measured dir- 
ectly, relative to annihilation radiation, with the 
help of a six-meter radius bent-crystal spectro- 
graph. The spectrograph is of the Cauchois type, 
in which a collimated but extended gamma-ray 
beam is incident on the convex side of an elasti- 
cally bent quartz crystal, is diffracted, and focused 
onto a focal circle defined by the radius of curva- 
ture of the crystal. The neutron-proton capture 
gamma rays are produced by placing a polyethyl- 
ene sample in the through port of the MIT research 
reactor. The (310) planes of quartz are used for 
diffraction, and the gamma-ray lines are recorded 
on glass-mounted 600-micron thick Ilford G-5 
emulsions. The value of B(D) obtained is 2225.5 
+1.5 kev, where the error is the standard devia- 
tion. This is the most precise direct measure- 
ment reported to date, and is in agreement with 
previous work. Using recent mass spectroscopic 
data, the mass of the neutron is found to be 
1.008 984 + 0.000002 amu. 

The efficiency of the spectrograph is low. At 
2225 kev, 6000 curie hours are required to re- 
cord a line; at 511 kev, 1200 curie hours are 
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necessary. The error in B(D) agrees with the 
estimated precision which varies from about 
0.01 % at 100 kev to 0.3% at 4000 kev. The latter 
energy is close to the practical upper energy 
limit of the instrument. 


*Present address: John Jay Hopkins Laboratory for 
Pure and Applied Science, General Atomic Division of 
General Dynamics Corporation, San Diego, California, 


NUCLEAR ORIENTATION OF Pm". C. A. Love- 
joy, J. O. Rasmussen, and D. A. Shirley, Law- 
rence Radiation Laboratory and Department of 
Chemistry, University of California, Berkeley, 
California (Received March 24, 1961). 


Promethium- 143 has been oriented in a crystal 
of neodymium ethyl-sulfate. The angular distri- 
bution of the 740-kev » ray was found to be W(@) 
= 1-(0.060+ 0.006)P,(cos9) at 0.02°K. Values for 
the mixing ratio, 5, of the 740-kev y» ray of Nd 
were obtained as a function of the magnetic mo- 
ment of the ground state of Pm™**. The spin of 
the excited state of Nd‘** was assigned as 9/2-. 
An absolute lower limit of || > 1.0 was set on 
the magnetic moment of Pm***. The mixing ratio 
of the y ray of Nd'** was found to lie in the range 
0.23 < 5(E2/M1)< 0.35. 


Ca**( p, y)Sc** REACTION. J. W. Butler, Nucleon- 
ics Division, U. S. Naval Research Laboratory, 
Washington, D. C. (Received March 24, 1961). 


Proton capture by Ca*® has been studied by 
means of observations on both the prompt gamma 
rays from resonance states in the compound nu- 
cleus and the delayed positrons from the decay of 
the ground state. Targets of CaO were prepared 
by the electrodeposition of Ca onto a Pt backing 
followed by the oxidation of the Ca and the purging 
of impurities by heat. These targets were bom- 
barded by protons from a 2-Mv Van de Graaff 
accelerator, producing the reaction Ca*(p, )Sc*. 
The gamma rays from this reaction were observed 
with the use of a 3-in. diam by 3-in. Nal(T1) crys- 
tal and a 256-channel pulse-height analyzer. Pos- 
itrons from the decay of Sc** were detected with 
the use of a thin plastic phosphor, 1.5 in. diam by 
0.012 in. thick. Two resonances in the reaction 
were observed at bombarding energies of 650+ 5 
kev and 1850+ 10 kev. Two other possible reso- 
nances were observed at 1550+ 15 kev and 
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1630+ 15 kev. The 650-kev resonance corre- 
sponds to an excited state in Sc** at 1.7234 0.011 
Mev, has an integrated cross section of 0.02 ev 
barn (factor of two uncertainty either way), has 
a width of less than 5 kev, and involves a gamma 
ray whose energy was measured to be 1.71+ 0.03 
Mev. The 1850-kev resonance corresponds to an 
excited state in Sc*! at 2.883+ 0.014 Mev, has an 
integrated cross section of 0.3 ev barn (factor of 
two uncertainty either way), has a width of less 
than 10 kev, and involves a gamma ray whose 
energy was measured to be 2.894 0.02 Mev. 


NUCLIDES Ar* and Cl°*. Nelson Jarmie and 

M. G. Silbert, Los Alamos Scientific Laboratory, 
University of California, Los Alamos, New Mex- 
ico (Received March 27, 1961). 


The masses of Cl** and Ar* and the energies 
of their first excited states have been determined 
by an investigation of the reactions Ar*(t, a)Cl°® 
and Ar*%(¢,p)Ar*?. Charged reaction products were 
analyzed with a high-resolution magnetic spectrom- 
eter. The experimental results for Ar® are: 
mass excess (M-A) =-34.423+ 0.040 Mev (C’* =0) 
or -21.990+ 0.040 Mev (O'*=0); energy of first 
excited state 1.138+0.030 Mev. The experimental 
results for Cl°* are: mass excess (M-A)=-29.772 
+ 0.040 Mev (C’*=0) or -18.227+ 0.040 Mev 
(O**=0); energy of first excited state = 0.364 
+0.030 Mev. The Q values of Ar*°(¢, p)Ar*? and 
Ar*°(t, a)C1°® were found to be 7.0464 0.040 Mev 
and 7.259+ 0.040 Mev, respectively. 


ELASTIC SCATTERING OF IDENTICAL SPIN- 
ZERO NUCLEI. D. A. Bromley,* J. A. Kuehner, 
and E. Almqvist, Atomic Energy of Canada Lim- 
ited, Chalk River Laboratories, Chalk River, 
Ontario, Canada (Received March 27, 1961). 


Elastic scattering measurements have been 
carried out on the C’*+C™ and O'*+0"* systems 
in the energy range from 6 to 35 Mev using heavy- 
ion beams from the Chalk River tandem accelera- 
tor and Au-Si surface barrier detectors. At en- 
ergies below the Coulomb barriers the Mott scat- 
tering predictions are in excellent accord with 
the measurements as functions of both angle and 
energy. At energies above the barrier the 0+O 
excitation curve drops exponentially below the 
Mott predictions to a value of 10 mb/sr at E, m. 
= 16.5 Mev and remains approximately constant 


thereafter at that value; in contrast, the C+C 
excitation curve shows marked resonant inter- 
ference structure. The states involved have 
T~10~" sec, large compound elastic branching 
ratios T c/T, and appear to correspond to resonant 
absorption of high order partial waves; hence they 
have high angular momentum. It is suggested that 
these states reflect a quasi-molecular interaction 
mechanism which is involved in grazing collisions 
of these nuclei and which is critically dependent 
upon the structure of the nuclei involved. 


*Present address: Department of Physics, Yale 
University, New Haven, Connecticut. 


“BREAKS” IN THE ACTIVATION CURVES OF 
THE P**(y,n)P®° REACTION. D. Sadeh,* Labora- 
toire de Synthése Atomique, Centre National de 

la Recherche Scientifique, Ivry, France (Received 
March 17, 1961). 


Discontinuities in the slope of the activation 
curve of the P*"(y,n)P®° reaction were found while 
using two large Nal(T1) crystals to detect the an- 
nihilation gammas from P**. These breaks cor- 
respond well with resonances found recently in 
the Si°°(p,)P°° reaction. 

The correspondence between breaks and known 
resonances in some light nuclei is discussed. 
Such a correspondence was found to exist in the 
case of N“, F®, and P* and is in doubt in the 
case of C™, 

The detection system used allowed accurate 
measurements of the thresholds of the P*"(y,n)P*° 
reaction -(12.23+ 0.04 Mev) and the Cl**(y,”)C1™* 
reaction (12.66+0.04 Mev). 


*Present address: Atomic Energy Commission of 
Israel, P.O.B. 7056, Hakirya, Tel-Aviv, Israel. 


T=3/2 LEVEL IN F’® AT 9.07 Mev. G. Amsel 
and G. R. Bishop, Ecole Normale Supérieure, 
Paris, France (Received March 24, 1961). 


Contrary to previous reports we have detected 
a@ particles emitted from the 9.07-Mev state of 
F’* formed by the reaction O**(p, a)N’®. The a- 
particle reduced width indicates that the isotopic 
spin of this state is T=3/2 with an admixture of 
about 5% of T=1/2. The spin of the state is de- 
termined to be 7/2 from the a-particle angular 
distribution, while odd parity is favored by the 
proton reduced width. 
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NEUTRON CAPTURE CROSS SECTION OF GOLD 
AT 30 kev AND 64 kev. L. W. Weston and W. S. 
Lyon, Oak Ridge National Laboratory, Oak Ridge, 
Tennessee (Received January 27, 1961; revised 
manuscript received May 8, 1961). 


The capture cross section of gold has been 
measured with kinematically collimated neutrons 
from the Li’(p,)Be’ and T(p,n)He® reactions. 
The cross sections at 30.2 kev and 63.9 kev were 
found to be 0.767+ 0.060 and 0.456+ 0.040 barn, 
respectively. 


3.14-hr ISOMERIC LEVEL IN Y™. R. L. Heath, 

J. E. Cline, C. W. Reich, E. C. Yates, and E. H. 
Turk, National Reactor Testing Station, Phillips 
Petroleum Company, Idaho Falls, Idaho (Received 
March 24, 1961; revised manuscript received 
May 5, 1961). 


A previously unreported 3.14-hr isomeric ac- 
tivity has been produced in Y” by thermal neutron 
irradiation of Y®*. The thermal capture cross 
section for the production of the isomer was meas- 
ured to be 1.0+0.2 millibarn. The decay of the 
isomeric level is characterized by the emission 
of two cascade gamma rays of 0.482 and 0.203 
Mev. Conversion coefficient measurements in- 
dicated that the isomeric level at 0.685 Mev de- 
cays by M4 or E5 radiation to an intermediate 
level at 0.203 Mev followed by a predominantly 
M1 transition to the ground state of Y*®. Gamma- 
gamma directional correlation measurements and 
internal conversion measurements indicated that 
the level at 0.203 Mev has spin 3 with odd parity 
and that the level at 0.685 Mev has spin 7 or 8 
with even parity. A shell-model configuration of 
(Zy24,2) has been assigned to the isomeric state. 


(n,d) and (2,p) REACTIONS NEAR Z =50. R. A. 
Peck, Jr., Brown University, Providence, 
Rhode Island (Received March 10, 1961). 


An emulsion study is reported of charged par- 
ticles produced by 14-Mev neutron bombardment 
of Rh’, In'45, Sn4®, Sn’°, Sb, and Te. For all 
but Te (no detectable yield) cross sections and 
spectra are presented, with distributions over 
the first 40° of laboratory angle of energy groups 
from Rh, In, and Sb. Contrary to an assumption 
common in earlier work, there is strong evidence 
that the (z,d) reaction contributes strongly. Five 
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peaks among the Rh, In, and Sb spectra are 
identified with pickup transitions to known re- 
sidual levels, four of their angular distributions 
fitting the predicted Butler forms and the fifth 
providing the new assignment 0-, 1- for the 1.88- 
Mev level in Ru'®. The wide (z,np) group is 
found at the expected energy in the Rh, In, Sn"*, 
and Sb spectra; its angular distribution is anoma- 
lous for Rh but displays the expected isotropy in 
the other three cases. Up to at least 6-Mev ex- 
citation the (2,p) gross structure is dominated by 
single-particle effects, the uncontaminated (n, p) 
yield obeying predictions of the Nilsson model as 
to spectral concentration and angular distribution; 
the low collective levels excited in (p,p’) are not 
observed. Systematic behavior of the direct-in- 
teraction radius for (n,d) and (n,p) and of the re- 
duced width for pickup are found to be reasonable. 
It is inferred that the parent state for proton 
pickup with low residual excitation is almost 
purely a single-particle state in the case of Sb, 
has a strong single-particle character in Rh, and 
a very weak one in In. 


PHOTOALPHA REACTION IN Sb™?. John H. 
Wolfe and John P. Hummel, Department of 
Chemistry and Chemical Engineering and De- 
partment of Physics, University of Illinois, 
Urbana, Illinois (Received March 20, 1961). 


The Sb!**(y, @)In"’-""" reaction has been stud- 
ied by determining the radioactivities of the prod- 
uct nuclei in samples that were irradiated with 
bremsstrahlung of maximum energy varying from 
15.5 to 24 Mev. Excitation functions for the total 
cross section and the cross sections to each iso- 
mer were obtained. The total cross section rises 
steadily over the energy range studied, reaching 
a value of 360 microbarns at 24 Mev. The ratio 
of the cross section for the direct production of 
the ground state to the cross section for the pro- 
duction of the isomer is constant over the energy 
range studied at a value of 2.60+0.40. The total 
excitation function up to 18 Mev agrees well with 
cross sections calculated on the basis of a sta- 
tistical theory for compound-nucleus decay. The 
cross-section predictions could not be made for 
higher energies. The observed ground to isomer 
cross-section ratio is consistent with that ex- 
pected from a compound-nucleus mechanism and 
probably inconsistent with that expected from a 
direct mechanism. Thus, the reaction appears 
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to involve compound-nucleus processes for the 
most part. 


SHORT -LIVED ISOMERS OF Ge”, As™, Br”, 
AND Tc. A. W. Schardt* and Albert Goodman, f 
Los Alamos Scientific Laboratory, University of 
California, Los Alamos, New Mexico (Received 
March 20, 1961). 


Short-lived isomers have been produced with 
the pulsed beam of a Van de Graaff generator and 
observed between pulses with scintillation detec- 
tors. The results are summarized in the follow- 
ing table: : 





Isomeric Observed 
transition gamma ray 
Isotope Half-life (kev) (kev) 
Ge™ 20.340.3 msec 23(M2) 175 
As" 8.0+0.3 sec 283(M3) 283+5 
Br® 118.0+1.5 usec 149(M2) 149+2, 32+2 
Tc 8.15 +0.20 usec 177 +4 
Te 15.5 +0.8 usec 43 +3 





*Present address: Advanced Research Project Agency, 
Washington, D. C. 

*+Present address: Sandia Corporation, Albuquerque, 
New Mexico. 


ELECTRON-SCATTERING STUDY OF NUCLEAR 
LEVELS IN COBALT, NICKEL, LEAD, AND 
BISMUTH. H. Crannell, R. Helm, H. Kendall, 

J. Oeser, and M. Yearian,* High-Energy Physics 
Laboratory, Stanford University, Stanford, Cali- 
fornia (Received March 24, 1961). 


We have observed inelastic scattering of 183- 
Mev electrons through angles of 40°-90° in the 
laboratory, leading to excitation of discrete nu- 
clear excited states in Ni**, Co, Ni®, Pb?%, 
and Bi?°*. The excitation energies were below 
8 Mev. Born-approximation analysis of the 
measured inelastic form factors was used to 
deduce the multipolarities ; (when not previously 
known), and, by extrapolation, the transition 
rates, for 15 corresponding gamma transitions. 
A number of groups of electric transitions for 
\=2, 3, and 4 were observed, each group having 
strikingly similar form factors. In all but one of 
these groups the ratios G of the observed gamma 
transition rates to the single-particle predictions 





were greater than 15, and for some transitions 
from 30 to over 100. One of the groups, in co- 
balt and the nickels, contains the 1.33-Mev E2 
transition to the first excited state of Ni®. An- 
other group consists of fast E3 transitions, seen 
in all five niclei, from states known as the 
“anomalous levels.” They included the transition 
to the first excited state in Pb?°* (G=31) anda 
transition in Bi?” identical in energy and form 
factor. Among three slow £4 transitions in co- 
balt and the nickels was the 2.50-Mev 4+ —0* 
transition in Ni®. The E4, E3, and an E2 tran- 
sition in Co™ identify states analogous to the 4", 
3”, and 2* seen in the neighboring even-even 
nuclei. The last two transitions are strongly 
enhanced. A pair of fast 4.30-Mev £4 transi- 
tions was observed in Pb?” and Bi®”; their speed 
(G = 37) indicates that they may constitute the 
lowest energy configuration of 16-pole mode of 
excitation of the nuclear surface. Values of the 
collective vibrational parameters C, and B, 

and the degree to which some of the transitions 
exhaust ordinary sum rules support the con- 
clusion that the inelastic scattering process 

is strongly exciting nuclear collective excita- 
tions. Some of the observed results are ex- 
pected on the basis of the theory of collective 
vibrational excited states; some are the con- 
sequence of unidentified configurations. 


ISOMER RATIO FOR THE Sn""*(y,p) REACTION. 
John P. Hummel, Department of Chemistry and 
Chemical Engineering and Department of Physics, 
University of Illinois, Urbana, Illinois (Received 
March 29, 1961). 


The relative yields of In” and In” from the 
Sn'4%(y,p) reaction have been measured by count - 
ing the gamma rays emitted by a sample of SnO, 
enriched in Sn™* that had been irradiated with 
24-Mev bremsstrahlung. The measured ground 
to isomer yield ratio was 0.65+0.15, and is con- 
sistent with predictions made on the basis of both 
the Wilkinson direct mechanism and the compound 
nucleus mechanism. However, the values of the 
isomer ratios for the series of even-even tin iso- 
topes (Sn'"*, Sn’#°, Sn’*?, and Sn’, the latter 
three from the work of Yuta and Morinaga) are 
not at all consistent with the predictions of the 
compound nucleus mechanism but can be ration- 
alized on the basis of the Wilkinson model, lend- 
ing support for that theory for direct photoproton 
reactions. 
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ENERGY LEVELS IN Na”'. R. E. Benenson” and 
L. J. Lidofsky, Columbia University, New York, 
New York (Received March 17, 1961). 


The Na” nucleus has been studied through two 
reactions: the Ne?°(d,n)Na** reaction using a fast- 
neutron spectrometer, and the Ne?°(p,7)Na”! re- 
action using both single-crystal and Hoogenboom 
techniques. E , for the neutron spectroscopy for 
the major part of the work was centered about 
4.87 Mev, while the capturing resonance for the 
gamma-ray work was at Ep =1.17 Mev. Enriched 
neon was employed as target gas. Energy levels 
in Na”! are found at excitations of 0.334 0.03 Mev, 
1.77+ 0.05 Mev, 2.42+0.04 Mev, 2.80+ 0.06 Mev, 
and 3.61+0.06 Mev. From the neutron work level 
parameters can be assigned to the 0.33- and 2.42- 
Mev levels, while from the gamma-ray work spin 
limits and parity may be assigned to the 3.61- Mev 
level. This latter level corresponds to the 1.17- 
Mev capture resonance, and a study of the de- 
excitation cascades has been partially completed. 
A discussion of the level scheme of Na** in terms 
of the collective model is given. 


*Permanent address: City College, New York, New 
York. 


ANTISHIE LDING OF NUCLEAR ELECTRIC 
HEXADECAPOLE MOMENTS. R. M. Stern- 
heimer, Brookhaven National Laboratory, Upton, 
New York (Received March 27, 1961). 


The antishielding factor 7. for a possible nu- 
clear electric hexadecapole moment has been 
calculated for the Cu*, Ag*, and Hg** ions, 
using the Hartree-Fock wave functions for the 
3d, 4d, and 5d electrons involved. It was found 
that n..(Cu*) =- 1200, n,.(Ag*) = - 8050, and 
n.AHg*t*)=-63 000. The implication of these re- 
sults is discussed. 


CONFIGURATION MIXING AND THE EFFECTS 
OF DISTRIBUTED NUCLEAR MAGNETIZATION 
ON HYPERFINE STRUCTURE IN ODD-A NUCLEI. 
H. H. Stroke and R. J. Blin-Stoyle,* Physics De- 
partment, Laboratory for Nuclear Science and 
Research Laboratory of Electronics, Massachu- 
setts Institute of Technology, Cambridge, Mas- 
sachusetts, and V. Jaccarino, Bell Telephone 
Laboratories, Murray Hill, New Jersey (Received 
January 16, 1961). 


The theory of Blin-Stoyle and of Arima and 
Horie, in which the deviations of the nuclear 
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magnetic moments from the single-particle 
model Schmidt limits are ascribed to config- 
uration mixing, is used as a model to account 
quantitatively for the effects of the distribution 

of nuclear magnetization on hyperfine structure 
(Bohr -Weisskopf effect). A diffuse nuclear charge 
distribution, as approximated by the trapezoidal 
Hofstadter model, is used to calculate the re- 
quired radial electron wave functions. A table 

of single-particle matrix elements of R? and R‘* 

in a Saxon-Woods type of potential well is in- 
cluded. Explicit formulas are derived to permit 
comparison with experiment. For all of the avail- 
able data satisfactory agreement is found. The 
possibility of using hyperfine structure measure- 
ments sensitive to the distribution of nuclear 
magnetization in a semiphenomenological treat - 
ment, in order to obtain information on nuclear 
configurations, is indicated. 


*On sabbatical leave from the Clarendon Laboratory, 
Oxford, England. 


ISOTOPIC-SPIN SELECTION RULE VIOLATION 
IN THE B?°(d, a)Be® REACTION. J. R. Erskine* 
and C. P. Browne, Department of Physics, Uni- 
versity of Notre Dame, Notre Dame, Indiana 
(Received March 28, 1961). 


All previous tests of the isotopic-spin selection 
rule in (d, a) reactions were obscured by statis- 
tical weight factors because in each case the ini- 
tial and final nuclear states had spin and parity 0*. 
The B’°(d, a)Be® reaction provides a test of the 
isotopic-spin selection rule free from this re- 
striction. The energy levels of Be® near the low- 
est T=1 level were studied with the Li®(He’, p)Be* 
and Be*(He*, a)Be® reactions as well as with the 
B’°(d, a)Be® reaction. Energy levels in Be® were 
found at 16.623+ 0.010 Mev, 16.921+0.010 Mev, 
and 17.637+ 0.006 Mev. The widths are 95+ 20 
kev, 85+ 20 kev, and <15 kev, respectively. The 
first of these levels is the lowest T=1 state where- 
as the second is 7 =0 and the third probably T =1. 
Energy levels at 16.08 Mev and a J=2 level at 
17.7 Mev, reported by other laboratories, were 
not observed. The ratio of the differential cross 
sections for formation of the 16.62- and 16.92- 
Mev levels was measured over a range of angles 
and bombarding energies. The ratio is about 1.4 
and is roughly constant for both the Li®(He’, p)Be’ 
and B’°(d, a)Be® reactions. This implies complete 
violation of the selection rule because the latter 
reaction should not go to the T=1 level. Argu- 
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ments are given which indicate that the T=0 im- 
purities in the 16.62-Mev T =1 level are no more 
than a few percent. Consequently, the failure of 
the selection rule probably results from the com- 
plete intermixing of T=0 and T=1 states in the 
Cc’ compound nucleus near 28-Mev excitation. 
Groups from the C**(He*, p)N™ reaction were seen 
corresponding to levels in N“ at 5.691 + 0.008, 


5.834+ 0.008, 6.203+ 0.008, and 6.440+ 0.008 Mev. 


*Present address: Laboratory for Nuclear Science, 
Massachusetts Institute of Technology, Cambridge, 
Massachusetts. 


PROTONS FROM ALPHA-INDUCED REACTIONS. 
Wayne Swenson and Nikola Cindro,* Department 
of Physics and Laboratory for Nuclear Science, 
Massachusetts Institute of Technology, Cam- 
bridge, Massachusetts (Received December 21, 
1960). 


The results of proton energy spectra meas- 
ured at several angles from 30.5-Mev alpha- 
particle-induced reactions on Al*’, V™, Co*, 
As”, Nb®, Rh’™, In'®, and Ta’®’ were analyzed 
using the statistical model. The analysis yielded 
the differential cross section d*o/dQdE and the 
relative level density w(E) of the residual nucleus 
as a function of proton and excitation energy of 
the residual nucleus. The nuclear temperature 
1/T =d(lnw)dE and the level density parameter 
aof w=C exp[(aE)”] were obtained. The energy 
and angular dependence of the spectra are ad- 
equately described by the statistical model at 
back angles, with the indication of the presence 
of a direct reaction mechanism contribution at 
forward angles, which extends to high excitation 
energies. 


*On leave of absence from the Institute “Rudjer 
Boskovic,” Zagreb, Yugoslovia. 


POLARIZATION CORRELATION MEASURE- 
MENTS ON Eu’. C. V. K. Baba, Atomic Energy 
Establishment, Trombay, Bombay, India, and 

S. K. Bhattacherjee, Tata Institute of Fundamen- 
tal Research, Bombay, India (Received March 
15, 1961). 


The spins and parities of the 1400-kev and 
1723-kev levels in Gd™™ are discussed on the 
basis of the polarization-directional correlation 


measurements on the 1277-123 kev and 725-998 
kev gamma-gamma cascades following the 6 de- 
cay of Eu’™, In the case of the 1277-123 kev cas- 
cade, the directional correlation measurements 
and the polarization correlation measurements 

of the 1277-kev gamma ray are made with liquid 
and solid sources. An attenuation in the aniso- 
tropy of directional correlation and in the degree 
of linear polarization correlation has been ob- 
served in both the sources, the attenuation in the 
solid source being more. In the 725-998 kev cas- 
cade, the polarization of either of the two radia- 
tions has been measured in two separate experi- 
ments. On the basis of these measuremeits, the 
spins and parities of both the 1400-kev and the 
1723-kev levels in Gd™™ are assigned as 2~. The 
performance of the apparatus has been checked 
by measuring the polarization correlations of the 
known cascades in Ni®, Ti*®, Pd', and Sr®. 


THEORY OF AVERAGE NEUTRON REACTION 
CROSS SECTIONS IN THE RESONANCE REGION. 
P. A. Moldauer, Argonne National Laboratory, 
Argonne, Illinois (Received February 2, 1961). 


The scattering matrix for compound-nucleus 
processes is studied in the R-matrix formalism, 
using a series expansion which is due to Thomas. 
It is shown that this series generally converges 
when (a) the average total resonance width is less 
than the average resonance spacing, (b) the num- 
ber of important channels is not too large, and 
(c) the width amplitudes have random signs. The 
treatment also suggests strongly that the series 
does not converge in the continuum region. In the 
region of convergence the exact relationship be- 
tween the channel transmission factor T; and the 
ratio of partial width to level spacing is found in 
the absence of direct-scattering reactions to be 
Te =2n(Tyc)/D - PI yc)*/D*. The quadratic term 
is shown to be important in the vicinity of optical- 
model maxima. Correction terms to the Hauser- 
Feshbach relations for average reaction cross 
sections arising from the higher order terms of 
the series are obtained and are found to depend 
on the statistical properties of both resonance 
widths and resonance spacings. The effect on 
average neutron inelastic, compound elastic, 
and capture cross sections is discussed and an 
example of a calculation is presented. 
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RESONANCE MODEL OF A-K° PRODUCTION. 
Akira Kanazawa,* Department of Physics, Purdue 
University, Lafayette, Indiana (Received March 
30, 1961). 


A resonance model is proposed to explain the 
excitation function, the angular distribution, and 
the large polarization of A in the reaction p+77 
~A+K°. It is assumed that there exists a low 
angular momentum resonance in the channel p+z7 
~ (resonant state)—A+K°. There are five real 
parameters in this model. Two of these are the 
coupling constants of the usual interactions. The 
other three are the position, half-width, and 
height of the assumed resonance. With reason- 
able choices of parameters a fairly good fit is 
obtained, for both a scalar and a pseudoscalar 
K meson, to the experimental data in the interval 
910 to 1300 Mev of the pion kinetic energy in the 
laboratory system. 


*On leave from Hokkaido University, Sapporo, Japan. 


PHOTOPRODUCTION OF K* MESONS FROM 
DEUTERIUM. R. L. Anderson, F. Turkot, and 
W. M. Woodward, Cornell University, Ithaca, 
New York (Received February 13, 1961). 


The differential cross section for the reaction 
y+n-+K*+2~ has been measured. The data are 
obtained by observing the charged K particle with 
a magnet spectrometer alternately from hydrogen 
and from deuterium. For y-ray energy of 1122 
Mev and for 0¢, m, = 82°, we get do, -/dops0 
=1.6+0.7. 


PION PRODUCTION AND THE SECOND PION- 
NUCLEON RESONANCE. Charles J. Goebel* 
and Howard J. Schnitzer, Department of Physics 
and Astronomy, University of Rochester, Roch- 
ester, New York (Received March 16, 1961). 


A model for the reaction 7+N ~27+WN at low 
energies, which includes pion-pion interaction 
and final-state interactions in the (3,3) state, 
is discussed. The theory involves 2 parameters 
which are related to the S-wave and P-wave n-7 
scattering lengths. These parameters are chosen 
from a fit to the total cross section for 1~+p— 
Meson production is predicted to be 
primarily in the T=1/2 state. Predictions are 
made for the total cross sections of the various 
channels (e.g., 1*+p+a*++n*t+n, etc.) in the 


Tt ++ +n. 
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energy range from threshold to = 500 Mev, in 
good agreement with experiments. Angular dis- 
tributions are predicted. These are in qualita- 
tive agreement with the 7* angular distribution 
for 7~+p+a~+7*+n. From these data it is sug- 
gested that the S-wave m-m scattering length has 
opposite sign to the P-wave scattering length. 

A conjecture concerning rapidly rising inelastic 
cross sections in a single partial wave is made 
to connect the large T=1/2, D,. production 
cross sections with the T=1/2, D,,. pion-nucleon 
resonance. The 1-7 scattering length found are 
a,=-0.290u", a, =0.122u"", and a,4,=2/5 by 
hypothesis. 


*On leave of absence 1960-61 at the Institute for 
Advanced Study, Princeton, New Jersey. 


ANNIHILATION PROCESS OF NEUTRINO PRO- 
DUCTION IN STARS. H. Y. Chiu, Institute for 
Advanced Study, Princeton, New Jersey (Re- 
ceived March 31, 1961). 


The rate at which the energy of a black-body 
radiation is converted into neutrinos by the pair 
annihilation process e~+e* + v+7 has been cal- 
culated. At T~6x10°°K, the relaxation time 
for such conversion process is around 100 sec- 
onds for pure radiation. Since neutrinos have 
very long mean free path (> stellar dimensions) 
they will escape, thus carrying away the energy. 
This process therefore will be of astrophysical 
importance. The rate of energy loss dU/dt is 
tabulated, as a function of temperature and den- 
sity, together with the chemical potential, the 
pressure, and the electron-positron energy. 
This table should be useful for numerical in- 
tegrations of stellar structure equations in the 
temperature range (0.5—10)x10°°K, and the den- 
sity range 0—10° gram/cc. 


MEAN LIFETIME OF THE NEUTRAL PION. 

R. G. Glasser, N. Seeman, and B. Stiller, Nu- 
cleonics Division, Naval Research Laboratory, 
Washington, D. C. (Received March 30, 1961). 


An estimate of the mean lifetime of the 7° 
meson has been obtained from an experiment 
employing a direct time-of-flight technique 
first attempted by Harris et al. in 1957. This 
method is based upon the observation in nuclear 
emulsion of the decay of the K,9* meson (Kt + 
1*+7°) and the subsequent decay of the 1° via 
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the Dalitz mode, 1°+e*+e~+y. In the present 
experiment we were able to utilize a new fine- 
grained emulsion (Ilford L.4) that yielded mark- 
edly improved resolution. The availability of 
the separated K* beam from the bevatron at 
Berkeley permitted detection and measurement 
of 76 Dalitz decays. We obtain for the mean 
lifetime of the 7°, 7 =(1.9+0.5) x107** sec. 


NUCLEON POLARIZABILITY CORRECTION TO 
HIGH-ENERGY ELECTRON-NUCLEON SCAT- 
TERING. N. R. Werthamer,* University of Cal- 
ifornia, Berkeley, California, and M. A. Ruder- 
man,t New York University, New York, New 
York (Received March 23, 1961). 


The contribution of nucleon polarizability to 
ultrarelativistic electron-nucleon scattering 
cross sections is estimated and found to be 
small for nonforward scattering angles at all 
energies. 


*Present address: University of California, La 
Jolla, California. 

+tOn leave from the University of California, Berkeley, 
California. 


BROKEN SYMMETRIES AND BARE COUPLING 
CONSTANTS. Murray Gell-Mann and Fredrik 
Zachariasen, California Institute of Technology, 
Pasadena, California (Received March 20, 1961). 


There are known cases of symmetry laws valid 
for one kind of interaction but broken by another. 
Each symmetry is then supposed to be exact for 
bare masses and coupling constants but only ap- 
proximate for the renormalized quantities, like 
neutron and proton masses. We ask how the equal- 
ity of unrenormalized constants can be rephrased 
as a statement about measurable quantities. This 
question is particularly important in connection 
with proposed strong-interaction symmetries that 
are supposed to be badly broken. The answer ap- 
pears to involve the limits of ratios of experi- 
mental quantities at very high momenta. We dis- 
cuss first the connection between wave function 
renormalizations and weak and electromagnetic 
form factors. Then we take up the measurement 
of strong-interaction vertex renormalization fac- 
tors by the study of scattering amplitudes at en- 
ergies and momentum transfers large compared 
to all masses. The last part of the work is based 


in part on indications from the perturbation de- 
velopment of pseudoscalar meson theory, but we 
hope it will point the way to similar results ina 
better theory. 


NOTE ON REARRANGEMENT COLLISIONS. 
T. B. Day, L. S. Rodberg, G. A. Snow, and 

J. Sucher, Physics Department, University of 
Maryland, College Park, Maryland (Received 
March 29, 1961). 


The conventional Born approximation formula 
for rearrangement collisions is used extensively 
in both atomic and nuclear physics. This formula 
contains a direct contribution from the heavy 
particle or “core” interaction. A straightforward 
demonstration shows that for the usual case of 
a massive core this contribution does not appear, 
so that the only effect of this interaction is to 
distort the incident and outgoing waves. Such 
problems as the “post-prior” discrepancy are 
clarified. 


ASYMPTOTIC BEHAVIOR AND SUBTRACTIONS 
IN THE MANDELSTAM REPRESENTATION. 
Marcel Froissart,* Department of Physics, Uni- 
versity of California, Berkeley, California (Re- 
ceived March 27, 1961). 


It is proved that a two-body reaction amplitude 
involving scalar particles and satisfying Mandel- 
stam’s representation is bounded by expressions 
of the form Cs In’s at the forward and backward 
angles, and Cs™ 1n™s at any other fixed angle in 
the physical region, C being a constant, s being 
the total squared c.m. energy. This corresponds 
to cross sections increasing at most like In’s. 
These restrictions limit the freedom of choice of 
the subtraction terms to six arbitrary single spec- 
tral functions and one subtraction constant. 


*On leave of absence from Centre d’ Etudes Nucléaires 
de Saclay, Gif sur Yvette, Seine et Oise, France. 


CORRELATIONS BETWEEN CUSPS IN DIFFER- 
ENTIAL CROSS SECTIONS AND IN POLARIZA- 
TIONS. M. Nauenberg and A. Pais, Institute for 
Advanced Study, Princeton, New Jersey (Received 
March 31, 1961). 


Experiments have shown that partial waves with 
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l> 1 appear in 7~+p~A+K° at the DK thresholds. 
This necessitates a reconsideration of the criteria 
sufficient to determine the A parity P(ZA) by the 
method of cusps. In this paper we start from the 
usual assumption that contributions which show 

a cusp in either the differential cross section or 
the polarization may be ignored beyond some op- 
timal power in cos@. On this sole basis, previ- 
ously stated criteria are rendered inadequate due 
to the occurrence of Minami and other ambigui- 
ties. It is shown that under suitable circumstan- 
ces there exist unambiguous correlations between 
certain properties of cross-section cusps and of 
polarization cusps. These correlations could pos- 
sibly be of use to determine P(ZA) and give infor- 
mation as to which states contribute significantly 
to the AK production at ~900 Mev. The finite sep- 
aration between D°K° and ©“K* thresholds is taken 
into account. The results are summarized ina 
table of cusp properties. 
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ANALYTIC PROPERTIES AND RESCATTERING 
CORRECTION TO THE BORN APPROXIMATION 
FOR TRANSITION MATRIX ELEMENTS. B. 
Bosco,* Institute of Theoretical Physics, Depart- 
ment of Physics, Stanford University, Stanford, 
California (Received March 22, 1961). 


In this paper the analytic properties of a matrix 
element of a general operator between a bound 
state and a scattering state are studied in the 
framework of the Schrédinger theory. It is shown 
that the singularities of such a matrix element 
are easily inferred from those of the Born ap- 
proximation. Finally, using the fact that the 
possible singularities which are not contained 
in the Born approximation are located far apart 
from those included in the lowest approximation, 
a simple formula is derived which allows one to 
obtain the rescattering correction to the Born 
approximation using the phase shifts explicitly. 

*Permanent address: Istituto di Fisica dell’ Univer- 
sita, Torino, Italy. 














